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ABSTRACT 


Resolution of both the extent and mechanism of lateral heterogeneity in the upper 
mantle constrains the nature and scales of mantle convection. Oceanic regions are of 
particular interest as they are likely to provide our closest glimpse at the patterns of 
temperature anomalies and convective flow in the upper mantle because of their young 
age and simple crustal structure relative to continental regions. Our objectives in this 
thesis are to determine lateral variations in the seismic velocity and attenuation structure 
of the lithosphere and asthenosphere beneath the oceans, and to combine these 
seismological observations with the data and theory of geoid and bathymetry anomalies 
in order to test and improve current models for seafloor spreading and mantle 
convection. We concentrate on determining variations in mantle properties on a scale of 
about 1000 km, comparable to the thickness of the upper mantle. Seismic velocity, 
geoid, and bathymetry anomalies are all sensitive to variations in upper mantle density, 
and we formulate inversions to combine quantitatively these different data and search for 
a common origin. V ariations in mantle density can be either of ther mal or compositional 
origin and are presumably related to mantle convection and differentiation. 

By means of a large data base of digital seismograms and waveform cross- 
correlation and spectral ratio techniques, we have measured SS-S differential travel time 
residuals and differential attenuation in order to determine lateral variations in upper 
mantle structure beneath the Mid- Atlantic Ridge and East Pacific Rise. Differential 
travel times of such phases as SS and S with identical source and receiver have the 
advantage that residuals are likely to be dominated by contributions from the upper 
mantle near the surface bounce point of the reflected phase (SS). Under this assumption, 
differential SS-S travel time residuals are mapped at the SS bounce points as a means of 
delineating lateral variations in mantle structure. After removing the signature of 
lithosphere age, we find evidence for long-wavelength variations in SS-S residuals along 
the Mid-Atlantic Ridge. The dominant wavelength of these variations is 1000 to 2000 
km. These travel time anomalies correlate qualitatively with along-axis va ria tions in 
bathymetry and geoid height. We formulate a joint inversion of travel time residual. 
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geoid height, and bathymetry under the assumption that all arise from variations in upper 
mantle temperature or bulk composition (parameterized in terms of Mg#). The inversion 
employs geoid and topography kernels which depend on the mantle viscosity structure. 
Inversion for temperature perturbations alone provides good fits to travel time and geoid 
data. The fit to topography, which is likely dominated by unmodeled crustal thickness 
variations, is not as good. The inversions for temperature favor the presence of a thin 
low viscosity layer in the upper mantle and temperature perturbations concentrated at 
depths less than 300 km. Compositional variations alone are unable to match the travel 
time and geoid or bathymetry data simultaneously. A joint inversion for temperature 
and composition provides good fits to both geoid and travel time anomalies. 

Temperature variations are i 50 K and compositional variations are i 0.5-3 % Mg# for 
models with the temperature variations uniformly distributed over the uppermost 300 km 
and the compositional variations either distributed uniformly over the same interval or 
concentrated at shallower depths. The magnitudes of these variations are consistent with 
the chemistry and geothermometry of dredged peridotites along the Mid-Atlantic Ridge. 

Differential travel times of SS-S pairs in the east central Pacific show several 
differences from the north Atlantic. The most obvious difference is that the travel time 
residuals are significantly larger than in the Atlantic, even at a fixed age. The travel time 
- age relation is weaker in the Pacific, although this may be partially atmbutable to the 
fact that we have not sampled a large range of plate ages in the eastern Pacific. In the 
Atlantic our results are not consistent with the presence of a simple pattern of azimuthal 
anisotropy, while in the Pacific the data are consistent with the presence of weak 
anisotropy in the upper mantle. It has been suggested that anisotropy may be more 
pronounced at fast spreading rates than at slow spreading rates both in the lithosphere 
(due to a rate dependence of the mechanism for orienting olivine crystals in the 
lithosphere) and the asthenosphere (because the asthenospheric flow beneath fast moving 
plates is likely to take the form of a progressive simple shear which can induce a lattice 
preferred orientation of olivine crystals), and our results are consistent with this 
suggestion. There is substantial ambiguity in our anisotropy measurements for the 
Pacific, however, due to a poor sampling of azimuths, so that it is also possible that 
lateral heterogeneity rather than azimuthal anisotropy is producing the observed 
azimuthal pattern. Sampling at a more uniform distribution of azimuths should make this 
result less ambiguous, and as more seismic stations are deployed at new geographic 
locations our chances of resolving this issue will improve. 

Inversion of travel time residuals, geoid, and bathymetry data for the eastern Pacific 
indicates that compositional variations alone are inadequate to match all of the data 
simultaneously, similar to our results for the north Atlantic. Temperature variations 
alone, however, produce significant variance reduction. The inversion solutions indicate 
excess temperature in the vicinity of the Galapagos hotspot in the range 50- 150 K. 
Further analysis is needed to determine the effects of subduction zone structure and 
possible crustal thickening in the eastern Cocos plate region. 

As a complement to the study of travel times, we have measured SS-S differential 
attenuation in the north Atlantic region. Mapping seismic Q in the upper mantle is an 
important tool for assessing mechanisms of lateral heterogeneity because the attenuation 
of seismic waves is sensitive to variations in temperature and to partial melting. 
Differential attenuation is positively correlated with SS-S travel time residual. Both 
differential attenuation and travel time residual decrease with increasing seafloor age. 
The age dependence of SS-S travel time residual can be explained entirely by the cooling 
of the oceanic lithosphere, i.e., contributions from the asthenosphere or from a mantle 
melt fraction are not required. On the assumption that plate cooling also dominates the 
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variation of differential attenuation with age, we derive an empirical Q* 1 -temperature 
relation for the oceanic lithosphere. The variation of Q* 1 with temperature that we derive 
is not as strongly dependent on temperature as that observed in laboratory studies. 
Systematic long-wavelength (1000-6000 km) variations in upper mantle differential 
attenuation are evident along the axis of the Mid-Atlantic Ridge. These variations 
correlate approximately with long-wavelength variations in shear wave travel time 
residuals and are attributed to along-axis differences in upper mantle temperature. 


Thesis Supervisor: Sean C. Solomon 
Title: Professor of Geophysics 
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Chapter 1 


Introduction 

It is important to measure lateral variations in the Earth's mantle because of the 
key role of this information in the formulation and testing of theories of the interior 
dynamics of the Earth. The mantle convection associated with the escape of heat from the 
Earth’s interior gives rise to global tectonic activity, but the exact forms and scales of 
convection and the interaction of convective processes with the lithosphere are still 
controversial. Mid-ocean ridges are of particular interest as they are known to be sites of 
upwelling and elevated temperature. The excess temperatures have expressions in the 
seismic velocity and density fields. The focus of this study is to examine lateral 
variations in upper mantle properties near mid-ocean ridges on a scale of about 1000 km, 
comparable to the thickness of the upper mantle. The broad questions we seek to address 
include the following: What are the variations in mantle convection at this scale? Are 
they observable as anomalies in temperature-sensitive physical properties? Are such 
differences in temperature manifested as differences in melt production and thus in the 
chemistry of the mantle residuum? If so, are these compositional variations observable? 

Learning about the large-scale structure of the Earth has become a 
multidisciplinary effort. Consequently in this thesis we adopt a multidisciplinary 
approach to the problem of resolving lateral variations in mantle properties at the 1000-km 
scale. Specifically we combine seismic data with geoid and bathymetry anomalies. It is 
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important when combining these different data to consider dynamic effects, possible 
compositional variations, and the spatial scales of convection. 

Seismology has made great contributions to the resolution of lateral velocity 
differences in the mantle, and advances in techniques and a steady improvement in the 
capabilities of global seismic networks have greatly increased our knowledge of the 
Earth's interior. No longer can the Earth be treated as a rigid body possessing radial 
symmetry; large lateral variations have been observed in the upper mande as well as in the 
crust from both regional and global studies. Global tomographic models such as M84C 
[Woodhouse and Dziewonski, 1984] utilize long-period surface waves to describe the 
Earth’s upper mande heterogeneity on a global scale, up to spherical harmonic degree and 
order 8 (X= 5000 km). In this thesis we seek to examine lateral variations on a somewhat 
shorter scale (wavelengths of about 1000 to 5000 km). To do this we make use of body 
waves rather than surface waves, and rather than performing a global analysis we 
concentrate on ocean ridge environments. 

The body waves used in this study are direct and surface-reflected long-period S 
waves. Differential travel times and attenuation are measured using SS and S wave pairs 
from the same source and station. Much progress has been made in resolving upper 
mande heterogeneity through differential travel times [Sipkin and Jordan , 1976, 1980a; 
Stark and Forsyth , 1983; Butler , 1979; Kuo etal., 1987; Woodward and Masters, 1991]. 
Among the advantages of differential rather than absolute times is that source and receiver 
effects are approximately common to both phases and are thus largely eliminated by 
differencing. If we assume that the lower mande is relatively homogeneous and that the 
portions of the wave paths in the upper mande are steep, the differential travel time 
anomaly is associated with upper mande heterogeneity centered beneath the surface 
bounce point of the reflected (SS) phase. This technique is thus well suited to the 
investigation of horizontal variations in upper mande structure. 
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Other data sensitive to variations in mantle properties at depth include oceanic 
bathymetry and geoid height Variations in mantle density can be either of thermal or 
compositional origin and, like seismic velocity, are presumably related to mantle 
convection and differentiation. Geoid (or gravity) and topography have become the most 
commonly used tools for mapping out and constraining models of upper mantle 
convection [e.g., Anderson et al., 1973, McKenzie and Bowin, 1976; McKenzie , 1977; 
McKenzie et al., 1980, Parsons and Daly , 1983; Buck and Parmentier, 1986; Craig and 
McKenzie, 1986]. Several workers [Dziewonski et al., 1977; Nakanishi and Anderson, 
1984; Tanimoto and Anderson, 1984; Stark and Forsyth, 1983; Dziewonski, 1984; Kuo 
et al., 1987] have noted correlations of geoid and travel time (or velocity structure) at a 
number of different wavelengths, although only a few [Hager et al., 1985; Hager and 
Clayton, 1989; Hager and Richards, 1989] have combined observational seismology with 
geoid anomalies in a quantitative and dynamically consistent manner. 

Since very different convective flows can produce the same geoid and surface 
topography, the inversion of these data alone for the thermal or compositional source 
function is nonunique. Because this inverse problem is not well posed, most studies 
have concentrated on forward modelling, i.e., varying the parameters of a simple model 
until a good fit to the data is achieved. With this approach, there is no guarantee that the 
set of parameters which give the best fit to the data is unique and that the correct solution 
has been isolated. Including seismic data provides additional constraints which are 
sufficient to allow us to formulate simple one-dimensional inversions. 

In Ch ap ter 2 we describe a comprehensive study of differential travel times in the 
north Atlantic region. We present details of the measurement procedure we use to obtain 
SS-S measurements from seismograms recorded by the GDSN network, as well as other 
aspects of the collection and reduction of these data. We examine the relationship of these 
measurements to such factors as seafloor age and upper mantle anisotropy, and we 
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examine the correlation with geoid and bathymetry. We develop techniques for 
quantitatively combining the observed travel time residuals with geoid and bathymetry 
anomalies in a joint inversion for upper mantle thermal and compositional variations. 
Since the inversion includes dynamic geoid and topography kernels which depend upon a 
given viscosity model, inversions can be carried out with a range of different viscosity 
models in order to determine which are the most consistent with the data. We invert 
separately and jointly for distributions of thermal perturbations and compositional 
variations which best produce the observed travel time, geoid, and bathymetry. 

In Chapter 3 we examine differential attenuation between the phase pairs SS and 
S, with SS bounce points in the north Atlantic, and we examine the possible relations 
between Qs* 1 and the temperature of the upper mantle. Earlier studies with body waves 
suggest that the upper mantle beneath mid-ocean ridges should display a significantly 
greater than average level of S-wave absorption [e.g., Molnar and Oliver, 1969; 

Solomon, 1973]. We examine differential shear wave attenuation in the north Atlantic 
and its variation with lithosphere age, and we utilize thermal models to interpret the 
observations in terms of an empirical relation between differential attenuation and 
temperature. This work is a complement to the study of SS-S differential travel times 
presented in Chapter 2. The data show clear evidence for a decrease in Qs' 1 in the upper 
mantle with increasing plate age. There are also systematic along-axis variations in 
differential attenuation, and we address whether these variations might be the result of 
along-axis differences in the characteristic temperature of the upper mantle, such as those 
obtained from the inversions of travel time residuals, geoid, and bathymetry described in 
Chapter 2. 

In Chapter 4 we examine SS-S differential travel times, geoid, and bathymetry in 
the region of the East Pacific Rise in the east-central Pacific. Comparision of the results 
from the East Pacific Rise with those from the Mid- Atlantic Ridge provides us with the 
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opportunity to compare and constrast upper mantle properties beneath fast and slow 
spreading ridges. Several authors [e.g., Duschenes and Solomon, 1977; Stark and 
Forsyth, 1983; Zhang and Tanimoto, 1990b] have found that shear wave velocities in the 
eastern Pacific are significantly lower than those observed in either the Atlantic or Indian 
Oceans. In addition, the presence of anisotropy has been found to be more pronounced 
in the Pacific upper mantle than in the Atlantic upper mantle [Montagner and Tanimoto, 
1990], consistent with predictions of models of shear-induced alignment of olivine 

crystals [McKenzie, 1979; Ribe, 1989]. 

In Chapter 5 we discuss some conclusions of this thesis relating to the 
mechanisms of heterogeneity in the oceanic upper mantle. We raise several questions 
generated by the analyses in this work. Finally, we suggest directions for future research 
designed to help answer these questions. We feel that further application of the 
techniques developed in this thesis will allow constraints to be placed on upper mantle 
viscosity structure, the presence of partial melt, and the mechanisms of lateral 
heterogeneity on both global and regional scales. These techniques will allow the theories 
of mantle convection to be tested and sharpened with abundant and diverse data. 



Chapter 2 


Joint Inversion of Shear Wave Travel Time Residuals and Geoid 
and Depth Anomalies for Long-Wavelength Variations in Upper 
Mantle Temperature and Composition along the Mid-Atlantic 

Ridge 


Introduction 


Seismic velocity and density of upper mantle material are expected to be functions 
of temperature and composition. The delineation of long wavelength variations in these 
physical properties thus provide important constraints on mantle convection, crust-mantle 
differentiation, and mantle chemical heterogeneity. In this study we determine lateral 
variations in upper mantle temperature and composition along the Mid- Atlantic Ridge 
through the combined inversion of shear wave differential travel times, geoid height, and 
bathymetric depth anomalies. 

The advent of seismic tomography has led to a number of three-dimensional maps 
of lateral variations in seismic velocity in the upper mantle, and several such models of the 
north Atlantic region have been developed, both as parts of global studies [e.g.. 

Wood ho use and Dziewonski, 1984; Nakanishi and Anderson, 1984; Tanimoto, 1990] and 
through regional investigations of long-period surface waves [e.g., Honda and Tanimoto, 
1987; Mocquet et ai, 1989; Mocquet and Romanowicz, 1990]. With surface wave 
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methods each wave samples the average vertical variation in upper mantle structure along 
its path, but because of the long wavelengths involved the inversion of phase or group 
velocity from many paths tends to smooth out lateral variations. Body wave travel times 
can provide independent information about upper mantle heterogeneity at potentially shorter 
horizontal scales than surface waves can resolve, and progress has been made in the 
determination of lateral heterogeneity in the north Atlantic through the use of both 
differential and absolute travel times of body waves [Kuo et al., 1987; Grand, 1987, 

1989]. 

The travel times used in this study are differential times of the body wave phase 
pair SS-S (Figure 2. 1). Differential travel times of shear wave pairs are well suited to the 
study of upper mantle heterogeneity [Sipkin and Jordan, 1976, 1980a; Stark and Forsyth, 
1983; Butler, 1979; Kuo et al., 1987; Woodward and Masters, 1991] and have the 
advantage that source and receiver effects are approximately common to both phases and 
are thus largely eliminated by differencing. Under the assumption that the lower mantle is 
relatively homogeneous and that the portions of the wave paths in the upper mantle are 
steep, the differential travel time anomaly can be associated with upper mantle structure 
within a small volume centered beneath the surface bounce point of the reflected (SS) 
phase. This technique is thus well suited to the investigation of horizontal variations in 
structure, but the resolution of variations with depth is poor. 

Oceanic bathymetry and geoid height data are sensitive to variations in mantle 
density at depth. Such variations can be either thermal or compositional in origin and, like 
seismic velocity, are presumably related to mantle convection and differentiation. Geoid 
(or gravity) and topography have become the most commonly used tools for mapping out 
and constraining models of upper mantle convection [e.g., Anderson et al., 1973, 
McKenzie and Bowin, 1976; McKenzie, 1977; McKenzie et al., 1980, Parsons and Daly, 
1983; Buck and Parmentier, 1986; Craig and McKenzie, 1986]. In addition, measurement 
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of the admittance (the spectral ratio of geoid to topography) has been widely utilized to 
estimate the depth and mode of compensation of oceanic swells and plateaus [e.g.. Watts et 
al., 1985, Cazenave et al., 1988; Sandwell and MacKenzie, 1989; Sheehan and McNutt, 
1989]. Several workers [Dziewonski et al ., 1977; Nakanishi and Anderson, 1984; 
Tanimoto and Anderson, 1984; Stark and Forsyth , 1983; Dziewonski, 1984; Kuo et al., 
1987] have noted correlations of geoid and travel time (or velocity structure) at a number of 
different wavelengths, although only a few [Hager et al., 1985; Hager and Clayton, 1989; 
Hager and Richards, 1989] have combined observational seismology with geoid anomalies 
in a quantitative and dynamically consistent manner. 

In this study we present the first formal inversion of geoid, depth, and travel time 
anomaly data for lateral variations in upper mantle temperature and composition along the 
Mid-Atlantic Ridge. Given a distribution of temperature or density perturbations in the 
upper mantle, the forward problem of calculating differential travel time, geoid, and depth 
anomalies is straightforward. This forward problem forms the basis for a joint linear 
inversion of these three types of observations under the assumption that all arise from 
parameterized long-wavelength variations in upper mantle temperature or composition. 
Results of a set of inversions carried out under different assumptions regarding the depth 
extent of lateral heterogeneity and the mantle viscosity structure are compared with other 
constraints on variations in mantle temperature and degree of melt removal. 

MEASUREMENT OF DIFFERENTIAL TRAVEL TIMES 

The seismic data used in this study consist of long-period S and SS phases 
obtained from the Global Digital Seismic Network (GDSN) [Peterson et al., 1976; 

Peterson and Hutt, 1982]; the Network of Autonomously Recording Seismographs 
(NARS), a linear broadband array in western Europe [Nolet and Vlaar, 1982]; and several 
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broadband stations from the global GEOSCOPE network [Romanowicz et al., 1984]. A 
list of stations used in this study is presented in Table 2.1. We use only transversely 
polarized (SH) seismograms (rotated from N-S and E-W components) to avoid interference 
from the SKS phase and contamination from P-SV conversions at the base of the crust and 
other near-surface discontinuities. Recent work by Gee and Jordan [1989] suggests that 
travel times depend on the frequency band used in the analysis. In order to maintain a self- 
consistent data set for our study, additional processing is applied to data from the NARS 
and GEOSCOPE arrays in order to mimic the instrument response of the longer period 
GDSN stations. This processing allows us to measure travel times from a set of 
seismograms that all have essentially the same frequency response. Data from the NARS 
and GEOSCOPE arrays are decimated (with a low-pass antialiasing filter) to a common 
sampling interval of 1 s. The data are further filtered using a noncausal 3-point 
Butterworth filter [Rader and Gold , 1967] with a frequency bandpass of 0.01 - 0.20 Hz. 
This additional filtering greatly improves the signal-to-noise ratio of the SS phase. 

A waveform cross-correlation method is utilized to determine the differential travel 
time between the phases S and SS [Butler, 1979; Stark and Forsyth, 1983, Kuo et al., 
1987]. The procedure involves the construction of a "synthetic” SS pulse from S and the 
evaluation of the cross-correlation function between the real and synthetic windowed SS 
phases (Figure 2.2). The synthetic SS pulse is created from S in the following manner. 

The S pulse is windowed and attenuated (with attenuation parameter t*= 3 s) [Grand and 
Helmberger, 1984; Kuo et al., 1987] to account for the additional time SS travels in the 
man tle, and then a jc/2 phase shift (Hilbert transform) is applied to the attenuated S pulse to 
simulate the frequency-dependent phase shift which the SS wave undergoes at an internal 
caustic [Choy and Richards, 1975]. The differential time is obtained from the peak of the 
cross correlation between the synthetic SS constructed from the S wave and the real SS. 
The residual SS-S times are obtained by subtracting the observed differential time from that 
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predicted by the PREM Earth model [Dziewonski and Anderson, 1981] and correcting for 
Earth ellipticity [Dziewonski and Gilbert, 1976] and SS bounce point bathymetry. Our 
convention is that negative residuals are indicative of either early SS or late S. 

Constant window lengths of 120 s are used for both the S and SS phases. In 
general, the observed differential travel times vary by as much as 1 s depending on how S 
and SS are windowed. Our modelling with synthetic seismograms indicates that 
emphasizing the onset of the SS waveform can lead to bias for bounce points in areas of 
oceanic sediments. The effect of sediments at long periods is to produce precursory 
arrivals from reflections at the base of the sediments and late arrivals from waves which 
travel through the low-velocity sediments and are reflected at the crust-water interface. The 
net effect, after convolving the crustal response with the long-period GDSN instrument 
response, is that the time center of the SS phase is effectively unchanged but the pulse is 
broadened both at the front and at the back. In our procedure the use of a constant window 
containing the entire SS pulse should yield differential travel times that are little affected by 
the presence of sediments. 


DATA 

The north Atlantic is an ideal area for conducting a differential travel time study in 
terms erf the geographic distribution of available events and stations at suitable distances. 
The range in source-receiver separation was taken to be 55* to 86* to ensure separation of S 
and ScS at the longer distances and to avoid triplication in SS at shorter distances. The SS 
and S phases bottom from about 670 km to 2300 km depth. We performed a search over 
all earthquakes in the Harvard centroid moment tensor (CMT) catalog for the years 1977- 
1987 [Dziewonski et al., 1981; Dziewonski and Woodhouse, 1983] and over all GDSN, 
NARS, and GEOSCOPE digital seismic stations in order to find event-station pairs of the 
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proper epicentral distance which provide SS bounce points in the North Atlantic region. 
Epicenters were obtained from the Bulletin of the International Seismological Centre (ISC) 
for events occurring before 1987 and from the "Preliminary Determination of Epicenters 
of the U.S. National Earthquake Information Service (NEIS) for events occurring in 1987. 
The final distribution of sources and stations used to measure SS-S differential travel times 
is shown in Figure 2.3. The majority of data in this study comes from records of 
equatorial fracture zone earthquakes at North American and European stations, north and 
central Atlantic events at North American stations. Central American events at European 
stations, and Mediterranean and European earthquakes at North American stations. 

This search yielded over 2000 event-station pairs with the proper epicentral 
separation. After winnowing the list because of station inoperation, poor signal to noise 
ratio for the p hase s of interest, and interfering events, the final data set consists of nearly 
500 SS-S differential travel time residuals with bounce points in the north Atlantic (Figure 
2.4). Uncertainties are determined for each measurement following the procedure outlined 
in Appendix 2.A. A tabulation of all residuals is given in Appendix I. 


RESULTS 


We interpret the variations in SS-S differential travel times in terms of lateral 
velocity variations within the crust and upper mantle beneath the surface reflection points of 
the SS wave path. Kuo etal. [1987] and Woodward andMasters [1991] tested the validity 
of this assumption by plotting absolute S and SS residuals against SS-S residuals. They 
found that S and SS-S residuals are uncorrelated while SS and SS-S residuals are strongly 
correlated, indicating that the assumption is justified. The validity of this assumption is 
further supported by the strong correlation of SS-S times with surface tectonic features in 
the vicinity of the SS bounce point The residuals are further interpreted in terms of such 
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upper mantle processes as lithospheric aging, flow-induced anisotropy, and along-axis 
heterogeneity in mantle structure. 

Lithospheric Aging 

Cooling and thickening of the lithosphere should yield a tendency toward an 
increase in seismic velocity with increasing lithospheric age. A linear regression 
experiment was performed to examine the correlation of the SS-S residuals with seafloor 
age. A gridded map of seafloor ages was constructed for the north Atlantic from the 
magnetic anomalies of Klitgord and Schouten [ 1986] and ages assigned according to Kent 
and Gradstein [1986] and Klitgord and Schouten [1986]. The isochrons of Sclater et at. 
[1981] were used in a few regions which were not covered by the Klitgord and Schouten 
[1986] data set To obtain a representative age value for the region spanning approximately 
one horizontal wavelength of the incident (SS) wave, an average seafloor age was 
estimated for a 1* x 1* box centered on each SS bounce point. To reduce scatter, 
measurements whose bounce point depths differed by more than 2500 m from the depth 
predicted by the Parsons and Sclater [1977] plate cooling model were excluded from the 
final age regression. Although each SS wave samples the upper mantle at a finite range of 
lithosphere ages, we expect that the different travel time anomalies contributed by the SS 
path segments on the younger and older sides of the bounce point approximately cancel so 
that the age at the SS bounce point is appropriate to the associated SS-S residual. 

The SS-S residuals for the north Atlantic are consistent with the expectation of an 
increase in seismic velocity with seafloor age. For bounce points between 0* and 60*N 
latitude, the coefficient derived by linear regression of residual with square root of age is 
-0.68 ± 0.08 s My* 1/7 from 0 to 100 My, with a linear correlation coefficient of -0.85 
(Figure 2.5). However, residuals from 60 - 90*N do not seem to be strongly correlated 
with lithospheric age. This may be due to the fact that this area is more tectonically 
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complicated than ‘normal’ oceanic lithosphere [e.g.. White , 1988; Zehnder and Mutter, 
1990], includes several ridge jumps, is in close proximity to continental regions, and does 
not closely follow the age-depth relation of Parsons and Sclater [1977]. Compared with 
the residuals for 0-60*N, those from 60-90'N are anomalously negative at young ages and 
anomalously positive at older ages. The slope of SS-S residual vs. square root of age for 
data from 0 to 60 *N is smaller than that inferred from S delays of intraplate earthquakes in 
the Atlantic by Duschenes and Solomon [1977] (two-way S delay = -1.2 s My and that 
reported by Kuo et al. [1987] (-1 s My l/2 ). It is larger, however, than the global average 
obtained by Woodward and Masters [1991] (-0.51 s My ^). We find that the residual-age 
relation is not constant over the entire north Atlantic, so that some of these variations in 

slope may reflect real geographic differences. 

We may compare the variation of SS-S residual versus age with that due only to 
lithospheric cooling. For a lithospheric structure given by the plate cooling model of 
Parsons and Sclater [1977], we may convert temperature variations to differences in shear 
velocity v s by adopting a value for cH' s /9T, which we take to be uniform and equal to -0.6 
m/s K* 1 [McNutt and Judge, 1990]. For a horizontal slowness typical of the teleseismic S 
and SS waves of this study (0.1375 s/km), the slope of the line best fitting the SS-S travel 
time delay versus age given by the plate cooling model over 0-100 My is then -0.64 ± 0.07 
s My' 1 ^ 2 , a result indistinguishable from the observed slope. This agreement indicates that 
the dependence of travel time residual on plate age can be explained entirely by lithospheric 
cooling. 

The trend of the travel time residual versus lithospheric age relation changes at 
about 100 My. After 100 My, the residuals appear to flatten out (Figure 2.5), in the same 
sense as the plate cooling model of Parsons and Sclater [1977]. Such a pattern may reflect 
the unmodeled effect of increased sediment or crustal thickness, or, as suggested by 
Parsons and Sclater [1977] may be partially the result of secondary convection which 



22 


supplies heat to the base of the plate at older ages. To avoid possible biases associated with 
any of these effects we shall restrict our analysis to data with bounce points on lithosphere 
less than 100 My. To look for other systematic variations in the residuals, we correct for 
age by removing the linear relation shown by the solid line in Figure 2.5. This correction 
is effectively a normalization of residuals to 22-My-old lithosphere (the zero crossing of the 
regression line). 

Anisotropy 

Another systematic velocity variation that has been suggested as a possible 
contributor to residual SS-S travel times is azimuthal anisotropy. Kuo et al. [1987] 
examined this phenomenon in detail and concluded that alignment of olivine cystals in the 
asthenosphere created a significant pattern of azimuthal anisotropy in SS-S residuals 
measured in the Atlantic region. We have also searched for evidence of azimuthal 
anisotropy with our data set. 

Backus [1965] and Cr ampin [1977] demonstrated, from the general form of body 
wave anisotropy in a weakly anisotropic medium, that the linear form of the azimuthal 
variation of velocity is given by 

V 2 = Ao + Aj cos 20 + A 2 sin 20 + A 3 cos 40 + A 4 sin 40 (2.1) 

where V is the body wave velocity, the A n are linear functions of the elastic moduli, and 0 
is an azimuth, defined for our problem by the angle between the great circle path and the 
direction to geographic north measured at the SS bounce point. Equation (2.1) was further 
simplified by Kuo et al. [1987] and parameterized in terms of travel time residuals: 


R = Ro + R! cos 20 + R 2 sin 20 + R 3 cos 40 + R 4 sin 40 


( 2 . 2 ) 
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where R is the travel time residual and the R„ are constants. By fitting a function of this 
form to our age-corrected measurements we can determine if our data are consistent with 
the presence of anisotropy. 

We have conducted several tests of azimuthal anisotropy with our travel time data. 
We performed least squares inversions to determine 20 and 40 patterns which provide best 
fits to the age-corrected SS-S residuals. The anisotropy indicated by our regression 
experiments differs significantly from the preferred model of Kuo et al. [1987] both in 
magnitude and in phase (Figure 2.6). Our results indicate that for the 20 model the slow 
direction for SS-S is N4*W and the peak-to-peak magnitude of the effect is less than 1 s; 
for the 40 model the slow directions are N32*W and N58*E and the magnitude is 2.5 s; for 
the joint 20 and 40 model the slow direction is N32*W and the magnitude is just under 3 s. 
Kuo et al. [1987] obtained a peak-to-peak variation with azimuth of 5-7 s and a slow 
direction at N13°W. The slowest residuals in the Kuo et al. [1987] study were from north- 
south paths, i.e., nearly along the ridge, and the fastest residuals were from northeast- 
southwest-trending paths with bounce points north of the Azores-Gibraltar plate boundary 
(an area noted to be anomalously fast in their study), so their reported anisotropy may have 
been at least partly the result of unmodelled upper mantle heterogeneity. Our inversion for 
a 20 pattern of anisotropy provided a variance reduction of only 2%, compared with 20% 
for a 40 pattern, and 22% for a combined 20 and 40 pattern. On the basis of these values 
of variance reduction and the number of free parameters involved, our results suggest that 
there is no single coherent pattern of upper mantle anisotropy in the north Atlantic. The 
latest anisotropic upper mantle models obtained from surface wave tomography [Montagner 
and Tanimoto, 1990] also show a complex pattern of anisotropy in the region. Any 
azimuthal anisotropy in the asthenosphere induced by plate motions in the north Atlantic 
may be heterogeneous because the three plates in the region are slow-moving and the return 
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flow is not closely related to plate divergence [Hager and O'Connell, 1979, 1981; 
Parmentier and Oliver, 1979]. 

Spatial Patterns of Age-corrected Residuals 

After removal of the dependence chi seafloor age, a plot of SS-S travel time 
residuals at the SS surface reflection point (Figure 2.7) shows several interesting features. 
Perhaps the most striking is that residuals in the western Atlantic north of about 35* N are 
on average nearly 4 s more negative than those to the south. This feature is also noticeable 
in Figure 2.4 but is more obvious after age-dependence is removed. A similar change at 
approximately this latitude was noted for SS-S residuals with bounce positions on the 
eastern side of the Mid-Atlantic Ridge by Kuo etal. [1987] and was attributed to a change 
in upper mantle structure across the Azores-Gibraltar plate boundary. The signal we 
observe is predominantly from data with bounce points on the western side of the ridge. A 
map view of the azimuthal distribution is shown in Figure 2.8 and serves as an aid to 
assess qualitively the geometry of wave paths to the south and north of 35*N. We 
examined the possibility that this signal may be from the Caribbean anomaly, a region of 
anomalously high velocity in the mantle between 600 and 1400 km depth beneath the 
Caribbean originally reported by Jordan and Lynn [1974] and further confirmed by Grand 
[1987]. If the first leg of the SS rays propagating to western Europe were to bottom in the 
high velocity Caribbean region, the result would be early SS-S residuals. This would 
produce a feature of opposite sign from that observed, so we discount it as an influence 
here. Another possible explanation for the long-wavelength signal could be azimuthal 
anisotropy, but the examination above of possible patterns of azimuthal anisotropy does not 
support this suggestion. 

Another distinctive feature of the residuals in Figure 2.7 is a row of negative 
values which trends northwest to southeast along the trend of the New England Seamounts 
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and across the ridge to the vicinity of the Great Meteor Seamount This feature comes from 
event-station pairs at a number of different azimuths and distances so cannot be attributed to 
a source or receiver effect. We do not observe distinctive anomalies in the vicinity of the 
Bermuda, Azores or Canary Islands hotspots. The data density is poor for the Bermuda 
region, however, and any signal associated with the Canary Islands may be obscured by 
the ocean-continent transition. Recently active hotspot islands might be expected to display 
strong positive (late) residuals, such as Stewart and Keen [1978] observed for PP-P 
residuals at the Fogo Seamounts. In contrast. Woodward and Masters [1991] found 
mostly negative (early) SS-S residuals in the vicinity of the Hawaiin hotspot, and Jordan 
[1979] and Sipkin and Jordan [1980b] have suggested that the net effect of hotspots may 
be to produce early arrivals because of the presence of high velocities in a depleted mantle 
residuum. 

There is a systematic variation of SS-S residual with latitude, i.e., effectively along 
the direction of the Mid-Atlantic Ridge axis. Age-corrected SS-S residuals with SS bounce 
points on lithosphere younger than 100 My are shown versus latitude in Figure 2.9. The 
along-axis variations show a variety of scales, notably at wavelengths of about 1000 - 2000 
km in the region from 15* to 35*N, and at about 6000 km wavelength from late (positive 
residuals) in the south (20-35*) to early (negative residuals) farther north (45-55*N). The 
largest of these variations are robust with respect to selective removal of portions of the 
Hata The Iceland region appears as a local maximum (positive SS-S delay) on the profile, 
but the Azores hot spot does not have a distinct seismic signal. 

JOINT INVERSION OF TRAVEL TIME RESIDUALS AND GEOID AND DEPTH ANOMALIES 

Long- wavelength variations in shear wave velocity of the sort depicted in Figure 
2.9 pres umab ly are a consequence of some combination of variations in temperature and 
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composition of the upper mantle. Such lateral variations should also have signatures in 
other physical quantities measurable at these wavelengths, notably gravity (or geoid height) 
and topography (or residual bathymetry), because of the dependence of these quantities on 
bulk density. Travel time residuals, geoid anomalies, and residual depth anomalies are 
independent quantities dependent in different ways on temperature, bulk composition, and 
their variation with depth. We therefore seek a quantitative procedure for treating travel 
time residuals jointly with geoid and bathymetry data and in particular for a combined 
inversion of all three quantities for horizontal variations in upper mantle temperature and 
composition. 

To ensure complementarity of data sets, bathymetry and geoid height values are 
obtained at each SS bounce point, and both are corrected for subsidence with seafloor age 
by means of the plate cooling model [Parsons and Sclater, 1977; Parsons and Richter, 
1980]. In this manner we effectively normalize all observations to zero age. Bathymetric 
data are obtained from the corrected Digital Bathymetric Data Base (5’ grid) [C/.5. Naval 
Oceanographic Office, 1985]. Geoid data are taken from a combined set of Seasat and 
GEOS3 altimeter data [Marsh et al., 1986]. Data north of 70° N were not included in the 
Marsh et al. [1986] data set due to the high probability of being over sea ice, so our 
analysis below is confined to latitudes less than 70*N. We find that the correlation of SS-S 
residuals with the low order geoid is negative, but that at high order the correlation is 
positive (Figure 2.10). This relationship may indicate a depth dependence of contributions 
to geoid and travel time (e.g., the long wavelength signal may be a lower mantle effect). 
Low degree harmonics are likely linked to deep-seated density heterogeneities and 
subducting slabs [Hager, 1984; Hager et al., 1985]. Since we are interested in upper 
mantle processes, we filter out the long-wavelength component of the geoid by subtracting 
a reference field [Lerch et al., 1979] expanded in spherical harmonics to degree and order 7 
and tapered to degree and order 1 1. To provide a comparable bathymetric data set, 
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bathymetry is high-pass filtered (comer at 4000 km, cutoff at 6000 km) to remove long- 
wavelength trends. Along- axis profiles are constructed from the age-corrected and filtered 
geoid and bathymetry data. 

Profiles of age-corrected travel time residuals, geoid, and bathymetry are compared 
in Figure 2. 1 1 . While qualitative correlations among profiles are apparent, we seek to 
quantify possible models of temperature and compositional variations that can match these 
observations. Oceanic bathymetry and geoid height are both sensitive to variations in 
man tle density at depth. Such variations can be either thermal or compositional in origin 
and, like seismic velocity, are presumably related to mantle convection and differentiation. 
For a given density change, the seismic signature of thermal and compositional 
heterogeneity are of opposite sign, so travel time residuals constitute key information for 
distinguishing between mechanisms of heterogeneity. 

Inversion for Thermal Structure 

We seek to form ula te an inversion for the distribution of temperature anomalies 
T(x,z) (where x is along-axis and z is depth) that can produce the along-axis geoid, 
bathymetry, and travel-time anomalies shown in Figure 2.1 1. Topography and geoid 
kernels were calculated for prescribed models of viscosity for an incompressible 
Newtonian mantle with free slip at the surface and the core-mantle boundary. The 
convccting region is assumed to be overlain by a high-viscosity layer 40 km thick. We 
performed calculations both for a mantle of constant viscosity and for a mantle with a 
shallow low-viscosity layer. Topography and geoid anomalies depend on the viscosity 
structure, but the predicted travel times do not Kernels were calculated using a method 
similar to that of Richards and Hager [ 1984] except that the solution was directly integrated 
across the layers instead of being obtained via propagator matrices [McNutt and Judge, 


1990]. 
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The inversion is best conducted in the horizontal wavenumber domain. The 
thermal anomalies AT(k,z) at depth are related to the predicted dynamic topography h(k) for 
wavenumber k via an integral of the form 


Ah (k) = 



H (k,z) AT (k,z) dz 


(2.3) 


[Parsons and Daly, 19831 where a is the volumetric coefficient of thermal expansion, po 
and p w are the densities of the mantle and of water at standard temperature and pressure, 
and Zmin and Zmax are the upper and lower boundaries of the layer in which temperatures 
are allowed to vary. Table 2.2 contains a summary of the constants adopted here. The 
depth and wavenumber-dependent topography kernel H(k,z) is calculated from the 
equations of continuity and motion given a set of boundary conditions, a viscosity model, 
and a constitutive relation between stress and strain [Parsons and Daly, 19831. Similarly, 
the kernel G(k,z) for the geoid relates the thermal anomalies to the geoid N(k) via 

z u*» 

27irp n a r 

AN(k) = I G(k,z) AT(k,z) dz (2.4) 

Z am 

[Parsons and Daly, 19831 where T is the gravitational constant, and g is the surface 
gravitational acceleration. 

Sample geoid and topography kernels calculated for different wavenumbers and 
viscosity structures are shown in Figures 2.12 and 2.13. Cartesian kernels are used 
throughout this study because of their computational efficiency and straightforward 
application to Fourier transform techniques. We have compared extrema of the upper 
mantle portions of the geoid and topography kernels for a layered cartesian Earth and a 
spherical Earth for a number of wavelengths and different viscosity structures (Figure 
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2.13), and wc note good agreement even at very long wavelengths (spherical harmonic 
order C= 6). This agreement suggests that the results presented here should be applicable 
to the spherical Earth without introducing unreasonably large errors. 

Temperature perturbations at depth can be converted to a seismic velocity 
perturbation by assuming a value for the partial derivative of shear wave velocity with 
respect to temperature, dvV9T. The resulting two-wave travel time perturbation is given by 


At (k) 


= 2^1 f- 

ar J „ 


Z . 
mm 


AT (k, z) d 2 
V s (z ) 2 (1 - p 2 v s (z) 2 ) 


1/2 


(2.5) 


where v s (z) is from the reference shear velocity model [Dziewonski and Anderson, 1981] 
and p is the ray parameter, generally taken to be a representative value for the range of 
epicentral distances considered here. We use a value of -0.6 m/s K*1 for dv$/dT. This 
value is higher than the values of Anderson et ai. [ 1968] and Kumazawa and Anderson 
[1969] at standard temperature and pressure but is similar to the value of -0.62 m/s K' 1 
determined by McNutt and Judge [1990] by a least squares fit of Love-wave phase 
velocities to predicted temperature of the lithosphere. Such a value is consistent with the 
change in P-wave velocity with temperature, dvp/0T = -0.5 m/s K* 1 , found from modeling 
wave propagation along subducting slabs [Creager and Jordan, 1986; Fischer et al. t 1988] 
if we assume that dvg/9T =1.1 dvj^3T [Woodhouse and Dziewonski, 1984]. Partial melt 
would increase the value of dvs/3T [Sleep, 1974; Sato and Sacks, 1989], but simultaneous 
analysis of both shear and compressional differential travel times by Woodward and 
Masters [1991] indicates that significant partial melting is not required to explain the 
differential travel time residuals in the north Atlantic region. 

The forward problem consists of calculating geoid, topography, and travel time 
residual profiles given a starting two-dimensional temperature structure T(x,z). The 
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inverse problem consists of finding a temperature structure that predicts (via equations 2.3 - 
2.5) geoid, topography, and travel time profiles which best fit those observed. The 
familiar matrix equation d = A m is formed from discrete versions of equations 2.3 • 2.5. 
The data vector d consists of the topography, geoid, and travel time residuals, the model 
vector m contains the temperature variations for which we are solving, and the matrix A 
contains the coefficients and kernels which relate the data to the model. As a check on our 
procedure, we constructed a forward problem for geoid and topography and found good 
agreement with the modelling results of McKenzie et al. [1980]. 

The bathymetry, geoid, and travel time profiles of Figure 2. 1 1 are interpolated to a 
constant spacing, demeaned, tapered at both ends with a 10% sine squared taper, and 
Fourier transformed. Since our profile extends from 10 to 72°N, the first and last 10% of 
the profile (10 - 16°N and 66 - 72°N) will be affected by the taper. The 3n x 1 data vector d 
is then constructed, using the complex (to retain both amplitude and phase) bathymetry, 
geoid, and travel time data sampled at n discrete wavenumbers: 

d = [ Ah(k 1 ),...,Ah(k n ),AN(ki),...,AN(k„),At(ki),..., At(k n ) ] T (2.6) 


where T denotes transpose, and n in this case is equal to 5, representing the first 5 
coefficients of the Fourier series expansion (wavelengths 7104 km, 3552 km, 2368 km, 
1776 km, and 1420 km). For the case where temperature perturbations are constrained to 
be in a single layer, the n x 1 model vector m is given by 

m = [ AT(ki),...,AT(k„)]T (2.7) 


For the more general case of a multi-layer system, the nj x 1 model vector m is given by 
m = [ AT(ki.zi) AT(ki,Zj), AT(k 2 zi),...,AT(k 2 Zj),AT(k„ zi)...,AT(k„,Zj)] T (2.8) 
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where zj is the layer index and j is the total number of layers. In this paper we perform 
inversions for single-layer models only. The “layers” of temperature variations are 
independent of the “layering” system erf lid, low viscosity zone, and mantle which we use 
for the calculation of kernels, although major changes in viscosity would tend to segment 
AT as well. The temperature layering simply refers to that region bounded by Zmin and 

Zhuk in the integrals of equations 2.1 - 2.3. 

The 3n x n matrix A contains the coefficients and kernels that relate the temperature 


perturbations to the observations, which for the single-layer case is given by 
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The matrix A contains both bathymetry and topography kernels and is thus viscosity 
dependent; i.e., a viscosity structure must be assumed. We solve the equation d = A m by 
least squares 

m = ( A A ) A d (2.10) 

where A is the complex conjugate transpose of the matrix A. Construction of the data 
covariance matrix R<Jd is discussed in Appendix 2.B. Equation (2.10) is solved for the 
solution vector m, and variance reduction is calculated via 


variance reduction = 1 - 


(d-Am) R^d (d-Am) 

dR^d 


( 2 . 11 ) 


The resulting model vector m is inverse Fourier transformed back to the spatial domain to 
produce an along-axis temperature profile. The solution m is also substituted into 
equations 2.3 - 2.5 to compare predicted geoid, bathymetry, and travel time residuals with 
those observed. 

Six inversion experiments few temperature structure were performed (Tables 2.3 
and 2.4). Table 2.4 contains spectral coefficients for both the observed data and the 
predicted models, and serves as a guide to how well the various spectral components of the 
daxa are being fit Since the geoid and bathymetry kernels do not include any phase 
information (except that a sign change produces a 180* phase shift), large phase differences 
between components of the observed geoid and bathymetry would indicate that these 
components cannot simultaneously be well fit by our models. Inversions were carried out 
for two different viscosity models and for three different thicknesses of the layer in which 
lateral temperature variations were assumed to occur. Because topography and geoid 
ano malie s depend only on the ratios of viscosity in different layers [Richards and Hager, 
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1984; Robinson et al., 1987; Hong et al ., 1990], we set the dimensionless viscosity of the 
layer representing the bulk of the mantle to unity. In one viscosity model, termed the 
“constant viscosity mantle,” a 40-km-thick high-viscosity lid overlies a unit viscosity 
mantle. We set the viscosity in the lid to 10 4 , which effectively mimics rigid behavior. In 
a second model, a 160-km-thick low-viscosity zone is present beneath a 40-km-thick lid; 
the viscosity in the low-viscosity zone is a factor of 100 less than in the underlying mantle. 
The thickness of the layer of temperature perturbations was taken variously to extend from 
0-150 km depth, 0-300 km depth, and 0-650 km depth. The matrix A is different for each 
of these cases, as it involves viscosity-dependent geoid and topography kernels and also a 
summation over depth. 

Inversion results for the constant- viscosity-mantle cases are shown in Figure 2.14. 
The “observed” profile is actually a filtered version of the observations, containing only the 
wavelengths used in the inversion (1400 to 7100 km). Predicted profiles were calculated 
from equation 2.5. For these solutions, the long-wavelength fit to geoid is better than at 
short wavelengths. The fit to bathymetry is poor. The predicted magnitude of the SS-S 
residuals range from a factor of 5 too small for the 650-km-thiek layer to a factor of about 
1.5 too small for the 150-km-thick layer. Increasing the temperature variations to improve 
the fit to the SS-S residuals leads to predicted geoid variations that are too large. The 
highest total variance reduction and best fit for the constant-viscosity cases come when 
lateral variations are constrained to shallow (0-150 km) depth. The variance reduction is 
25% for bathymetry, 79% for geoid, and 58% for travel times (Table 2.3). The total 
variance reduction is 53%. The variation in temperature is 180 K for the 150-km-thick 
layer, and only 60 K for the 300-km-thick layer. 

Figure 2.15 shows inversion results for the models with a thin low-viscosity zone. 
A good fit to both geoid and travel time is found, although the alignment in phase of 
predicted and observed geoid is not as good as for the constant-viscosity case. The fit to 
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bathymetry is again poor. The total variance reduction for the 150- km- thick and 300- km- 
thick layers are both 57%, although the shallow model provides slightly higher variance 
reduction for bathymetry (27% for 0*150 km deep layer, 24% for 0-300 km deep layer) 
and the 300- km- thick-layer model provides higher variance reduction for geoid (79% for 0- 
150 km deep layer, 85% for 0-300 km deep layer). The variation in temperature for the 
150-km-thick layer is 230 K and in the 300- km- thick layer is 1 10 K. 

We have explored the hypothesis that the lack of correlation of predicted and 
observed topography is an indication that the source of variations in the geoid and travel 
time anomalies is deep. To test this hypothesis, we performed inversions with temperature 
variations restricted to deeper layers and found that fits to topography were still poor. It is 
possible that the bathymetric signal is dominated by crustal thickness variations which are 
not included in our calculation of dynamic topography. An assessment of such thickness 
variations is discussed further in Appendix 2.B. 

Inversion for Compositional Variations 

A possible alternative to along-axis variation in mantle temperature is lateral 
variation in bulk mantle composition, due perhaps to a variable extent of melt extraction or 
different degrees of mixing of compositionally distinct volumes of mantle material. The 
dynamical effects erf compositionally induced density variations can be large [ O’Hara , 

1975; Boyd and McCallister, 1976; Oxburgh and Parmentier, 1977; Sotin and Parmentier, 
1989]. The fraction of mantle potentially extractable as basaltic melt is thought to be 15- 
25% [e.g.. Green and Liebermann, 1976]. Thus, for every volume of basalt removed 
from the mantle, a volume of residuum several times larger is left behind. The effect of 
basalt depletion is to increase the molar ratio Mg/(Mg + Fe) (or Mg#) in the residuum, 
which reduces the density and increases the seismic velocities (e.g., Liebermann, 1970; 
Akimoto, 1972]. For example, subtraction of 20 mole % olivine basalt from pyrolite can 
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decrease the density of the residuum by nearly 2%, equivalent to a thermal perturbation of 
nearly 500 K [Jordan, 1979]. Thus compositional changes need only be slight to produce 
effects on the order of 100 K, comparable to values obtained from the inversions for 
temperature variations. In this section we explore the effects erf compositional variations 
parameterized in terms of the variation in the Mg# in the upper mantle along the ridge. Our 
motivation for parameterizing compositional variations simply in terms of Mg# is that 
differences in this quantity yield significant variations in seismic velocity and density, in 
contrast to most other measures of degree of melt extraction. 

Partial derivatives of density and seismic velocity with respect to Mg# are obtained 
from Akimoto [ 1972]. These values were measured on a suite of samples ranging from 
pure forsterite (Mg 2 Si0 4 ) to pure fayalite (Fe 2 Si 04 ). While these partial derivatives are at 
standard temperature and pressure, it is expected that a change to elevated temperature and 
pressure will have only a second order effect, since temperature and pressure corrections 
work in opposite directions [Jordan, 1979]. Above the solidus temperature, however, the 
amount and distribution of partial melt, which may depend strongly on composition and 
particularly volatile content, is important The presence of melt is likely to have a larger 
effect on shear wave velocities than on bulk density. Calculations of melt migration, 
however, suggest that once created, melt segregates rapidly by a percolation mechanism 
[e.g. Scott and Stevenson , 1989], so that the melt fraction present in the mantle at any 
given time is probably small. Studies of mantle peridotites [Johnson et al., 1990] also 
support the importance of fractional melting. 

It is straightforward to convert equations (2.3) and (2.4) to relations between geoid 
or topography and a compositionally induced density perturbation by means of the relation 

Ap = - p 0 a AT 

Compositional anomalies at depth yield a dynamic topography h(k) given by 


( 2 . 12 ) 
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_ Zmax 

Ah(k) = J H(k,z) AMg(U) dz (2.13) 

z . 

min 

where AMg represents the fractional change in the Mg#. Compositional anomalies yield a 
geoid anomaly 


AN(k) = 


Ml 

gk dMg 


Z n*X 

J G(k,z) AMg(k,z) 
z ™ 


dz 


(2.14) 


For a compositional perturbation at depth the resulting two-wave travel time perturbation is 
given by 


At(k) = 2 


Z ran 


3v s f AMg(k,z) dz 

9M V v s (z) 2 (l-p 2 v s (z ) 2 ) 1/2 

nan 


(2.15) 


Using equations (2.13) - (2.15), an inversion scheme similar to that used for 
thermal perturbations is formed. The solution vector now has the form 

m = [AMg(ki), ..., AMg(k n )] T (2.16) 

The data vector remains the same as in equation (2.6), while the matrix of coefficients, A, 
changes to reflect the relation between the data and mantle composition, rather than 
temperature, as outlined in equations (2.13) - (2.15). 

The results of the inversions for compositional variations are summarized in Table 
2.3 and in Figures 2.16 and 2.17. We are unable to match simultaneously both SS-S travel 
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time residuals and geoid and bathymetric anomalies with solely mantle compositional 
variations for either a constant viscosity mantle or one with a low viscosity zone. This is 
not surprising, as the travel times are for the most part positively correlated with geoid and 
bathymetry, but compositional variations (at least for the Mg 2 Si 04 - Fe 2 Si 04 system 
examined here) have an opposite effect on travel time and geoid-ba thyme try. 

For the constant viscosity case, the fit to the geoid is excellent, and the fit to 
bathymetry is slightly better than in the inversion for temperature. The fit to SS-S residuals 
is so poor that the variance reduction is negative for travel time. Large compositional 
changes would be required to affect travel times, whereas only small compositional 
changes are needed to produce significant density contrasts to match the geoid signal. The 
total variance reduction for the constant viscosity case does not vary greatly (from 32-33%) 
for compositional changes constrained to be over different depth intervals, though the 
variance reductions for individual data sets (bathymetry, geoid, travel time) vary 
significantly from model to model (see Table 2.3). The range in Mg# is about 1% if the 
variation is constrained to the depth range 0-150 km and only 0.1% for the 0-650 km depth 
range. 

Figure 2.17 shows inversion results for the model with a low viscosity zone. A 
good fit to both geoid and bathymetry is found, although the alignment of predicted and 
observed geoid is not as good as in the constant viscosity case. The fit to bathymetry is the 
best of any models so far. The total variance reduction is still low (43 to 49%), due to the 
poor fit to travel times (negative variance reduction in all cases except the 0-650 km model). 
The range in Mg# is 2.4% if constrained to 0- 150 km depth, 1.3% over 0-300 km depth, 
and 0.5% over 0-650 km depth. 

Joint Inversion for Temperature and Composition 

We next explore whether a combination of temperature and compositional 
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variations can provide a good match to the observed geoid, travel time, and bathymetry. 
Joint inversions provide improved fits to all data at the expense of introducing additional 
free parameters. For these inversions the data vector remains the same as in equation (2.6), 
the solution vector is modified to include both temperature and composition, and the matrix 
of coefficients. A, includes the effects of both temperature and composition. The matrix- 
building equations become, for example, for topography, 

^ax z m*x 

AhOO- J H(k,z) AT(k,z) dz + f H(k,z) AMgrk.z) dz (2 . 17) 

z . z . 

mm mm 

which is simply a combination of equations (2.3) and (2.13). The new geoid equation 
comes from a combination of equations (2.4) and (2.14) and the travel time equation from a 
combination of equations (2.5) and (2.15). Cross terms, such as compositional changes 
induced by increases or decreases in temperature, are neglected. 

The results for the joint inversion for temperature and composition are summarized 
in Table 2.3 and Figures 2. 17 and 2.18. The travel time residuals are well-modeled in all 
cases, as are the geoid data. The topography is best fit for the case with a low viscosity 
zone. Resolution of the depth interval of the most important lateral variations is rather 
poor. The topography is fit marginally better for the case where temperature and 
compositional anomalies are constrained to be shallower than 300 km. For the constant 
viscosity mantle, the temperature variations range from 210 K, if constrained to 0-150 km 
depth, to 55 K if over 0-650 km depth; variations in Mg# range from 1.5% if over 0-150 
km depth to 0.4% for 0-650 km depth. For the case with an upper mantle low viscosity 
zone, the temperature variations are similar to those in the constant viscosity case, but the 
variations in Mg# are larger, from over 2% for 0-150 km depth to nearly 1% for 0-650 km 
depth. 

The travel time residuals are perfectly fit in the joint inversions for temperature and 



39 


composition (Tabic 2.3). This occurs because of the way the model parameters act in a 
similar manner on both gcoid and bathymetry, producing a singular matrix if only geo id 
and bathymetry data are inverted for both temperature and composition. Undamped least 
squares always provides perfect solutions when the number of equations is equal to the 
number of unknowns unless the matrix to be inverted is singular. If we perform an 
inversion including only travel time and geoid data, we have the same number of equations 
as unknowns, the matrix is nonsingular, and we obtain perfect fits to both travel time and 
geoid. Similarly, if we perform an inversion of travel time and bathymetry data, we again 
obtain perfect fits to both data sets. If we perform an inversion of geoid and bathymetry 
data, however, we are unable to obtain solutions without applying damping. In the joint 
inversion of travel time, geoid, and bathymetry, we have more equations than unknowns 
and the inversion is overdetermined. However, the travel times are perfectly determined in 
this case because of the nonuniqueness inherent with geoid and bathymetry. We have 
performed undamped inversions with various weightings on the geoid, bathymetry, and 
travel time data, and in all cases the travel times remain perfectly fit 

We have also performed joint inversion for temperature and composition with Mg# 
variations constrained to be in the upper 50 km of the lithosphere so as to mimic 
compositional variations due solely to variable melt extraction at the ridge. Temperature 
perturbations were allowed to remain within the depth ranges adopted earlier. The results 
for this inversion are summarized in Table 2.3 and Figures 2.20 and 2.21. The variance 
reduction was similar for the constant viscosity case and for the model with a low viscosity 
zone. In general, the geoid is fit very well, the predicted amplitudes are a bit low for travel 
Hme residuals, and the topography fit is slightly out of phase. For the constant viscosity 
mantle, the range in temperature is 210 K over 0-150 km depth and 25 K over 0-650 km, 
while Mg# variations constrained to be confined to 0-50 km depth were over 5%. For the 
case with a low viscosity zone, the temperature variations were not dramatically different 
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from those in the constant viscosity case, and variations in Mg# were about 4.5%. The 
inversion solution shows high temperatures near 30’N and low temperatures in the region 
from 50-60*N. Iceland also appears to be underlain by high-temperature mantle. Going 
from south to north along the ridge, compositional variations indicate low Mg# in the 
vicinity of 20-30*N, high Mg# in the Azores region (40*N), low values near 50*N, and 
high values near 60* N. 


DISCUSSION 


The temperature and compositional variations in Figures 2.14-2.21 are broadly 
consistent with observed travel time, geoid, and bathymetry anomalies in the north Atlantic 
region. Temperature variations alone can account for most of the observed anomalies. In 
contrast, compositional variations alone cannot match all anomalies simultaneously. We 
infer that a component of the observed anomalies is due to long-wavelength variations in 
upper mantle temperature. Joint inversions for temperature and composition provide better 
fits than single-variable models, but at the expense of introducing additional free 
parameters. 

It is difficult to select a ‘best’ model from the suite of inversions presented. The 
variance reductions in Table 2.3 serve as a guide, but independent criteria may allow us to 
reject some of the models, even those with high variance reductions. In particular, those 
models with large temperature variations (well in excess of 100 K) can be seriously 
questioned. Lateral temperature variations at upper mantle levels beneath oceanic ridges are 
thought to be no more than about 300 K globally [Klein and Langmuir, 1987], so a 
variation in temperature of 230 K (as in the inversion with a low-viscosity zone and a 150- 
lan-thick layer of temperature perturbations) solely within a section of the north Atlantic is 
probably unreasonably large. Further, as White and McKenzie [1989] have noted. 
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relatively small increases in mantle temperature above values typical for the mid-ocean ridge 
are sufficient to cause large increases in melt production. Their models indicate that, for 
fixed bulk composition, an increase of 100 K above normal doubles the amount of melt 
while a 200 K increase can quadruple it Such increased melt production should lead to 
appro xima tely corresponding increases in crustal thickness. Variations in oceanic crustal 
thickness away from fracture zones, however, are generally thought to be small, with 
thicknesses typically 6-7 km and ranging from 4.5 to 8.5 km [Spudich and Orcutt, 1980; 
White , 1984; Purdy and Detrick, 1986]. In the joint inversion for temperature and 
composition, temperature variations if confined to 1 50 km depth are excessive (over 200 
K) and if the variations extend over 0-650 km the fit to topography is poor, especially for 
the constant-viscosity mantle. On the basis of these results we prefer the models with 
temperature variations occurring over 0-300 km depth. For the constant viscosity mantle, 
the temperature variation is 1 10 K, and the variation in Mg# is 0.75%. For the case with 
an upper mantle low viscosity zone, the predicted temperature variation is 125 K, and the 
variation in Mg# is 1.1 %. The total variance reduction is greater in the model with a low 
viscosity zone. 

Even a temperature variation of about 100 K is high for a mantle of constant 
composition, since we do not observe increased crustal thickness in regions that our 
models indicate have high temperatures. The assumption of approximately constant upper 
man tle composition warrants discussion. In particular, lateral variation in trace amounts of 
man tle volatiles may have a large effect on seismic velocity at a given temperature. The 
presence of even a slight amount of water, for instance, is sufficient to cause a significant 
decrease in the initial melting temperature of peridotite [Wyllie, 1971]. Estimates of volatile 
contents and their lateral variations in the north Atlantic region have been made from 
measurements of abundances of halogens, Si02, K 2 O, and H 2 O in basalts and from the 
volumes of vesicles in basalts [Schilling et al., 1980, 1983; Schilling, 1986; Michael, 
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1988]. These studies indicate that Cl, Br, F, and H 2 O contents increase toward the Azores 
and Iceland and that H 2 O is two to three times more abundant in Mid- Atlantic Ridge basalts 
erupted over the Azotes platform than at adjacent normal ridge segments. The effect of 
volatiles on density and shear wave velocity will be slight at subsolidus temperatures but 
can be major over the melting interval [ Goetze , 1977]. The presence of melt will act to 
decrease significantly the seismic velocity [Duschenes and Solomon, 1977] and, to a lesser 
extent, lower the density of the mande. To the extent that seismic velocity depends on 
proximity of the temperature to the solidus temperature [Sato etal., 1988, 1989], volatile 
content can trade off with temperature in its effect on velocity at subsolidus conditions. 
Thus, variation in volatile content could lessen the variations in melt production implied by 
the inversion solutions. 

Even without significant variations in volatile content, it is clearly an 
oversimplification to parameterize mantle composition in farms of only a single quantity. 
Further we have assumed that the partial derivatives of bulk density and seismic velocity 
with respect to Mg# that are those for olivine [Akimoto, 1972]. The work of Jordan 
[1979] indicates that these derivatives remain nearly constant for many different mande 
compositions (i.e., pyrolite-type compositions with various amounts of olivine, 
orthopyroxene, clinopyroxene, spinel, and garnet), so the latter assumption is sound. 
However, at any given Mg#, orthopyroxene and clinopyroxene have lower velocities and 
are less dense than olivine, while garnet and spinel are seismically faster and denser than 
olivine [Jordan, 1979], so an increase in the weight percent of orthopyroxene and 
clinopyroxene or a decrease in the weight percent of garnet and spinel with respect to 
olivine in the mande could counteract some of the temperature variations obtained undo' the 
assum ption of effectively uniform mineralogy. Several studies [Wood, 1979; Jaques and 
Green, 1980; Dick et al., 1984] have suggested that compositional variations in the mande 
are plausible. Indeed a number of workers [e.g., Davies, 1984; Alligre et al., 1984] favor 
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dynamic models for the mantle in which dispersed heterogeneities of various sizes and 
shapes are passively embedded in a continually mixed, converting mantle. Variations in 
modal fractions of olivine, orthopyroxene, and clinopyroxene in peridotites recovered 
along the Mid-Atlantic Ridge have been reported in several studies [Dick et al., 1984; 
Michael and Bonatti, 1985]. These variations are typically attributed to different degrees of 
mrlt extraction but could also be partially due to intrinsic upper mantle heterogeneity. For 
example, the relative fractions of olivine, clinopyroxene, and orthopyroxene indicated by 
Michael and Bonatti [1985] at 26* N and 30* N, if extended to depth, could counteract a 
portion of the temperature differences indicated by the inversion solutions for these 
regions. 

Chemical analysis of dredged peridotites in the north Atlantic indicate a range of 
about 2.5% variation in Mg# [Michael and Bonatti , 1985]. This value is intermediate 
between what we find for models with compositional variations constrained to be shallow 
(4.5 to 6% variation) and those models with compositional variations in the same depth 
ranges as the thermal variations (1-2%). This suggests that compositional variations may 
be concentrated slightly shallower than the temperature variations. Michael and Bonatti 
[1985] present an along-axis profile of Mg# variations from dredged peridotites which can 
be compared with our calculated profile. The main feature in their profile is a zone of high 
values of Mg# in the Azores region, from 34-45*N, relative to the rest of the ridge, 
consistent with our modelling results. Their data sampling is too sparse to delineate other 
long-wavelength features. Their average value for 26*N also has a high Mg# relative to 
adjacent data. This is consistent with our observation of early SS-S travel times and low 
geoid in this region. This anomaly is of too short a wavelength (< 1000 km), however, to 
resolve in our inversions. We should note that comparisons merit caution, as small scale 
features, such as those due to ridge segmentation, can produce large differences in 
composition between peridotites over scales of tens of kilometers. In addition, dredged 
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peridodtes are mostly from fracture zone environments, which may not be representative of 
typical ridge mantle [Dick, 1989]. 

On the SS-S residual profile the Iceland region appears as a local maximum Gate 
SS) but the Azores hotspot does not show a distinct seismic signal. The inversion results 
for these two regions are also markedly different. The results erf the joint inversion for 
temperature and composition predict a high Mg# in the Azores region while indicated 
temperatures are not anomalously high. At Iceland, in contrast, high temperatures 
dominate. Wait by Schilling [1986] and Bonatti [1990] outlines the differences in 
geochemical signatures between the Azores and Iceland hot spots. These workers suggest 
that Iceland is a "traditional" plume ha spot, with a predominantly thermal origin, but that 
the Azores might be more aptly named a “wet spot” because of the presence of excess 
hydrous phases and the lack of a thermal anomaly. Bonatti [1990] suggests that because 
the Azores hotspot is rich in volatiles, enhanced melt production could occur with little or 
no increase in temperature. The high Mg# indicated in our inversions allows the region to 
be seismically fast (as we observe) but of low density (as geoid and bathymetry require). 
The results are consistent with the hypothesis that the Azores ha spot is na associated with 
a pl um e- like thermal anomaly. Inversion of surface wave dispersion data can potentially 
provide further tests of these ideas, but studies to date have yielded apparently conflicting 
results. Results of several such investigations [Nakanishi and Anderson, 1984; Tanimoto, 
1990; Zhang and Tanimoto, 1990a] suggest that the Azores region is seismically slow at 
depths less than 300 km but a study utilizing 50-200-s-period Rayleigh waves by Mocquet 
et al. [1989] does not These differences may be partially attributable to the differences in 
wave periods employed and mode of analysis from study to study. It may be possible that 
what appear to be low velocities at the Azores are a result of horizontal smoothing of the 
low velocities along the ridge and have little to do with the actual structure in the Azores 
region. None of these long-period surface wave studies resolve a distinctive anomaly at 
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Iceland. Clearly, more work is needed to resolve the upper mantle velocity structure of hot 
spot regions. 

Bonatti [1990] has constructed profiles of the equilibrium temperature of dredged 
peridotites along the Mid-Atlantic Ridge axis from 0 to 60*N by means of two different 
geothermometers [Wells, 1977; Lindsley , 1983]. Comparison of these profiles with the 
along-axis temperature variations obtained from our inversions reveals a number of 
qualitative correlations as well as a few discrepancies. The range of temperature variations 
in the profile based on the Lindsley [1983] geothermometer is about 150 K, neglecting high 
values termed “anomalous.” When the high values are included the range increases to 400 
K. The profile utilizing the Wells [1977] geo thermometer has a range of 100 K neglecting 
the anomalous values and 350 K including them. The highest temperatures in our 
inversions are near 30* N (Figures 2.18-2.19), a region showing a slight peak in Bonatti’ s 
temperature profile estimated according to Lindsley [1983] and a very weak rise in the 
profile utilizing the Wells [1977] geothermometer. There is a small dip in temperature at 
26*N (a region which we find to be seismically fast) in the Lindsley [1983] and Wells 
[1977] profiles, but the difference may not be significant considering the error bars. 
Bergman and Solomon [1989] also found the upper mantle near 26’N to be seismically fast 
from an analysis of tcleseismic P-wave travel time residuals from earthquakes in this region 
recorded by a local ocean-bottom seismic network. The lowest temperatures on the profiles 
of Bonatti [1990] are at 43*N. Temperatures from our inversion solutions are also low in 
this region, although the Bonatti [1990] profiles indicate an increase in temperature 
proceeding north from 43*N to 53*N, whereas our results favor continued low 
temperatures. Part of the difference between our results and the geochemical studies may 
be attributed to the fact that the depth sampled by basalts and pendotites is likely to be 
shallower than the layer thicknesses of most of our models. Assuming that the 6-km-thick 
oceanic crust was formed by 9 to 22 % partial melting of the mantle [Klein and Langmuir, 
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1987], then the volume of residual peridodte will extend from the base of the crust to 
somewhere between 30 and 70 km depth. The amount of depletion will vary with depth if 
we assume a fractional melting model. Our models with compositional variations confined 
to depths less than 50 km are most representative of shallow fractionation and 
differentiation. 

Several improvements in future studies of the type presented here may be 
envisioned. Our models thus far have been limited to simply parameterized one- 
dimensional variations in temperature and composition within a single layer. It is likely that 
these lateral variations are not constant within a given layer and that there are two- 
dimensional lateral variations independent of lithospheric aging. The techniques outlined in 
this paper can be generalized to a multilayer system and to two-dimensional wavenumber 
(see equation 2.8), but we do not feel that the resolution of our data can justify more 
complicated models at this time. Kernels for seismic surface waves are strongly peaked in 
the upper mantle, and such data would provide a useful constraint in future models. The 
inclusion of surface wave data would help to distinguish between lithospheric and 
asthenospheric effects and may allow for two or more independently resolved layers. 
Extension of the modelling to three dimensions would permit an assessment of the degree 
to which mantle anomalies beneath the ridge extend off axis. Implicit in our age-correction 
is the assumption that the anomalous properties of the ridge mantle are steady state on a 
time scale of 100 My. Recent seafloor surveys and theoretical studies [e.g., Pockalny et 
al., 1988; Scott and Stevenson, 1989] bring this assumption into question and suggest that 
at least on short time scales (< 1 My) and at slow spreading rates (as in the Atlantic) 
intermittent periods of melting and crustal formation may be separated by periods with little 
or no melt production. These temporation variations are likely to be averaged out, 
however, over the typical horizontal wavelength (100 km) of a long-period SS wave. 

Another limitation of our models is that they depend on the assumed values of 
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several physical constants. It is straightforward, however, to estimate the effect of 
choosing different values. The viscosity structures we employ are also quite simple but 
have been chosen to represent two models widely invoked in other studies - a constant or 
nearly constant viscosity mantle [e.g., Peltier, 1989] and a mantle with a thin low viscosity 
layer [e.g., Craig and McKenzie , 1986; Robinson etal., 1987]. The viscosity structure of 
the Earth may be temperature and pressure dependent or vary laterally, but we have not 
considered viscosity structures of this type. Nor have we modelled the effects of partial 
melting which could accompany the temperature variations we predict The effect of 
retained melt on the physical properties of the mantle depends critically on the melt fraction 
and geometry, characteristics presently poorly known. Sato etal. [1988, 1989] downplay 
the importance of partial melt and suggest that most mande seismic velocity anomalies can 
be explained by temperature variations at subsolidus conditions. The combined analysis of 
both shear and compressional differential travel times also suggest that significant partial 
melting is not required to explain the travel time residuals in the north Adantic region 
[Woodward and Masters, 1991]. 


CONCLUSIONS 


We have measured 500 SS-S differential travel times for paths in the north Adantic 
region. The SS-S travel time residual decreases linearly with square root of age, in general 
agreement with the plate cooling model to an age of 80-100 My [Parsons and Sc later, 
1977]. Azimuthal anisotropy is not clearly resolved, and the azimuthal patterns of our data 
are not consistent with the preferred upper mande anisotropy model of Kuo et al. [1987] 
for the north Adantic. An along-axis profile of age-corrected travel time residuals displays 
significant long-wavelength variations, notably at wavelengths of 1000-2000 km. The 
largest of these variations are robust with respect to selective removal of portions of the 
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data. 

We have formulated a joint inversion of travel time residuals and geoid and 
bathymetric anomalies for lateral variation in upper mantle temperature and composition. 

On the basis of variance reduction, inversion few temperature favors the presence of an 
upper mantle low viscosity zone and temperature anomalies concentrated at depths less than 
300 km. We are unable to match travel time residuals simultaneously with geoid and 
bathymetry solely with lateral variations in bulk composition (Mg#). Joint inversions for 
temperature and composition provide good fits to both travel time and and geoid regardless 
of viscosity structure or layer depth and thickness, but the best fits to bathymetry come 
from models with a low- viscosity zone and thermal or compositional variations confined to 
shallow depth. The Mg# variations predicted in the joint inversion for temperature and 
composition are comparable to those found by Michael and Bonatti [1985] in a study of 
dredged peridotites along the Mid-Atlantic Ridge and may be related to variations in melt 
production along the ridge. 

The preferred inversion solutions have variations in upper mantle temperature 
along the Mid-Atlantic Ridge of about 100 K. For a constant bulk composition, such a 
temperature variation would produce about a 7 km variation in crustal thickness [White and 
McKenzie, 1989], larger than is generally observed [Spudich and Orcutt, 1980; White, 
1984; Purdy and Detrick, 1986]. Introducing compositional variations as well as 
temperature variations in the inversions does not change the range of temperature 
appreciably. The presence of volatiles in the mantle can have a strong effect on 
temperatures required for melting, and variations in volatile content along the ridge may 
reduce the large variation in melt production implied by the lateral temperature variations 
indicated in our models. 
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APPENDIX 2 A: ESTIMATION OF ERRORS FOR SS-S DIFFERENTIAL TRAVEL TIMES 

It is imp ortant to quantify the uncertainties in the differential travel time 
measurements. After cross-correlation, the “quality” of each individual SS-S measurement 
is rated and a grade is assigned. The cross correlation coefficient, which describes the 
degree of fit between the synthetic and real SS phases, is used as an objective aid in the 
assignment of quality. However, our final assignment of quality is largely subjective and 
based upon visual inspection of the "synthetic” SS, real SS, and cross correlogram, taking 
into account the sharpness of the arrivals and their alignment, the clarity of the seismogram, 
and the appearance of a single clear peak in the cross correlation function. An "A" quality 
grade indicates an excellent fit, "B" quality indicates good phase alignment but only a fair 
fit, and a "C" quality grade indicates a poor fit or some ambiguity as to phase alignment. In 
addition to A, B, and C grades, there were data that were rejected due to poor signal to 
noise ratio for either the S or SS phases. 

Assuming that the uncertainty in an individual measurement comes from a 
combination of measurement error, unmodeled lower mantle structure, and epicentral error, 
we write, for example, for the measurement variance of an "A" quality datum: 

<*A 2 = °Am 2 + <*lm 2 + °epi 2 (2A 1 } 

where o A is the total uncertainty, is the measurement error, ojn, is the uncertainty due 
to unmodeled lower mantle structure, and is the epicentral error. We assume that Ojm 
and Ocpi are the same for A, B, and C quality measurements, but the measurement error is 
obviously a strong function of data quality. 
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Effect of Epicentral Error 

In general, epicentral errors affect the travel times only slightly. The events used 
in this study were well recorded by a large number of stations over a wide range of 
azimuths, and typical epicentral mislocations are probably less than 10 km (which would 
yield a differential travel-time error of 0.35 s at 75* distance). The travel times are even 
less sensitive to errors in focal depth; an error in depth of 25 km contributes only about 0.3 
s to the SS-S residual. Using the rule of thumb that one standard deviation is about one half 
of the estimated extremes, we adopt = 0.75 s as a conservative estimate of epicentral 
error. 

Effect of Unmodelled Lower Mantle Heterogeneity 

We estimate the likely magnitude of lateral variations in the shear wave velocity of 
the lower mantle from models of lower mantle heterogeneity in P wave velocity (such as 
model L02.56 of Dziewonski [1984]). The average variation in travel times of direct P 
waves bottoming in the lower mantle is in the range ± 0.5 s. Global tomographic studies 
by Dziewonski and Woodhouse [1987] indicate that the scaling ratio (SV/V^ /(5Vp/V p ) * 

2 in the lower mantle. Such a scaling is also suggested by comparison of lower mantle P 
wave models with the recent lower mantle S model of Tanimoto [1990]. Assuming such 
an S to P velocity anomaly scaling, the resulting variation in S wave arrival time 
contributed by the lower mantle would likely be about ± 1.5 s, a fraction of the observed 
range in SS-S residual. While the major features of lower mantle model L02.56 
[Dziewonski, 1984] and the lower mantle portions of Tanimoto’ s [1990] model are for the 
most part similar, enough differences exist that the application of a lower mantle 
“correction” to our data might add more uncertainty than it removes. Further, absolute S- 
wave travel times do not show enough variance for us to suspect large lower mantle effects 
[e.g., Randall, 1971; Girardin and Poupinet, 1974; Hart and Butler, 1978; Uhrhammer , 



51 


1978, 1979], and the work of Gudmundson et al. [1990] indicates that most of the variance 
from the ISC tables is attributable to the shallow mantle, i.e., most of the Earth’s 
heterogeneity is in the upper mantle and the lower mantle is fairly homogeneous. On the 
basis of the above information, we set Ojn, = 0.5 s for our study. 

Measurement Error 

As an objective means to obtain error estimates, we examine the scatter in A, B, 
and C quality picks in a small region. We measured the root mean square (rms) difference 
between travel time residuals of the same grade (A, B, or C) with bounce points separated 
by less than 80 km and with differences in path azimuth at the bounce point of less than 
10*. An 80-km distance is less than the horizontal wavelength of SS (which is about 180 
km at 25 s period) so we do not expect much contribution to the rms difference from actual 
lateral variations in structure. The rms difference for the 16 A quality residual pairs which 
were within 80 km of each other was 1.15 s. For B quality picks, an rms difference of 
2.08 s was measured using 20 residual pairs, and for C quality picks 44 residual pairs 
yielded an rms difference of 2.96 s. 

We interpret these estimates of the rms diffences as representing the average 
overall errors in the A, B, and C grade measurements. (Unmodelled lower mantle structure 
should be nearly identical for data with bounce points within 80 km and at similar 
azimuths). Under this interpretation we can write, for A-quality residuals, 

<*Arms 2 = °Am 2 + <*epi 2 (2,A2) 

Substituting values of GArms Ggpj into (A2) yields G Am = 0.87 s. Similarly, for B and 
C quality measurements, we find Og m = 1.94 s and Gem = 2.86 s. From (Al), the total 
uncertainty for A, B, and C quality is, respectively, a A = 1.25 s, o B = 2.14 s, and g c = 
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3.00 s. 

In the weighted regression experiments the A, B, and C quality measurements are 
weighted inversely by their measurement variance. 
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appendix 2.B: Errors in the along- axis profiles and construction of the 

data Covariance matrix 


Errors in Bathymetry, Geoid, and Travel Time Profiles 

Uncertainties in the along-axis profiles of geoid, bathymetry, and travel times are 
important information in the inversion. The gridded bathymetric data [U .S. Naval 
Oceanographic Office, 1985] include corrections for the deviation of water column acoustic 
velocity from the assumed value of 1500 m r 1 . The geoid data, provided in the form of a 
0.25* x 0.25* grid [Marsh et al., 1986], include corrections for orbit errors, instrument and 
atmospheric propagation effects, and solid Earth and ocean tides. 

We have averaged the bathymetry and geoid height values within a 1* x 1* box 
centered at each SS bounce point. The averaging yields a representative value for a region 
over approximately one horizontal wavelength of the SS wave and acts to smooth out 
short-wavelength variations. Both bathymetry and geoid are corrected for subsidence with 
seafloor age, using the plate cooling model [Parsons and Sclater , 1977; Parsons and 
Richter, 1980]. Error introduced into depth and geoid anomalies by isochron mislocation 
is difficult to estimate precisely, but for an error in age of 2 My, depth and geoid errors at 
80 My would be about 30 m and 0.2 m, respectively, while at 2 My, an error in age of 2 
My would have a much larger affect, giving depth and geoid errors of 350 m and 0.3 m. 
The magnitude of this error highlights the importance of accurate age determination, 
especially at young ages. 

The presence of oceanic sediments is another source of error. In the Atlantic 
Ocean, the sediment thickness increases regularly from less than 100 m along the Mid- 
Atlantic Ridge toward continental margins where it can exceed 1 km [Ewing et al, 1973; 
Tucholke, 1986]. A 1-km sediment thickness leads to corrections to residual depth and 
geoid of about 500 m and 0.3 m, respectively [Cazenave et al., 1988; Sheehan and 
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McNutt, 1989]. On Atlantic lithosphere of 100 My age or less the sediment thickness is 
less than 500 m in most areas. Hence, neglecting the sediment loading correction should 
not be crucial in this region. 

The along-axis profile of SS-S residual is a weighted moving average of 10 
adjacent data points grouped by latitude, using the weights discussed in Appendix 2.A. 

The same weights and moving average are applied to geoid and bathymetry values at a 
given SS bounce point, even though bathymetry and geoid data are presumed to be of equal 
quality, in order that these profiles will be consistent with the SS-S residuals. The standard 
error of the mean values for SS-S residual ranges from 0.2 s to 1.6 s. For bathymetry the 
range of standard deviations from the mean value is from 24 to 370 m, and for geoid, 0.08 
to 1.0 m. The largest variances in the bathymetry and geoid data come from the Iceland 
region (north of 60*N), and may be due to the more complicated tectonics of this region 
[White, 1988]. 

Before Fourier transforming, the along-axis profiles must be interpolated to a 
constant spacing. We use a simple linear interpolation scheme to estimate values at a 0.5* 
spacing. We estimate that the typical error in the interpolated data is comparable to that in 
the along-axis moving averages, which for bathymetry is on the order of 125 m, for geoid 
0.4 m, and for travel time 1 s. 

Effect of Crustal Thickness Variations 

Our poor fit to topography in the inversion experiments can be at least partially 
attributed to unmodelled effects such as crustal thickness differences. Variations in oceanic 
crustal thickness about the typical value of 6-7 km [Spudich and Orcutt, 1980; White, 

1984; Purdy and Detrick, 1986] are generally thought to be small at horizontal scales of 
100 km and greater. However, the crust beneath the Azores plateau is estimated to be 
between 8 and 9 km thick [Searle, 1976; Whitmarsh etal., 1982] and that beneath Iceland 
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is at least 8 to 14 km thick [Bjornsson, 1983]. By simple isostatic mass balance, the depth 
anomaly due to excess crustal thickness in the Azores region would be about 400 m, and at 
Iceland, 200 m to 1.6 km. In general, simple variations in crustal thickness are insufficient 
to produce a significant SS-S residual. For crustal and mantle S wave velocities of 3.5 and 
4.4 km/s, a 2-km variation in crustal thickness would contribute less than 0.2 s to an SS-S 
differential travel time corrected for differences in bathymetry. However, at Iceland, where 
the crust is estimated to be as much as 14 km thick, the additional SS-S travel time could be 
up to 0.8 s. 


Data Covariance Matrix 

The data covariance matrix Rdd is of the form 


/ <*h(kl) 2 


Rdd ~ 


\ 


0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 


0 0 0 

0 0 0 

Oh(k„) 2 0 0 

0 o N (ki) 2 0 

0 0 .. 

0 0 0 

0 0 0 

0 0 0 

0 0 0 


0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

a N (k„) 2 0 0 

0 oAO 2 0 

0 0 .. 

0 0 0 


0 

0 

0 

0 

0 

0 

0 

0 

OtCkn ) 2 


where Oh 2 , Ofo 2 , and o t 2 are the nominal variances of the bathymetry, geoid, and 
travel time data, respectively. We may choose to construct the data covariance matrix not to 
reflect the true variance of the data but rather to allow weighting between the different data 
sets. In this way, the data covariance matrix can be altered to test the relative contributions 
of different data sets to the inversion results. 
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In ail of our inversions, the covariance matrix is constructed to weight the three 
data sets approximately equally. For example, examination of Figure 2. 1 1 indicates that at 
3000- Ion wavelength the amplitude erf" the geoid signal is approximately 4 m, bathymetry 1 
km, and travel time 2 s. Thus if a value of 1 m is chosen for On, then a value of 0.5 s for 
o t and 0.25 km for Oh should yield approximately equal weighting of data sets. The 
corresponding l/o 2 values are then 1 for geoid, 4 for travel time, and 16 for bathymetry. 
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TABLE 2.1. Digital Seismograph Stations Used 


Station Code 

Network 

Latitude (°N) 

Longitude (°E) 

ALQ 

DWWSSN 

34.942 

-106.458 

ANMO 

SRO 

34.946 

-106.457 

ANTO 

SRO 

39.869 

32.794 

BCAO 

SRO 

4.434 

18.535 

BER 

DWWSSN 

60.387 

5.326 

BOCO 

SRO 

4.587 

-74.043 

COL 

DWWSSN 

64.900 

-147.793 

GAC 

CAN 

45.70 

-75.47 

GDH 

DWWSSN 

69.250 

-53.533 

GRFO 

SRO 

49.692 

11.222 

JAS1 

DWWSSN 

37.947 

-120.438 

KBS 

DWWSSN 

78.917 

11.924 

KEV 

DWWSSN 

69.755 

27.007 

KOND 

ASRO 

59.649 

9.598 

LQN 

DWWSSN 

46.750 

-121.810 

NB04 

NARS 

52.810 

6.670 

NE06 

NARS 

50.100 

4.600 

NE09 

NARS 

44.850 

0.980 

NE10 

NARS 

43.090 

-0.700 

NE11 

NARS 

41.480 

-1.730 

NE12 

NARS 

40.640 

-4.160 

NE13 

NARS 

38.690 

-4.090 

NE14 

NARS 

37.190 

-3.600 

NE15 

NARS 

50.810 

5.780 

NE16 

NARS 

45.763 

3.103 

NE17 

NARS 

39.881 

-4.049 

RSCP 

RSTN 

35.600 

-85.569 

RSNT 

RSTN 

62.480 

-114.592 

RSNY 

RSTN 

44.548 

-74.530 

RSON 

RSTN 

50.859 

- 93.702 

RSSD 

RSTN 

44.120 

-104.036 

SCP 

DWWSSN 

40.795 

-77.865 

SSB 

GEOSCOPE 

45.280 

4.540 

TOL 

DWWSSN 

39.881 

-4.049 

WFM 

GEOSCOPE 

42.610 

-71.490 

ZOBO 

ASRO 

-16.270 

-68.125 


ASRO = Abbreviated Seismic Research Observatory Network 


CAN = Canadian Seismic Network 

DWWSSN = Digital World-Wide Standardized Seismograph Network 
GEOSCP = Geoscope Network 

NARS = Network of Autonomous Recording Seismographs 
RSTN = Regional Seismic Test Network 
SRO = Seismic Research Observatory Network 
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TABLE 2.2. Adopted Constants 


Variable 

Description 

Value 

a 

volumetric coefficient erf thermal 

2.5 x 10-5 K-l (•) 


expansion 


Po 

average mantle density 

3300 kg m' 3 

Pw 

density of seawater 

1000 kg m' 3 

r 

gravitational constant 

6.67 x 10* 11 N m 2 kg* 2 

g 

surface gravitational acceleration 

9.8 m s' 2 

dvs/9T 

thermal coefficient erf 

-6.0 x 10’ 4 km/s K' 1 


shear velocity 


dvj/SMg 

variation of shear velocity with Mg# 

1.8 x 10- 2 km/s/Mg# < b 

dp/9Mg 

variation of density with Mg# 

-12 kg/m^/Mg# (b) 

P 

average SS ray parameter at 70* 

0.1375 s/km 


(a) Stacey [ 1977], Duffy and Anderson [1989] 

(b) Akimoto [1972] 
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TABLE 2.3. Inversion Models and Variance Reduction 


Model: Temperature variations only Variance reduction, % 


Layer 

thickness 

Viscosity 

structure 





ss-s 

0-150 km 

cvm 

180 K 

53 

25 

79 

58 

0-150 km 

lvz 

230 K 

57 

27 

79 

66 

0-300 km 

cvm 

60 K 

47 

21 

85 

41 

0-300 km 

lvz 

110 K 

57 

24 

85 

65 

0-650 km 

cvm 

20 K 

41 

14 

91 

25 

0-650 km 

lvz 

33 K 

49 

17 

83 

51 

Model: Compositional variations only 

Variance reduction, % 

Laver 

Viscosity 











Mwwm 

ss-s 

A 


0-150 km 

cvm 

1.1 

33 

46 

74 

-y 


0-150 km 

lvz 

2.4 

44 

75 

76 

-9 


0-300 km 

cvm 

0.4 

33 

29 

87 

-6 


0-300 km 

lvz 

1.3 

43 

73 

73 

-7 


0-650 km 

cvm 

0.1 

32 

19 

93 

-3 


0-650 km 

lvz 

0.5 

49 

65 

86 

+5 


Model: Thermal and compositional 

variations in same layer 

Variance reduction, % 

Laver 

Viscosity 

AT 

AMg# 






structure 

range 

HRTTTTSHi 

Kona 



0-150 km 

cvm 

210 K 

1.5 

75 

44 

78 

100 

0-150 km 

lvz 

235 K 

2.1 

86 

75 

80 

100 

0-300 km 

cvm 

110 K 

0.7 

73 

28 

89 

100 

0-300 km 

lvz 

125 K 

1.1 

84 

74 

76 

100 

0-650 km 

cvm 

55 K 

0.4 

71 

18 

94 

100 

0-650 km 

lvz 

60 K 

0.8 

85 

66 

88 

100 


Model: Thermal inversion in layers as noted, compositional variations 0-50 km only 


Layer 

thickness 

Viscosity 

structure 

AT 

ranae 

0-150 km 

cvm 


0-150 km 

lvz 


0-300 km 

cvm 


0-300 km 

lvz 


0-650 km 

cvm 

25 K 

0-650 km 

lvz 

35 K 


AMg# 

ranee 

Variance reduction, % 
total bathvm eeoid SS-S 

KKflHi 

84 

83 

77 

91 

4.5 

85 

85 

70 

96 

5.9 

80 

84 

86 

72 

4.7 

86 


77 

89 


73 

85 

92 

47 

4.6 

75 

82 

73 

71 


evtn = constant viscosity mantle 
lvz = mantle with low viscosity zone 
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TABLE 2.4. Spectral Coefficients of Travel Time, Geoid, 
and Bathymetry Data and of Inversion Solutions 

Spectral Coefficients of Observations 



3552 3.25 -1.59 2.15 -1.08 2.73 -1.81 

2368 3.25 2.12 0.85 -2.99 1.42 -1.12 

1776 1.63 0.64 0.58 0.96 1.55 -3.00 

1421 0.52 -1.32 0.23 -1.70 0.84 -2.62 


Spectral Coefficients of Models (from Inversion) 

Inversion for temperature variations only 



2368 0.49 2.71 0.83 2.71 0.53 2.71 0.91 2.71 

1776 0.30 1.68 0.48 1.68 0.33 1.68 0.56 1.68 


1421 0.35 -2.13 0.54 -2.13 0.40 -2.13 0.68 -2.13 


Mantle with Low Viscosity Zone, 0-150 km 



3552 2.55 -1.59 2.25 -1.59 2.97 -1.59 3.84 -1.59 

2368 0.48 2.46 0.21 2.46 0.61 2.46 0.79 2.46 

1776 0.39 2.44 0.14 2.44 0.52 2.44 0.67 2.44 

1421 0.48 -2.26 0.21 -2.26 0.67 -2.26 0.86 -2.26 


Constant Viscosity Mantle, 0-300 km 
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Mantle with Low Viscosity Zone, 0-300 km 



Geoid 

amp phase 

H I 

KSralSI 

Temperature 
amp phase 



4.21 1.44 

1.62 1.44 

2.26 1.44 

3552 

2.42 -1.59 

2.27 -1.59 

3.02 -1.59 

4.23 -1.59 

2368 

0.15 1.94 

0.07 -1.20 

0.23 1.94 

0.33 1.94 

1776 

0.39 3.11 

0.39 -0.03 

0.74 3.11 

1.03 3.11 

1421 

0.30 -2.45 

0.27 0.69 

0.65 -2.45 

0.91 -2.45 


Constant Viscosity Mantle, 0-650 km 



Wavelength, 




1.58 

4 Si 133 


MTU MW f 

3552 

0.65 

-1.33 

3.15 -1.33 

0.68 -1.33 

3.13 -1.33 

2368 

0.23 

2.98 

0.93 2.98 

0.28 2.98 

1.27 2.98 

1776 

0.16 

1.59 

0.53 1.59 

0.21 1.59 

0.98 1.59 

1421 

0.20 

-2.16 

0.57 -2.16 

0.30 -2.16 

1.39 -2.16 


Mantle with Low Viscosity Zone, 0-650 km 



Wavelength, 

km 

Bathymetry 
amp phase 

Geoid 

amp phase 

ss 

amp phase 

Temperature 
amp phase 


KHI11 

4.57 1.56“ 



3552 

1.33 -1.46 

3.39 -1.46 

1.73 -1.46 

3.79 -1.46 

2368 

0.08 -3.02 

0.01 -3.02 

0.15 -3.02 

0.34 -3.02 

1776 

0.36 -3.09 

0.40 0.05 

0.99 -3.09 

2.16 -3.09 

1421 

0.20 -2.54 

0.27 0.61 

0.70 -2.54 

1.54 -2.54 


Inversion for compositional variations only 


Constant Viscosity Mantle, 0-150 km 


Bathymetry 
amp phase 

Geoid 
amp phase 


Composition 
amp phase 

wmurwmm 

1 .41 1.68 

2.74 1.68 

0.31 -1.47 

2.78 1.68 

3552 

1.48 -1.27 

2.69 -1.27 

0.33 1.87 

2.93 -1.27 

2368 

1.18 2.39 

2.01 2.39 

0.27 -0.75 

2.36 2.39 

1776 

0.83 0.66 

1.34 0.66 

0.19 -2.48 

1.69 0.66 

1421 

0.24 -1.19 

0.37 -1.19 

0.06 1.96 

0.50 -1.19 
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Mantle with Low Viscosity Zone, 0-150 km 


|m t wun m.i.' j 

Bathymetry 
amp phase 

Geoid 
amp phase 

53 

amp phase 

Composition 
amp phase 

kHfiDIHH 

HK!9K*1di 

Km 

1.67 

0.33 -1.47 

i.4i f.6? 

3552 

2.35 -1.35 

2.06 

-1.35 

0.56 1.79 

2.42 -1.35 

2368 

3.07 2.19 

1.31 

2.19 

0.80 -0.96 

3.45 2.19 

1776 

1.87 0.56 

0.66 

0.56 

0.52 -2.58 

2.23 0.56 

1421 

0.46 -0.90 

0.20 

-0.90 

0.13 2.24 

0.57 -0.90 


Constant Viscosity Mantle, 0-300 tan 



Bathymetry 
amp phase 

Geoid 
amp phase 

33 

amp phase 



” 1.04 1.68 

3.6 1 1.68 

"0.22-1.46 " 

3.15 1.68 

3552 

0.89 -1.22 

2.70 -1.22 

0.19 1.92 

2.76 -1.22 

2368 

0.64 2.51 

1.72 2.51 

0.14 -0.63 

2.05 2.51 

1776 

0.51 0.70 

1.24 0.70 

0.12 -2.44 

1.72 0.70 

1421 

0.16 -1.25 

0.36 -1.25 

0.04 1.89 

0.58 -1.25 


Mantle with Low Viscosity Zone, 0-300 km 


Wavelength, 

km 

Bathymetry 
amp phase 

Geoid 
amp phase 

33 

amp phase 

Composition 
amp phase 

■BirilHI 

warn 

3.33 1.68 

0.2o -1.46 

1.15 f .68 

3552 

2.25 -1.34 

2.10 -1.34 

0.58 1.80 

2.51 -1.34 

2368 

2.72 2.01 

1.37 -1.13 

0.88 -1.13 

3.84 2.01 

1776 

0.88 0.37 

0.89 -2.78 

0.35 -2.78 

1.51 0.37 

1421 

0.28 -0.22 

0.25 2.93 

0.12 2.93 

0.54 -0.22 


Constant Viscosity Mantle, 0-650 km 



Bathymetry 

— 

33 

Composition 

km 

amp phase 

amp phase 


amp phase 

7104 



0.13 -1.45 

195 1.69 

3552 

0.53 -1.18 

2.58 -1.18 

0.12 1.97 

2.61 -1.18 

2368 

0.37 2.63 

1.46 2.63 

0.09 -0.51 

2.04 2.63 

1776 

0.34 0.72 

1.14 0.72 

0.09 -2.42 

2.14 0.72 

1421 

0.12 -1.24 

0.34 -1.24 

0.04 1.90 

0.84 -1.24 


Mantle with Low Viscosity Zone 

, 0-650 km 


Wavelength, 

Bathymetry 

Geoid 

SS 

IMi&l 

km 

amp phase 

amp phase 


amp phase 


0.72 1.69 

3.96 1.6$ 

0.15 -1.45 

0.52 1.69 

3552 

1.04 -1.22 

2.67 -1.22 

0.28 1.92 

0.97 -1.22 

2368 

3.29 2.12 

0.38 2.12 

1.32 -1.02 

4.55 2.12 

1776 

0.84 0.32 

0.94 -2.83 

0.48 -2.83 

1.65 0.32 

1421 

0.24 0.13 

0.33 -3.01 

0.18 -3.01 

0.61 0.13 
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Inversion for both temperature and compositional variations 


HKSMI 

Bathymetry 
amp phase 

Geoid 
amp phase 

S3 

amp phase 


isition 

phase 



illtl IBM 



2.61 

-2.58 

3552 

1.65 -1.33 

3.01 -1.33 

2.73 -1.81 

2.88 -1.76 

1.59 

0.69 

2368 

1.10 2.42 

1.88 2.42 

1.42 -1.12 

1.12 -1.20 

2.89 

2.21 

1776 

0.74 0.75 

1.19 0.75 

1.55 -3.00 

1.33 -3.07 

2.56 

0.35 

1421 

0.27 -1.46 

0.41 -1.46 

0.84 -2.62 

0.81 -2.56 

0.86 

0.15 


Mantle 

with Low Viscosity Zone, 0-150 km 

/* " 



km 

amp phase 

amp phase 

amp phase 

amp phase amp phase 

— tub — 




3.59 1.05 

1.63 -2.63 

3552 

2.84 -1.42 

2.50 -1.42 

2.73 -1.81 

3.13 -1.74 

0.96 -0.49 

2368 

2.92 2.21 

1.25 2.21 

1.42 -1.12 

0.63 -1.32 

3.74 2.15 

1776 

1.62 0.67 

0.57 0.67 

1.55 -3.00 

1.10 3.09 

2.65 0.45 

1421 

0.52 -1.37 

0.22 -1.37 

0.84 -2.62 

0.83 -2.47 

0.67 -0.32 


Constant Viscosity Mantle, 0-300 km 
— «— i * ,ij * rs 



— 0 

km 

amp 

phase 

amp 

phase 

amp 

phase 


HD1 





K!££J 

3552 

0.95 

-1.26 

2.91 

-1.26 

2.73 

-1.81 

2368 

0.60 

2.55 

1.62 

2.55 

1.42 

-1.12 

1776 

0.46 

0.78 

1.12 

0.78 

1.55 

-3.00 

1421 

0.18 

-1.49 

0.40 

-1.49 

0.84 

-2.62 


amp phase amp phase 

3.fe 1 fm 3.65 -2.1T 


ft 
2.80 
1.26 
1.42 
0.83 


r. 

-1.78 

-1.17 

-3.04 

-2.58 


2.18 

2.26 

2.25 

0.87 


-Z 
1.05 
2.19 
0.30 
0.26 


Mantle with Low Viscosity Zone, 0-300 km 


TIeoId 




Temperature Composition 


— , 

km 

amp 

phase 

amp 

phase 

amp. 

phase 

IriDIBH 


“T 53” 


HTST - 

3.53 

K!££] 

3552 

2.75 

-1.41 

2.57 

-1.41 

2.73 

-1.81 

2368 

2.56 

2.01 

1.29 

-1.13 

1.42 

-1.12 

1776 

0.61 

0.50 

0.61 

-2.64 

1.55 

-3.00 

1421 

0.21 

-1.13 

0.19 

2.01 

0.84 

-2.62 


Sr 


ami 

T 
3.18 
0.55 
1.25 
0.80 


phase 

U5T 
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Mantle with Low Viscosity Zone, 0-650 km 
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Mantle with Low Viscosity Zone, 0-300 km 
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0-150 km, 0-300 km, 0-650 km = depths over which temperature or composition is 

allowed to vary , 

amp = relative spectral amplitude. The amplitudes of both observed and predicted data 
profiles are normalized by dividing by the average spectral amplitude of the observed 
profile. The amplitudes of model parameters are normalized by dividing by the average 
spectral amplitude of the parameter profile, 
phase = phase in radians 
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Figure 2. 1 . 


Figure 2.2. 


Figure 2.3. 


Figure 2.4. 


Figure Captions 


Schematic paths of S and SS phases. For the range of epicentral distances 
considered in this study the bottoming depth for S ranges from 1450 to 
2370 km, and for SS from 670 to 900 km. 

An example of the measurement of SS-S differential travel time for the 
event of December 24, 1985 (10 km focal depth), at GDH (63* epicentral 
distance), (a) “Synthetic” SS pulse generated from S. The S pulse is 
windowed and attenuated to account for the additional time that SS travels 
in the mantle (t* * 3 s), and a tc/ 2 phase shift is applied, (b) Windowed SS 
wave pulse, (c) Cross-correlation of the trace in (b) with that in (a). The 
differential travel time residual is -5.04 s. 

Distribution of earthquakes (triangles) and seismograph stations (circles) 
used to measure SS-S differential travel times. Stations are from the 
GDSN, NARS, and GEOSCOPE digital arrays. Earthquakes are from the 
Harvard CMT catalogue (generally mb > 5.0) from the years 1977- 1987. 
Lambert equal area projection with pole of projection at 45*N, 40*W. 

(a) Map view of SS-S residuals relative to PREM [Dziewonski and 
Anderson , 1981], corrected for Earth ellipticity and seafloor bathymetry. 
Residuals are plotted at the SS bounce point The size of each symbol 
scales linearly with magnitude of the residual. Lambert equal area 
projection with pole of projection at 40*N, 60*W. Negative residuals 
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indicate either early SS or late S. Plate boundaries are from DeMets et al. 
[1990]. 

(b) Same as (a) but including data only for SS bouncepoints on lithosphere 
younger than 100 My. 

Figure 2.5. SS-S travel time residual versus square root of seafloor age for data from 
0-60*N. Each plotted point represents the weighted mean of 14 adjacent 
data points. Weights are constructed from variances determined as 
discussed in Appendix A. Horizontal and vertical bars are standard errors 
of the means of the travel time residuals and (age) 1 / 2 . Linear regression 
yields a slope of -0.68 ± 0.08 s/ (My) 1 / 2 for a 0-100 My age range (solid 
line) or -0.76 ± 0.09 s/ (My) 1 / 2 for a 0-80 My range (dashed line). 

Figure 2.6. Age-corrected SS-S residual (see text) versus azimuth 0. Each plotted 

point represents the weighted mean of 10 adjacent data points. The solid 
curve shows the best-fitting 40 variation derived from these data. The 
dashed curve shows the preferred model of Kuo et al. [1987], which 
corresponds to an alignment of the a axis of olivine in the approximate 
direction N13°W. 

Figure 2.7. (a) Map view of age-corrected SS-S residuals. 

(b) Same as (a) but including data only with SS bounce points on 
lithosphere younger than 100 My. 


Figure 2.8. 


Map view of the distribution of sampling azimuths. Lines indicate the 
wave path azimuth at the SS bounce point Mercator projection. 
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Figure 2.9. Age-corrected SS-S residual versus latitude along the Mid-Atlantic Ridge 
from 10 to 90*N. The residuals shown are moving averages (such that 
each point is used twice) of 12 adjacent data points. Bounce points on 
lithosphere of age 0-100 My are used. The approximate locations of 
several fracture zones (Fifteen-Twenty, Kane, Atlantis, Oceanographer, 
and Charlie-Gibbs, denoted by abbreviations) and of the Iceland and 
Azores hotspots are indicated. 

Figure 2. 10. Linear correlation, by highest harmonic degree removed from the geoid, of 
observed SS-S residual with geoid height measured at the corresponding 
SS bounce point. Both travel time and geoid residuals are age-corrected. 
First the raw [Marsh et al., 1986] geoid data are correlated with SS-S 
residuals and a slope and correlation coefficient determined. Then a geoid 
reference field [Lerch et al., 1979] up to degree and order 2 (with taper to 
degree and order 6) is calculated and removed from the geoid data, the 
slope and correlation coefficient with SS-S calculated, and so on for 
higher harmonic degrees C removed from the geoid data, with appropriate 
tapers (up to f+ 4). (a) Linear correlation coefficient between geoid and 
SS-S residuals vs. highest harmonic degree and order removed from the 
geoid. (b) Slope of the correlation between geoid and SS-S data, as a 
function of highest harmonic degree and order removed from the geoid. 
Extra points at degree and order 10 are obtained by using different tapers 
(no taper, taper to C = 14, and taper to l- 15). 


Figure 2.11. Comparative plots of age-corrected (a) bathymetry, (b) geoid, and (c) SS-S 
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residual along the Mid-Atlantic Ridge, 10-65’N. Bathymetry and geoid 
have been high-pass filtered (see text). All of the residuals shown are 
moving averages of 10 adjacent data points. Bounce points from 
lithosphere of age 0-100 My are used, except that data from the Labrador 
Sea region are omitted. 

Figure 2. 1 2. Upper mantle portion of the kernels for geoid and topography at two 

wavelengths X = 2n/k few two viscosity models. The convecting region in 
both models is overlain by a high-viscosity layer 40 km thick, with 
viscosity 10 4 that of the underlying mantle. 

(a) High-viscosity lid is underlain by a mantle of uniform viscosity and 
other physical parameters. 

(b) High-viscosity lid is underlain by a zone extending to a depth of 200 
km having a viscosity equal to 0.01 that of the underlying mantle. 

Figure 2. 1 3. Upper mande portion of the kernels for geoid and topography at two 

wavelengths in spherical versus cartesian coordinates. The convecting 
region is overlain by a high-viscosity layer 40 km thick in each of two 
models for viscosity structure. 

(a) Underlying mande is of uniform viscosity. Cartesian kernels are for 
4000 km wavelength (solid lines) and spherical kernels are for (= 10 
(dashed lines). 

(b) High-viscosity lid is underlain by zone extending to 200 km depth 
having a viscosity equal to 0.01 that of the underlying mande. Cartesian 
kernels are for 4000 km wavelength and spherical kernels are for C- 10. 

(c) Same as (a) but with Cartesian kernels for 6667 km wavelength and 
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spherical kernels for (=6. 

(d) Same as (b) but with Cartesian kernels for 6667 km wavelength and 
spherical kernels for C= 6. 

Figure 2. 14. Results of combined inversion of geoid, bathymetry, and SS-S travel time 
residuals for upper mantle temperature variations. The viscosity structure 
is taken to consist of a 40-km-thick high-viscosity lid overlying a constant- 
viscosity halfspace. 

(a) Three solutions for along-axis temperature variations: Dotted line: 
Temperature perturbations constrained to be uniform over 0-150 km 
depth. Long-dashed line: Temperature perturbations constrained to be 
uniform over 0-300 km depth. Short-dashed line: Temperature 
perturbations constrained to be uniform over 0-650 km depth. 

(b) Observed (solid line) and predicted along-axis profiles of SS-S travel 
time residual. The “observed” profile is actually a filtered version of the 
observations, containing only the wavelengths used in the inversion (1400 
to 7100 km). Predicted profiles were calculated from equation 5. Line 
types correspond to those of the temperature models. 

(c) Observed and predicted along-axis geoid profiles. Same treatment as 
in (b). 

(d) Observed and predicted along-axis bathymetry profiles. Same 
treatment as in (b). 

Figure 2.15. Same as Figure 2.14 except for that the viscosity structure includes a zone 
extending from the base of the lid to a depth of 200 km with a viscosity 
equal to 0.01 that of the underlying mantle. 
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Figure 2. 16. Results of combined inversion of geoid, bathymetry, and SS-S travel time 
residuals for variations in upper mantle composition (Mg#). The viscosity 
structure is taken to consist of a 40-km-thick high-viscosity lid overlying a 
constant-viscosity halfspace. 

(a) Three solutions for along-axis composition variations: Dotted line: 
Composition perturbations constrained to be uniform over 0-150 km 
depth. Long-dashed line: Composition perturbations constrained to be 
uniform over 0-300 km depth. Short-dashed line: Composition 
perturbations constrained to be uniform over 0-650 km depth. 

(b) Observed (solid line) and predicted along-axis profiles of SS-S travel 
time residual. 

(c) Observed and predicted along-axis geoid profiles. 

(d) Observed and predicted along-axis bathymetry profiles. 

Figure 2.17. Same as Figure 2. 1 6 except for that the viscosity structure includes a zone 
extending from the base of the lid to a depth of 200 km with a viscosity 
equal to 0.01 that of the underlying mantle. 

Figure 2.18. Results of combined inversion of geoid, bathymetry, and SS-S travel time 
residuals for both upper mantle temperature and composition variations. 
The viscosity structure is taken to consist of a 40-km-thick high-viscosity 
lid overlying a constant-viscosity halfspace. 

(a) Three solutions for along-axis temperature variations: Dotted line: 
Composition perturbations constrained to be uniform over 0-150 km 
depth. Long-dashed line: Composition perturbations constrained to be 
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uniform over 0-300 km depth. Short-dashed line: Composition 
perturbations constrained to be uniform over 0-650 km depth. 

(b) Three solutions for along-axis composition variations: Dotted line: 
Composition perturbations constrained to be uniform over 0- 150 km 
depth. Long-dashed line: Composition perturbations constrained to be 
uniform over 0-300 km depth. Short-dashed line: Composition 
perturbations constrained to be uniform over 0-650 km depth. 

(c) Observed (solid line) and predicted along-axis profiles of SS-S travel 
time residual. 

(d) Observed and predicted along-axis geoid profiles. 

(e) Observed and predicted along-axis bathymetry profiles. 

Figure 2. 19. Same as Figure 2. 18 except for that the viscosity structure includes a zone 
extending from the base of the lid to a depth of 200 km with a viscosity 
equal to 0.01 that of the underlying mantle. 

Figure 2.20. Same as Figure 2.18 but compositional variations constrained to be from 
0-50 km depth only. 

Figure 2.21. Same as Figure 2. 19 but compositional variations constrained to be from 


0-50 km depth only. 
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Figure 2.3 
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Figure 2.4a 
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Figure 2.4b 






Figure 2.6 
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Figure 2.7a 




Figure 2.7b 
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Figure 2.8 
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Geoid/SS Correlation Coefficient* 



Geoid/SS slope 



Figure 2.10 
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Chapter 3 


Differential Shear Wave Attenuation and Its Lateral Variation 

in the North Atlantic Region 

Introduction 

The rate of attenuation of seismic waves provides important information comple- 
mentary to that provided by seismic velocities. Seismic attenuation can be strongly 
affected by variations in the physical state and temperature field of the Earth's interior. 
Recent analysis of global variations in Q structure using long-period surface waves 
suggests that the most significant variations are confined to the upper mantle low-velocity, 
low-Q zone [Ritzwoller et al., 1989] and that the main contributors to these variations are 
the mid-ocean ridges (M.H. Ritzwoller, pers. comm., 1989). Earlier studies with body 
waves suggest that the upper mantle beneath mid-ocean ridges should display a 
significantly greater than average level of S-wave absorption [e.g., Molnar and Oliver , 
1969; Solomon , 1973]. A more detailed analysis of mantle attenuation beneath ridges and 
its vari atio n with lithospheric age would do much to sharpen the inferences from global 
models and to assess whether the observed lateral variations reflect differences in the 
percentage of partial melt [Ritzwoller et al., 1989] or only in temperature. 

There is considerable difficulty in relating laboratory experiments on seismic 
properties to propagation characteristics in the mantle because of differences in pressure 
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and temperature conditions beetween the two environments. Moreover, most laboratory 
experiments are conducted at frequencies much greater than those in the seismic wave 
band. Direct quantitative measurement of variations in Q in the Earth are therefore 
important In this paper we present measurements of differential shear wave attenuation in 
the north Atlantic and its variation with lithosphere age, and we utilize thermal models to 
interpret the observations in terms of an empirical relation between differential attenuation 
and temperature. Our approach is to measure the spectral amplitude ratio of long-period 
SS and S waves and to attribute variations in this ratio to differential attenuation of the SS 
waves in the upper mantle near the bounce points of these surface-reflected phases. This 
work is a complement to our study of SS-S differential travel times in the north Atlantic 
[Sheehan and Solomon, 1991; Chapter 2]. Sip kin and Jordan [1980b] and Revenaugh 
and Jordan [1987] used the ScS phase as an effective tool to map lateral variations in 
QScS> which measures the vertically-averaged Q of the entire mantle. The SS-S differential 
attenuation measurement does not sample the significant lateral heterogeneity at the base of 
the mantle [Dziewonski, 1984] as does multiple ScS. Such measurements thus provide an 
important complement to ScS studies and can serve to isolate the upper mantle contribution 
to Q. The SS path segment in the upper mantle is not as nearly vertical as ScS, but many 
more paths are possible for SS because shallow sources can be employed, in contrast to 
ScS n , for which deep sources are normally required to avoid interference from surface 
waves. The most important result of this study is the documentation of an increase in Qs 
in the upper mantle with increasing plate age. We also examine along-axis variations in 
differential attenuation, and we test whether they might be produced by along-axis 
variations in temperature, such as those derived from the inversion of along-axis variations 
in SS-S travel times, geoid, and bathymetry [Sheehan and Solomon, 1991; Chapter 2]. 
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Measurement of differential attenuation 

The seismic data used in this study consist of long-pen od S and SS phases recorded 
at digital stations in the Global Digital Seismic Network (GDSN) [ Peterson et al., 1976; 
Peterson and Hutt, 1982], the Network of Autonomously Recording Seismographs 
(NARS), [Nolet and Vlaar, 1982], and the GEOSCOPE network [Romanowicz et al ., 
1984]. We use only transversely polarized (SH) seismograms (rotated from N-S and E-W 
components) to avoid interference from the SKS phase and contamination from P-SV 
conversions at the base of the crust and other near-surface discontinuities. 

A spectral ratio method is utilized to determine the differential attenuation between 
the phases SS and S [e.g., Teng, 1968; Schlue , 1981]. The spectral amplitude A of an S- 
phase can be written as the product of terms for source Ao(o>,0,<J>), instrument Aj(co), 
crustal layering at the source Acs(to) and receiver Ac, (to), and geometrical spreading G and 
attenuation t* along the path, as 

A(g>) = Ao(o),e,4>) Ai(co) Acs(cd) Ac,(to) G exp (-K ft*) (3.1) 

where 


t* 



ds 

Q(s)V(s) 


(3.2) 


and Q and V are the quality factor and wave velocity along the path s. For a given source 
and station, the SS phase will have the same values of Ao, Aj, A^, Ac as the S phase, so 
taking the natural logarithm of the amplitude ratio. 
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MAss/As) = -* f (t^s -$ = - 71 f5t* (3 3) 

where the differential attenuation 5t* = t*ss * t*s- Thus an estimate of 5t* may be 
obtained from the negative of the slope of the log of the amplitude ratio. 

In practice, the digital long-period seismograms are rotated, and the transverse (SH) 
component is windowed to isolate the S and SS pulses (Figure 3. la). Constant window 
lengths of 100 s are used for both the S and SS phases. A Kaiser-Bessel [Harris, 1978] 
taper with a = 3 (the parameter a controls the width of the central window versus the 
sidelobe amplitude) is applied to each windowed pulse. This taper is an effective tool for 
reducing amplitudes at the front and the tail of the window and acts to reduce the effects of 
signals not associated with the phase of interest We find that it is crucial to form narrow 
windows around the S and SS pulses, as the differential attenuation values can vary by as 
much as 1 s depending on how S and SS are windowed. The amplitude spectrum of the 
isolated phase is then obtained for the frequency band 0.01 to 0.15 Hz. The log of the 
amplitude spectrum is calculated and then smoothed by taking running averages over a 
0.04- Hz band (Figures 3.1b-c). An amplitude ratio spectrum (Figure 3. Id) is formed by 
subtracting the log spectrum of the S phase (Figure 3.1b) from that of the SS pulse 
(Figure 3.1c). The log amplitude ratio spectrum is smoothed by a moving average (0.02 
Hz wide) before a linear tit is derived. In general the amplitude ratios decrease 
systematically with frequency over the approximate band 0.01 to 0.08 Hz but above a 
frequency between 0.07 and 0. 1 1 Hz the amplitude ratios tend to increase. We attribute 
this increase to noise in both the S and SS phases. We measure the slope of the log 
amplitude ratio spectrum between a constant lower frequency of 0.016 Hz and three 
different upper frequencies: 0.08, 0.095, and 0.1 1 Hz. The upper frequency cutoff that 
gives the most negative (steepest) slope is that used to calculate 5t*. Choosing the steepest 
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slope may bias us toward large values of 8t* but our motivation is to obtain the most 
robust measure of attenuation and to avoid the part of the spectrum where the slope is 
likely to be susceptible to noise. Uncertainties involved with the measurement of St* are 
discussed further in Appendix 3 .A. 

As an alternative technique for measuring differential attenuation we have also 
constructed stacks of amplitude ratio spectra for groups of SS-S pairs. Stacks are 
constructed by calculating a mean value of In (Ass/As) with corresponding standard error 
at each frequency. A single slope is measured for the stacked spectrum, over the 
frequency band 0.016 to 0. 10 Hz. If we were to retain the use of variable cutoffs with the 
spectra going into the stacks the high frequency values of In (Ass/As) would not have a 
contribution from each seismogram and would therefore not be representative of all of the 
seismograms included in the stack. Since we choose a single cutoff frequency of 0. 10 Hz 
instead of picking the steepest slope, the values of St* from the stacks are consistendy 
lower than those measured individually. However, we find that the trends observed using 
stacked data are nearly identical to those obtained from the individual measurements. The 
analysis presented in the remainder of the paper is based on individual measurements. 


DATA 


The north Atlantic region has a good distribution of events and stations at suitable 
epicentral distances. Source-receiver separations were restricted to lie in the range 55* to 
86* to ensure clear separation of S and ScS at greater distances and to avoid triplication in 
SS at shorter distances. The SS and S phases bottom in the mantle between about 670 km 
and 2300 km depth. Digital seismograms selected for spectral analysis are from the travel 
time study of Sheehan and Solomon [1991]; see also Chapter 2. The distribution of 
sources and stations used to measure SS-S attenuation is shown in Figure 3.2. The 
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majority of data used in this study conies from equatorial fracture zone earthquakes 
recorded at North American and European stations, north and central Atlantic events 
recorded at North American stations, Central American events recorded at European 
stations, and Mediterranean and European earthquakes recorded in North America. The 
data set consists of over 150 6c* values with SS bounce points in the north Atlantic 
(Figure 3.3). A tabulation of all 5t* values is given in Appendix II. 

RESULTS 

Contributions to 5t* can arise from any portion of the wave path. It is expected, 
however, that 8t* will be controlled primarily by the upper mantle, where values of Q" 1 
are known to be large and variable [Solomon, 1972; Sipkin and Jordan , 1980b] and where 
seismic velocity and other physical properties also show significant lateral variation. We 
interpret the variations in 5t* in terms of lateral variation in Q within the crust and upper 
mantle beneath the surface reflection points of the SS wave path. The validity of this 
assumption is supported by the correlation, discussed further below, of 5t* with surface 
tectonic features (and with sediment thickness) in the vicinity of the SS bounce point In 
addition, the 5t* values correlate with the SS-S travel time residuals (Figure 3.4), and 
regression experiments by Kuo et al. [1987] and Woodward and Masters [1991] have 
shown that the SS-S travel time residuals can be attributed to lateral variations in upper 
mantle structure in the vicinity of the SS bounce point The 8t* values are further 
interpreted in terms of such upper mantle processes as lithospheric aging and along-axis 
heterogeneity in mantle structure. 
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Epicentral Distance Dependence 

The value of 8t* increases with epicentral distance A (Figure 3.5). Presumably this 
effect is at least partly due to the S phase being increasingly less attenuated relative to SS 
as a greater fraction of its path is in the high-Q deep mantle [Brune, 1977]. A simple 
correction for epicentral distance A was made by fitting a straight line, 5t* = (0.054 ± 
0.019) A - 0.265, to the observed values of 5t* versus A (Figure 3.5). The slope of this 
line is in approximate agreement with the predictions of several radial Q models [Anderson 
and Hart, 1978; Anderson and Given, 1982] (Figure 3.6). The form of the correction for 
epicentral distance dependence applied to the individual measurements is 

5t*corr = 5t*old - (b A + a) + 5t*av (3-4) 

where b = 0.054, a = -0.265, and 5t* av = 3.67 is the mean of all values of 5t*. 
Application of this correction does not substantially affect the relation between 8t* and 
travel times (Figure 3.7). 

Lithospheric Aging 

Cooling and thickening of the lithosphere should yield a tendency toward decreasing 
seismic attenuation with increasing lithospheric age. A linear regression experiment was 
p erfo r med to examine the correlation of 8t* with seafloor age. To obtain a representative 
age value for the region spanning approximately one horizontal wavelength of the incident 
(SS) wave, an average seafloor age was estimated for a 1* by 1* box centered on each SS 
surface bounce point [Sheehan and Solomon, 1991]. Measurements from zones with 
sediment thickness in excess of 1 km were excluded from die final age regression to avoid 
potentially large contributions from the low-Q sediments. Although the SS wave samples 
the upper mantle over a range of ages, we expect that the differential attenuation 
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contributed by the SS path segments on the younger and older sides of the bounce point 
approximately cancel so that the age at the SS surface bounce point is representative. 

The 8t* values for the north Atlantic are consistent with the expectation of a 
decrease in attenuation with increasing seafloor age. The slope derived by linear 
regression of 8t* with square root of age is -0.20 ± 0.07 s My**^ from 0 to 100 My, with 
a linear correlation coefficient of -0.97 (Figure 3.8a). This slope is consistent with 
observations of increasing Qs with lithosphere age [Cartas and Mitchell , 1981; Sipkin and 
Jordan , 1980b; Revenaugh and Jordan, 1987]. The trend of decreasing attenuation with 
increasing seafloor age is in the same sense as that for SS-S travel times for this region 
[ Sheehan and Solomon, 1991] (Figure 3.8b). 

To look for other systematic variations in the differential attenuation 
measurements, we correct for age t by removing the linear relation shown by the dashed 
line in Figure 3.8a. The form of the age correction for the individual measurements is 

5t*corr = 8t* 0 id - (b t ^ + a) + 8t* av (3.5) 

where b = -0.20, a = 4.71, and 8t* av = 3.42 is the average of all data after application of 
equation 3.4. 

Spatial Patterns of Age-corrected St* 

After removal of the age-dependence of 8t*, we may then search for systematic 
variations along the Mid- Atlantic Ridge. Such variations are seen in SS-S differential 
travel times and arc attributable to along-axis differences in upper mantle temperature and 
composition [Sheehan and Solomon, 1991; Chapter 2]. A north-south profile of 8t* 
values, constructed with averages of 10 adjacent data points grouped by latitude, is shown 
in Figure 3.9. In order to maintain consistency with our earlier travel time study [Sheehan 
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and Solomon , 1991] and to avoid possible biases associated with deviations from the 
simple plate cooling model, the along- axis 5t* profile is constructed using data with 
bounce points on lithosphere less than 100 My in age. Systematic variations of 8t* are 
evident with latitude, i.e., along the direction of the Mid-Atlantic Ridge axis. The along- 
axis variations show a variety of scales, notably at wavelengths of about 1000 - 2000 km 
in the region from about 15 to 35*N, and at about 6000 km wavelength from large 8t* 
(high attenuation, low Q) in the south (20-35’N) to lesser 5t* (low attenuation, high Q) 
farther north (45-60*N). These variations are qualitatively s i m ila r to the along-axis pattern 
of travel time residuals [Sheehan and Solomon , 1991; Chapter 2] and are also plausibly 
the result of thermal or compositional variations along the axis of the Mid- Atlantic Ridge. 

VELOCITY AND Q AS FUNCTIONS OF LITHOSPHERIC TEMPERATURE 


Because of the strong temperature dependence of both shear wave velocity and 
attenuation, the most straightforward hypothesis to explain the variations of SS-S travel 
time residual and differential attenuation with age is that both variations are due entirely to 
lateral variation in the thermal structure of the lithosphere, i.e., cooling of the oceanic 
plate. We have tested this hypothesis against the observed SS-S travel time delays as a 
function of age by means of the plate cooling model of Parsons and Sclater [1977]. 
Geotherms for different ages, according to this model, are shown in Figure 3. 10. 
Temperature perturbations AT can be converted to a seismic velocity perturbation by 
assuming that velocity varies linearly with AT and adopting a value for the partial 
derivative of shear wave velocity with respect to temperature, 8V s/3T. The resulting two- 
way travel time perturbation, as a function of horizontal wave number k, is given by 



(3.6) 


dv s f 

At(k) = 2 1rrJ 


AT(k,z) dz 


V s (z) 2 (l-p 2 V s (z) 2 ) 


1/2 


where V s (z) is from the reference shear velocity model [Dziewonski and Anderson, 1981], 
z is depth, p is the ray parameter, generally taken to be a representative value (0. 1375 
s/km) for the range of epicentral distances considered here, and Zmin and Zmax are the 
upper and lower boundaries of the region of significant lateral variations in temperature. 
We choose 50 My as a reference age, so that temperature anomalies are obtained by 
subtracting the geotherm at 50 My from the geotherm at an arbitrary age t : AT t (z) = T t 
(z) - T 50 (z). We use a value of -0.6 m s * 1 K * 1 for dVs/9T. This value is a factor of 1.5 
higher than the experimental values of Anderson et al. [1968] and Kwnazawa and 
Anderson [1969] for olivine at standard temperature and pressure but is comparable to the 
value of -0.62 m s' 1 K' 1 determined by McNutt and Judge [1990] by a least squares fit of 
Love- wave phase velocities to a temperature model for the lithosphere. 

The SS-S travel time delays predicted by the plate cooling model (Figure 3. 1 1 ) are 
in excellent agreement with observed values (Figure 3.8b). The slope derived by linear 
regression of observed SS-S residual with square root of age is -0.68 ± 0.08 s My 1 / 2 for 
bounce points between 0* and 60*N latitude and for ages between 0 and 100 My, with a 
linear correlation coefficient of -0.85 (Figure 3.8b). The value obtained using only that 
subset of the travel time data in common with the differential attenuation measurements is 
-0.74 ± 0. 10 s My 1 / 2 . For comparison, the slope calculated from the plate cooling model 
for 0-100 My age is -0.64 ± 0.01 s My 1 / 2 . The trend of the travel time residual versus 
age flattens out at about 80 to 100 My for both the observed and predicted residuals. 

These results support the hypothesis that most, if not all, of the dependence of travel time 
residual on plate age is due to plate cooling and does not require an additional contribution 
from below the lithosphere. 
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Wc next compare the age dependence of the observed St* values with lithospheric 
temperatures predicted by the plate cooling model. Specifically, we infer possible Q- 1 - 
temperature relations through the simultaneous use of the observed 5t*-age relation and 
geotherms as functions of lithospheric age as predicted by Parsons and Sclater [1977]. 
Assuming that the differential attenuation 5t* arises solely from the upper mantle portion 
of the SS path, we can express St* as 

St* = 2 f — 7 3 — i « 

(3.7) 

where Qs and Vs are the shear wave Q and velocity along the path, respectively, and p is 
the ray parameter. 

We choose to parameterize Q further in terms of temperature, following the 
Arrhenius law 

Q-l=A*lexp[-E/(RT(z))] (3-8) 

where the constant E is an activation energy and R is the gas constant. Substituting this 
relation into equation 3.7, we obtain 

8 “ = iAc“ T «V s (zH.-p^(z))‘« 

Thus an estimate of St* can be obtained by integration, given a velocity model, a 
geotherm, and the constants A and E. From t* and T(z) at a number of lithospheric ages, 
we can formulate an inversion for the parameters A and E by taking the derivatives of St* 
with respect to the model parameters and solving for perturbations to initial estimates of A 
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and E that minimize the difference between observed and calculated values of 8t*. The use 
of a constant activation energy might not be appropriate, as there may be several processes 
acting to produce the observed attenuation, each with its own activation energy. Despite 
these drawbacks, it is useful to try to parameterize Q in this way as it allows us to test 
whether or not a single thermally activated process can describe the observations. 

Since differential travel time versus age is well modelled if temperature variations 
are confined to the lithosphere, we solve for Q as a function of temperature in the 
lithosphere only. We assume a constant value of Q for the asthenosphere; specifically we 
adopt a value of 1 30 for Q in the depth range 125 - 500 km, which is the average value for 
this region given in the PREM Q model [Dziewonski and Anderson, 1981]. This 
asthenospheric contribution to 8t* obtained from integrating equation 3.7 over the depth 
range 125 to 500 km is approximately 1.7 s. 

The inversion for constants A and E indicates that equation 3.8 and the 
assumption that lateral temperature validations are confined to the lithosphere provide a 
good fit to observed values of 8t* (Figure 3. 12). The best fitting values of A and E are 
1.67 and 35 kJ mol' 1 , respectively. 

We also tested whether the laboratory-derived Q-temperature relation of Sato and 
Sacks [1989] also matches the observed 8t* values. From equations (1) and (2) of Sato 
and Sacks [1989] and assuming Qp/Qs = 2.25, we obtain 


v-l 


Q. 1 = 2.25 [ 3.5 + 


P(z) 

0.073 


] exp { -g [ 


T„(z) 

T(z) 


-a] } 


(3.10) 


where P(z) is pressure as a function of depth z, T m is the solidus temperature, T(z)/T m (z) 
is the homologous temperature, and g and a are piecewise constant functions of T/T m . We 
substitute this relation into equation 3.7 and integrate over depth to estimate the 
lithospheric contribution to 8t* (Figure 3.12). Including the asthenospheric contribution 
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to t* would add a constant of about 1.7 s to the plotted St* values. The variation of t 
versus a gc predicted by equation 3. 10 is much greater than that observed. 

We also examined several other parameter! zations of Qs' 1 in an attempt to assess 
the importance of pressure and homologous temperature dependence. A best fit for the 
parameterization Q -1 = A* 1 exp [-c Tm(z)/T (z)] produced constants of 2.86 for A and 2.05 
for c, and an equally good fit to the observed t* values as that obtained using equation 3.8. 
The parameterization Q- 1 = [A + b P(z)]* 1 exp [-c T m (z/T (z)] proved to be more 
problematic, as we were inverting for three parameters (A, b, and c) with a limited amount 
of data, and strong tradeoffs between these different parameters contribute to a 
nonuniqueness in the solution. We were able to obtain good fits to the data using widely 
different values of A, b, and c, and we had difficulty converging to physically reasonable 
solutions. For example, the values 2,1, and 1.7 for A, b, and c produce an rms misfit of 
0.12 s, whereas the values 3.5 x 10*2, -7.7 x 10*^, and 7.61 for A, b, and c produce an 
rms misfit of only 0.08 s. We would not expect a negative pressure dependence on Q as 
this last relation predicts, as Q generally increases with depth in the Earth. Some form of 
pressure dependence is warranted to avoid Q values that are too low at depth. A 
parameterization in terms of an activation volume will be attempted in a future study. 

Predicted distributions of Qs versus depth obtained using the different Q- 
temperature parameterizatxons are shown in Figure 3.13 for two different lithospheric 
ages. The Saxo and Sacks [ 1989] relation predicts Q values which are in general higher 
(except at very young ages) than the parameterizations more consistent with our data. The 
three Qs profiles shown as dashed lines are all nearly equally consistent with our St* 


observations versus age. 
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Discussion 

The differences between the Q* 1 -temperature relations we derived for oceanic 
lithosphere and those reported from laboratory measurements warrant discussion. The 
Q -1 -temperature relation of Sato et al. [1989] was obtained from laboratory measurements 
of P wave attenuation in peridotite at high temperature and pressure (but ultrasonic 
frequencies). Sato and Sacks [1989] argue that relationships between Q and temperature 
at ultrasonic frequencies can be extrapolated to the seismic band on the basis of the 
observed frequency independence of Q at ultrasonic frequencies and the prediction of 
reasonable mantle temperatures with their model. Their model actually tends to 
underpredict mantle temperatures, and to account for this difference they invoke a relation 
Qp'/Qp = 1.7 to 2.5 where Qp’ and Qp are quality factors of compressional waves for 
seismic and laboratory studies, respectively. Thus, they suggest that seismic Q will be a 
factor of about 2 greater than laboratory Q. We find that our seismically derived Q' 1 - 
temperature relation predicts smaller relative variations in upper mantle Q"* resulting from 
a given temperature difference than predicted by the experimentally obtained relations. 
Many factors could contribute to these discrepancies, including differences in frequency, 
pressure, temperature, and mineralogy. It is likely that different mechanisms for Q operate 
under these different conditions, each process having its own activation energy. Possible 
mechanisms for Q include partial melting, viscous grain boundary migration, and 
dislocation motion [e.g., Guegen et al ., 1981]. The observation that laboratory Q scales 
with homologous temperature [Sato and Sacks, 1989] is purely empirical. We have found 
that the St* values predicted by the Sato and Sacks [ 1989] relation can vary significantly 
upon application of slightly different solidus temperatures. For example, we found that 
the St* versus age relation predicted from the Sato and Sacks [1989] relation and 
employing the Wyllie [1971] dry peridotite solidus differed significantly from that 
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predicted "«« e the Takahashi [1986] dry peridotite solidus. This effect is most 
pronounced at temperatures near to and exceeding the solidus. 

We have not included the possibility of a small degree of melt fraction on 8t*. The 
influence of partial melt on seismic attenuation depends on the melt fraction and the 
geometry of its distribution [Cooper and Kohlstedt, 1986]. Melt affects velocities through 
a direct effect on the elastic moduli and a dispersive effect of relaxation or attenuation. The 
volume fraction and geometrical distribution of any melt in the mantle, however, are not 
known. Sato et at. [1989] and Kampfmann and Berckhemer [1985] downplay the 
importance of the melt phase on anelasticity, at least for temperatures not significantly 
above the solidus, and instead suggest that variations in temperature alone are sufficient to 
produce observed Q variations in the Earth. In their experiments with peridotite Sato et al. 
[1989] find that there is no large change in Q on first melting but rather a steady change 
with temperature. Kampfmann and Berckhemer [1985] similarly observe no large effect at 
slightly super-solidus temperatures. However, the Kampfmann and Berkhemer [1985] 
samples included significant melt, so it is difficult to separate clearly the effects of melt 
from those due solely to temperature. Goetze [ 1977] suggests that extrapolated laboratory 
creep data on unmelted olivines are compatible with geophysical evidence regarding the 
rheology of the upper mantle. Solomon [1973] reported a narrow (no wider than 100 km 
and shallower than 50 to 150 km deep) zone of low Q centered along the Mid-Atlantic 
Ridge axis, which he attributed to the presence of partial melt, and Molnar and Oliver 
[1969] found that Sn does not propagate efficiently in the immediate vicinity of the Mid- 
Atlantic Ridge, consistent with the presence of a narrow low-Q zone at the ridge. The 
observations in this work lack the resolution to discern the presence of such a narrow low- 
Q feature, but our data are consistent with the hypothesis that low Q beneath young 
oceanic lithosphere is due simply to elevated temperature. 
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The values that we obtain for Q* 1 averaged over depth in the lithosphere are high, 
with Q' 1 values ranging from 0.016 to 0.033 (Q= 30 to 60), with the range in Q' 1 values 
associated with the range of lithospheric ages from 0-100 My. This range is consistent 
with the values (Q* 1 = 0.018 to 0.052, Q =19 to 55) found for Q* 1 in the upper 220 km 
beneath the Lau back-arc spreading center [Flanagan and Wiens , 1990], but higher than 
the value of 0.014 for Q* 1 (Q = 70) at depths between 75 and 150 km beneath the East 
Pacific Rise region obtained by Ding and Grand [1987]. The Ding and Grand [ 1987] 
study did not report Q* 1 values for depths shallower than 75 km. Canos and Mitchell 
[1981] examined attenuation of Rayleigh waves in the north Atlantic and determined Q" 1 as 
a function of depth and age; their Figures 6 and 7 can be compared with our Figure 3.13. 
They found that a low-Q zone is prominent in the upper mantle of all regions younger than 
65 My in age but that such a zone is poorly developed in other regions. The Q* 1 values in 
their low-Q zone range from approximately 0.010 to 0.020 (Q = 50 to 100). 

We may also estimate the average difference in Q* 1 between S and SS in the upper 
mantle by means of the simple equation of 8t* to the product of differential travel time and 
average Q- 1 , where the travel time is the differential travel time of SS and S and t* has the 
value of 8t*. We obtain Q- 1 = 0.012 (Q = 82). This value likely represents an average 
Q-l across the oceanic upper mantle (0-650 km). This value is consistent with the average 
upper mantle Q* 1 value of 0.012 (Q = 82) obtained in a study of multiple ScS waves by 
Revenaugh and Jordan [1987]. Thus our 5t* measurements are consistent with other 
observations of seismic attenuation in the lithosphere and asthe nosphere in general and of 
low Q near active spreading centers in particular. 

As a further test of our Q -1 -temperature relation, we take the upper mantle 
temperature variations obtained from the inversion of SS-S differential travel time 
residuals, geoid, and bathymetry along the axis of the Mid-Atlantic Ridge [Sheehan and 
Solomon , 1991; Chapter 2] and predict the corresponding pattern of 8t* along the ridge 
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axis. In Figure 3. 14 we compare the observed pattern of along-axis variations in 8t* with 
those produced assuming the Q" ^-temperature relation from equa tio n 3.8. The predicted 
value of St* is obtained from equation 3.9, where T(z) is the sum of the Parsons and 
Sclater [ 1977] geotherm at a reference age and the average differential temperature AT 
obtained from the inversion. Below the lithosphere, T(z) is obtained by adding AT to an 
adiabatic temperature gradient We found that the predicted St* variations depend on the 
choice of reference geotherm. To obtain an average representative St* variation for the 
range of lithospheric ages that we sample, we calculate 5t* using the geotherm for each of 
the seven ages of the data groups in Figure 3.12, and we perform an arithmetic average to 
obtain the final St*. Along-axis profiles were calculated for all four Q- temperature 
relations examined in this study. 

Fits between the observed and predicted profiles of St* shown in Figure 3. 14 are not 
as good as we might have expected after a visual inspection of the apparent qualitative 
correlation between travel time and St*. We obtain variance reductions of about 25% for 
the models shown. A large part of the misfit in phase is due to an offset in the latitude at 
which the peak residuals in travel time and St* are observed in the vicinity of 30*N. Both 
the observed and predicted St* profiles indicate low attenuation in the region near 50*N. If 
the observed St* profile had been included in the inversion it is likely that a better fit could 
have been obtained for the St* data but at the expense of the fits to the geoid, travel time, 
and bathymetry profiles. The along-axis St* profiles obtained using the Sato and Sacks 
[1989] Q- temperature relation are somewhat different (Figure 3.15), with relative 
variations in St* larger than the observed. 

Profiles of St* were also constructed from temperature variations obtained from the 
joint inversion of travel time, geoid, and bathymetry for both temperature and composition 
[Sheehan and Solomon, 1991; Chapter 2]. These temperature profiles are influenced more 
strongly by the travel time residuals (relative to geoid and bathymetry) than are those 
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produced by the inversion for temperature only. No attempt has been made to include a 
compositional dependence of Q' 1 , so the model 8t* values are those predicted from the 
thermal perturbations only. We find that the 8t* profiles predicted from these temperature 
variations are generally in somewhat better agreement with the 8t* observations than those 
constructed using the temperature perturbations resulting from the inversion for 
temperature variations only. 


Conclusions 

We have measured SS-S differential attenuation 5t* in the north Atlantic region from 
about 150 wave pairs. After correcting 8t* for epicentral distance we find that 8t* 
decreases with increasing seafloor age. We do not observe evidence for a narrow region 
of very low Q along the axis of the Mid-Atlantic Ridge as reported by Solomon [1973], 
although this may be due to the limited spatial resolution of the long-period shear waves 
utilized in our study. 

We have derived empirical Q -1 -temperature relations by comparing measurements of 
8t* with values predicted under the assumptions that St* variations arise from lithospheric 
cooling and that the temperature structure as a function of lithospheric age is that given by 
the plate cooling model [Parsons and Sclater, 1977]. The Q" ^ - temperature relation that 
best fits our observations predicts smaller variations of Q^with temperature, especially at 
young ages, than the Q-tempcrature relation of Sato and Sacks [1989] derived from 
laboratory measurements of Qp' 1 in a spinel lherzolite at ultrasonic frequencies. 

Systematic variations of 8t* along the axis of the Mid-Atlantic Ridge are also seen. These 
variations are broadly consistent with the along-axis variations in temperature derived from 
an inversion of differential travel time residuals and geoid and bathymetry anomalies. 
Theoretical along-axis profiles of 8t* are constructed from the along-axis variations 
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derived in Chapter 2 of this thesis and the Q-temperature relation derived here. 
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APPENDIX 3 Jil ESTIMATION OF ERRORS FOR SS-S DIFFERENTIAL ATTENUATION (5t*) 

It is important to quantify the uncertainties in the measurements from individual and 
stacked spectra of 5t*. After windowing and calculating spectra, the “quality” of each 
individual St* measurement is rated and a grade is assigned. Our assignment of quality is 
largely subjective and based upon visual inspection of the waveforms and amplitude 
spectra of S and SS waves, and of samples of noise background, taking into account the 
clarity of the seismogram and shape of the amplitude spectral ratio. An "A" quality grade 
indicates that both S and SS pulses are well above the noise level and the amplitude 
spectral ratio shows a smooth decrease with frequency. "B" quality indicates good signal 
to noise ratios but an amplitude spectral ratio not as smooth, and a "C" quality grade 
indicates lower signal to noise ratios or log spectral ratios poorly fit by a straight line. In 
addition to A, B, and C grades, there were data that were rejected because of a poor signal 
to noise ratio for either the S or SS phase, or an irregular spectral ratio. 

As an objective means to obtain estimates of error, we examine the scatter in 
measurements of various quality within a small region. We measured the root mean 
squared (rms) difference between St* of the same grade (A, B, or C) with bounce points 
separated by less than 100 km and with differences in path azimuth at the bounce point of 
less than 10*. A 100-km distance is less than the horizontal wavelength of SS (which is 
about 180 km at 25 s period), so we do not expect much contribution to the rms difference 
from lateral variation in structure. The rms difference for 8 residual pairs of A quality 
which were within 100 km of each other was 3.4 s. For B quality picks, an rms 
difference of 3.8 s was measured using 17 residual pairs, and for C quality picks 20 
residual pairs yielded an rms difference of 3.3 s. The rms values are strongly affected by 
a few outliers, so rms values were also calculated without these estimates. The resulting 
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r ms values were 1.9 s for 6 A-quality St* pairs, 2.8 s for 13 B-quality pairs, and 2.6 s for 
17 C -quality pairs. Part of the reason that the rms values for the C-quality data are less 
than those for the B-quality measurements is that the mean values of the first group of 
measurements are smaller and thus the rms differences are also smaller. The C-quality 
measurements are often smaller because of increased noise at high frequencies which 
results in a flattening of the slope of the spectral amplitude ratios. 

Because of the bias of C-quality estimates towards low values of 8t*, we use the 
estimates of the rms differences only as an approximate guide for estimating the average 
overall errors in the A-, B-, and C- grade measurements. Our final choice of measurement 
errors used for relative weighting of the differently graded measurements is — 1.9 s, 
a B = 2.7 s, and cj c = 3.5 s. In the weighted regression experiments the A-, B-, and C- 
quality measurements are weighted inversely by their measurement variance. 
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Figure 3.1. 

Figure 3.2. 

Figure 3.3. 
Figure 3.4. 

Figure 3.5. 


Figure Captions 

An example of the measurement of SS-S differential attenuation for the 
event of December 24, 1985, at GDH. (a) Displacement seismogram with 
S and SS phases windowed, (b) Natural log of the amplitude spectrum of 
the S phase, (c) Natural log of the amplitude spectrum of the SS phase, 
(d) Natural log of the ratio of the amplitude spectra of the SS and S 
phases. Solid line is a least squares fit to the ratio from 0.01 to 0.08 Hz. 

Distribution of earthquakes (triangles) and seismograph stations (circles) 
used to measure SS-S differential attenuation. Stations are from the 
GDSN, NARS, and GEOSCOPE digital networks. Earthquakes are from 
the Harvard CMT catalogue (generally mb > 5.0) from the years 1977- 
1987. Lambert equal-area projection with the pole of projection at 45'N, 
40*W. 


Distribution of SS bounce points. Map projection as in Figure 3.2. 

SS-S differential attenuation (8t*) versus SS-S travel time residual. Each 
point represents the weighted mean of 38 data points adjacent in SS-S 
residual (x-axis). Weights are constructed from variances as discussed in 
Appendix 3.A. Horizontal and vertical bars are standard errors of the 
means. Linear regression yields a slope of 0. 10 ± 0.04. 

SS-S differential attenuation versus epicentral distance. Each point 
represents the weighted mean of 38 points adjacent in epicentral distance 
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(x-axis). Horizontal and vertical bars are standard errors of the means. 
Linear regression yields a slope of 0.05 i 0.02 s/ degree (line). 

Figure 3.6. Comparison of observed SS-S differential attenuation versus epicentral 
distance with predictions from several radial Q models [Anderson and 
Hart, 1978; Anderson and Given, 1982]. 


Figure 3.7. Distance-corrected SS-S differential attenuation (5t*) versus SS-S travel 
time residual. Each point represents the weighted mean of 38 adjacent 
data points. Linear regression yields a slope of 0.09 ± 0.04. 

Figure 3.8. (a) Distance-corrected SS-S differential attenuation (8t*) versus square 

root of seafloor age. Each point represents the weighted mean of 38 
adjacent da ta points. Horizontal and vertical bars are standard errors of 
the means. Linear regression yields a slope of -0.20 ± 0.07 s/ (My) 1 / 2 for 
an age range of 0-100 My (solid line). 

(b) SS-S travel time residuals versus square root of seafloor age in the 
north Atlantic, from Sheehan and Solomon [1991]. Each point represents 
the weighted mean of 14 adjacent data points. Horizontal and vertical bars 
are standard errors of the means. Linear regression yields a slope of -0.68 
± 0.08 s/ (My) 1 / 2 for a 0-100 My age range (solid line) or -0.76 ± 0.09 s / 
(My) 1/2 for a 0-80 My range (dashed line). 

Figure 3.9. (a) SS-S differential attenuation, corrected for epicentral distance and 

lithospheric age, versus latitude along the Mid- Atlantic Ridge, 10-70*N. 
Values shown are moving averages (such that each point is used twice) of 
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Figure 3.10. 


Figure 3.11. 


Figure 3.12. 


12 adjacent data points from lithosphere of age 0- 100 My. See the text for 
data reduction procedures. The location of the Iceland and Azores hotspots 
are indicated. 

(b) Age-corrected SS-S travel time residual versus latitude along the Mid- 
Adantic Ridge, 10-70‘N. The residuals shown are moving averages of 10 
adjacent data points from lithosphere of age 0-100 My. 

Geotherms at several different lithospheric ages predicted by the plate 
cooling model of Parsons and Sclater [1977]. Solidus for peridotite 
(dashed line) after Takahashi [1986]. 

SS-S differential travel time residual versus plate age, as predicted by the 
plate cooling model at a number of discrete ages. Linear regression yields 
a slope of -0.64 ± 0.01 s / My 1 / 2 for a 0-100-My age range (solid line). 
Note that a 50-My age corresponds to zero residual by convention. 

Observed (solid circles, with error bars) and predicted (triangles and 
squares) St* values versus age. Triangles show St* values for the 
lithospheric portion of the SS wave path calculated assuming plate cooling 
geotherms, the Takahashi [1986] peridotite solidus, and the Q- temperature 
relation of Sato and Sacks [1989]. Squares indicate St* values for entire 
upper mantle (0-500 km) portion of the SS wave path from nonlinear 
iterative inversion of the observed St* versus age data for the constants A 
and E in the Q- temperature relation Q" 1 = A’ 1 exp [-E / RT(z)]. This 
relation is assumed to hold over 0-125 km depth; a constant Q of 1 30 is 
assumed for 125-500 km depth. 
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Figure 3.13. (a) Qs versus depth (logio scale) at a lithospheric age of 20 My, 

calculated from nonlinear least squares inversion of observed 5t* versus 
age assuming different parameteri zations of Q in terms of temperature and 
pressure (dashed lines), and from forward modelling with the Sato and 
Sacks [ 1989] Q-temperature relation (solid line). Middle dashed line, Q* 1 
= 0.6 exp [-35/RT(z)]. Short dashed line, Q 1 = 0.35 exp [-2.05 
T m (z)/T(z)]. Long dashed line, Q* 1 = [2.0 + P(z)]* 1 exp [-1.7 
T m (z)/T(z)]. 

(b) Same as (a) but at 73 My age. 

Figure 3.14. Observed (solid line) and predicted (dashed lines) 8t* along the Mid- 

Atlantic Ridge, 10-70'N. “Observed” profile is actually a filtered version 
of the observations, containing only the wavelengths 1400 to 7100 km for 
comparison with the along-axis pattern of temperature variations inferred 
from a joint inversion of travel-time residuals, geoid heights, and residual 
bathymetry [Sheehan and Solomon, 1991; Chapter 2]. The predicted 5t* 
profiles are calculated assuming one of the Q~ ^ - temperature relations 
derived in this study (Q‘* = A'^ exp [-E/RT]) and the upper mantle 
temperature variations from models of Sheehan and Solomon [1991]; see 
also Chapter 2. Two of the predicted profiles were calculated using the 
temperature profiles [Sheehan and Solomon, 1991; Chapter 2] for models 
with temperature variations constrained first to be in the upper 150 km, 
then in the upper 300 km; the viscosity structure includes a high viscosity 
lid over a constant-viscosity mantle. The other predicted profiles were 
calculated from the temperature variations of Sheehan and Solomon 
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[1991] in which a joint inversion is performed for variations in both 
temperature and composition in the upper 150 km, and then upper 300 km 
for the same viscosity model. Mean values have been subtracted from 
both the observed and model profiles. 

Figure 3.15. Observed (solid line) and predicted (dashed lines) 5t* along the Mid- 
Atlantic Ridge, 10-70*N. The predicted 8t* profiles are calculated 
assuming the Sato and Sacks [1989] Q- temperature relation (equation 
3. 10) and the upper mantle temperature variations from models of 
Sheehan and Solomon [1991]; see also Chapter 2. Two of the predicted 
profiles were calculated using the temperature profiles [Sheehan and 
Solomon, 1991; Chapter 2] for models with temperature variations 
constrained to be in the upper 150 km, then in the upper 300 km; the 
viscosity structure includes a high viscosity lid over a constant-viscosity 
mantle. The other predicted profiles were calculated from the temperature 
variations of Sheehan and Solomon [1991] in which a joint inversion is 
performed for variations in both temperature and composition in the upper 
150 km, and then upper 300 km for the same viscosity model. 
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Figure 3.1a 


124 



Figure 3.1bc 
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Figure 3.2 
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CHAPTER 4 

Upper Mantle structure in the vicinity of the east 
Pacific Rise inferred from Shear Wave differential 

TRAVEL TIMES, GEOID, AND BATHYMETRY 

INTRODUCTION 

The formation and evolution of the lithosphere and the interaction of mantle 
convection with the surface plates can, in principle, be constrained by the determination of 
variations in shear wave velocity in the upper mantle. Oceanic regions, and in particular, 
mid-ocean ridge systems, are most amenable to examination of the processes which drive 
the plates due to their simple crustal structure and young age relative to the continents. In 
this chapter we investigate lateral variations in the velocity structure beneath the East 
Pacific Rise by measuring differential travel times for the phase pair SS-S, and relate these 
variations in travel time to depth and geoid anomalies in the region. Body wave methods 
are well suited for resolving lateral variations in upper mantle structure at relatively short 
(<1000 km) wavelengths. Differential travel times of shear wave phase pairs, in 
particular, are an effective tool in the study of upper mantle heterogeneity [Sipkin and 
Jordan, 1976, 1980a; Stark and Forsyth, 1983; Butler, 1979; Kuoet al., 1987; 
Woodward and Masters, 1991; Sheehan and Solomon, 1991] and have the advantage that 
source and receiver effects are approximately common to both phases and are thus largely 
eliminated by differencing. In addition, use of differential travel times eliminates any 
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absolute tuning errors. Differential times also have the ability to constrain structural 
variations in regions of the world where there are no sources or receivers. Under the 
assumptions that the lower mantle is relatively homogeneous and that the paths in the 
upper mantle are reasonably steep, the differential travel time anomaly can be associated 
with upper mantle structure in the vicinity of the surface bounce point of the reflected (SS) 
phase. 

A number of previous studies have brought important contributions to our 
understanding of upper mantle properties of the eastern Pacific. Surface wave dispersion 
studies in the Pacific Ocean region [e.g., Leeds , 1975, Forsyth, 1975, Yu and Mitchell, 
1979; Nishimura and Forsyth, 1985, 1988, 1989; Zhang andTanimoto, 1990b] show a 
clear relationship between the age of the seafloor and the seismic properties of the mantle. 
Woodward and Masters [1991], in a global study of SS-S and PP-P differential travel 
times, found the east Pacific region to be particularly slow. In the present study, we use 
shear wave differential travel time residuals, in combination with geoid and bathymetry 
data, to elucidate further the structure of the east Pacific region. Our data set allows us to 
search for velocity anomalies associated with the East Pacific Rise and with the Galapagos 
Spreading Center and Galapagos hotspot. In addition, we can compare the results from 
this intermediate to fast spreading ridge system with our earlier study of the northern Mid- 
Atlantic Ridge, a slow spreading ridge system [Sheehan and Solomon, 1991]. Several 
authors [e.g., Stark and Forsyth, 1983; Zhang and Tanimoto, 1990b, Woodward and 
Masters, 1991] have found that shear wave velocities in the east Pacific are significantly 
slower than those observed in either the Atlantic or Indian Oceans. In addition, the 
presence of anisotropy has been found to be more pronounced in the Pacific upper mantle 
than in the Atlantic upper mantle [Montag ner and Tanimoto, 1990] as is predicted by 
models for flow-induced (mentation of mantle olivine crystals [McKenzie, 1979; Ribe, 
1989]. Woodward and Masters [1991] found that their Pacific data were consistent with 
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the presence of weak anisotropy, although they suggested that lateral heterogeneity could 
also have produced the observed azimuthal pattern. They did not resolve significant 
azimuthal variations in the Atlantic portion of their data set 

Our analysis technique, developed in Chapter 2, involves the use of shear wave 
differential travel time residuals in combination with geoid and bathymetry data . The 
combination of geoid (or gravity) data with seismic data is an effective way to determine 
Earth structure on scales ranging from global [Hager and Richards , 1989; McNutt and 
Judge, 1990] to local [Lines et al., 1988; Lees and VanDecar, 1991]. Several authors 
[Watts et al., 1985; Haxby and Weissel, 1986] have suggested that bathymetry and 
especially gravity (or geoid) data in the central and east Pacific regions show evidence for 
small-scale (wavelength of a few hundred kilometers) convection. With our combined 
data set and analysis we can search for evidence of small scale convection with the added 
constraint of shear wave travel times. 

GEOLOGIC SETTING 

The East Pacific Rise lies to the west of Central and South America and stretches 
over 6000 km from the Juan Fernandez triple junction at 34*S to the Gulf of California. 
Segments of the East Pacific Rise spread at rates as high as 162 mm/yr (full spreading 
rate) at about 20*S (Figure 4.1). In comparison, spreading rates on the Mid-Atlantic 
Ridge range from just 5 to 40 mm/yr. The East Pacific Rise separates the large Pacific 
plate on the west from the Nazca and Cocos plates on the east. The triple junction between 
these three plates lies at 2'N, 102*W. The Cocos plate is bounded on the west by the East 
Pacific Rise, on the south by the Galapagos Spreading Center, and on the northeast by the 
Middle America Trench. The Cocos plate includes the aseismic Cocos Ridge, which 
trends NNE from 1*N, 90*W to 7*N, 84*W. A related feature is the Carnegie Ridge on 
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the Nazca Plate from 0*N, 90’W to O’N, 82*W. A rough-smooth bathymetric boundary is 
located at a distance of about 250 km to the west of the Cocos and Carnegie Ridges. The 
Cocos plate is shallower and its gravity is higher than that predicted by plate cooling 
models [Parsons and Sclater, 1977; Parsons and Richter, 1980]. Spreading rates on the 
north-south spreading Galapagos Spreading Center average around 70 mm/yr 
[Macdonald, 1989]. The Galapagos hotspot, just south of the Galapagos Spreading 
Center at about 100*W longitude, is a site of active volcanism. Monnereau and Cazenave 
[1988] found that the geoid-to- topography ratio in the region of this hotspot is not 
consistent with that found for swells in the center of a plate. They attributed this 
difference to the proximity of the Galapagos hotspot to the Galapagos Spreading Center. 

DATA 


Travel Time Data 

The east Pacific was chosen as the study area because the very high spreading rates 
on the East Pacific Rise offer an interesting complement to our study of the slow- 
spreading Mid-Atlantic Ridge (Chapter 2) and because there is a convenient distribution of 
earthquake sources and digital seismic stations to allow for a sufficiently dense sampling 
of SS bounce points. We performed a search over all earthquakes in the Harvard centroid 
moment tensor (CMT) catalog for the years 1980-1987 [Dziewonski et al., 1981; 
Dziewonski and Woodhouse, 1983] supplemented by a listing of intraplate earthquakes 
not included in the CMT catalog (E. A. Bergman, personal communication, 1991) and 
over all GDSN digital seismic stations in order to find event-station pairs of the proper 
epicentral distance which provide SS bounce points in the East Pacific Rise region. 
Searches were also performed to find event- station pairs recoverable using the 
GEOSCOPE broadband stations in the time periods for which we have data (January to 
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December of 1986 and January to May of 1988). The range in source-receiver separation 
was taken to be 55* to 86* to ensure separation of S and ScS at the longer distances and to 
avoid triplication in SS at shorter distances. A list of stations used in this study is 
presented in Table 4.1. We use only transversely polarized (SH) seismograms (rotated 
from N-S and E-W components) to avoid interference from the SKS phase and 
contamination from P-SV conversions at the base of the crust and other near-surface 
discontinuities. Epicenters were obtained from the "Preliminary Determination of 
Epicenters" of the U.S. National Earthquake Information Service (NEIS) for all events. 
The final distribution of sources and stations used to measure SS-S differential travel times 
is shown in Figure 4.2. Our study area extends from 25’S to 30*N latitude and 140*W to 
80*W longitude, and includes a dense sampling of the Cocos plate. The majority of data 
in this study comes from records of southern East Pacific Rise, Chile Rise, and Peru-Chile 
subduction zone earthquakes recorded at North American stations. Other data include 
Peru-Chile and Middle America trench events recorded at Hawaiian stations, Hawaiian 
events and events from western North America recorded at South American stations 
(BOCO and ZOBO), Central and North American earthquakes recorded at the DWWSSN 
station AFI (Afiamalu, Western Samoa), intraplate events on the Nazca plate propagating 
to North America and Hawaii, and central Pacific intraplate events propagating to North 
American stations. 

This search yielded over 1500 event-station pairs with the proper epicentral 
separation. Data with SS bounce points on regions of thick sediments (> 0.5 km) and 
thick crust off the coast of Central America [Ludwig and Houtz, 1979; Winterer , 1989] 
were omitted. After winnowing the list because of station inoperation, poor signal to 
noise ratio for the phases of interest, and interfering events, the final data set consists of 
over 600 SS-S differential travel time residuals with bounce points in the east Pacific 
(Figure 4.3). A total of 21 digital stations (Table 4. 1) and 342 different earthquakes 
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(Table 4.2) were used in the final analysis. Uncertainties for each measurement are 
adopted from the analysis in Chapter 2, Appendix 2.A. A tabulation of all residuals is 
given in Appendix HI. 

Bathymetry Data 

Bathymetric data are obtained from the corrected Digital Bathymetric Data Base (S’ 
grid) (DBDB5C) [US. Naval Oceanographic Office , 1985], DBDB5C data are corrected 
for the deviation of water column acoustic velocity from die generally assumed value of 
1500 m s' 1 . We obtain an average depth value for each SS surface bounce point by 
averaging all of the bathymetric data within a 1* x 1* square centered on each point. 

The bathymetry of the East Pacific Rise region is shown in Figure 4.3. The rise axis 
is marked by depths generally shallower than 3000 m. To the west of the rise the seafloor 
depth gradually increases with increasing age of the seafloor. At the western edge of the 
map the depth is about 5000 m. To the east of the East Pacific Rise, however, the 
bathymetry does not closely follow simple depth-age relations and is complicated by the 
presence of the Galapagos hotspot and by the Middle America subduction zone. The 
depth along the crest of the Galapagos Spreading Center, about 2500 m, is generally 
shallower than that of the East Pacific Rise. 

Geoid Data 

Geoid data are taken from a combined set of Seasat and GEOS3 altimeter data 
[Marsh et al., 1986]. The sea surface elevations are referenced to the IUGG 1980 
Geodetic Reference System defined by an Earth semimajor axis of 6,378,137 m and a 
flattening coefficient of 1/298.257 [Moritz, 1980]. The mean sea surface is presented in 
the form of a 0.25* grid. The data have been corrected for orbit errors, instrument and 
atmospheric propagation effects, and solid Earth and ocean tides. Since we are only 
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interested in signatures related to upper mantle processes, we attempt to remove from the 
data all other effects. Therefore, we retain only the components of the field with 
wavelengths less than about 6000 km by subtracting the observed field up to degree and 
older 7 assuming the GEM-T2 coefficients [Marsh et al., 1990] for the spherical harmonic 
expansion of the Earth’s geoid and tapering the spectrum to degree and order 1 1 to avoid 
truncation effects [Sandwell and Renkin, 1988] (Figure 4.4). The geoid field over the 
East Pacific Rise and Cocos plate is characterized by positive anomalies over the East 
Pacific Rise that fall off gradually with increasing plate age to the west To the east of the 
rise, however, the geoid does not follow the simple relation predicted by the plate cooling 
model [Parsons and Richter , 1980] but rather is characterized by positive values over 
much of the Cocos plate and a large positive anomaly presumably related to the slab 
subducting under Central America [Hager, 1984]. The Galapagos hotspot is marked by a 
short (400 km) wavelength geoid high superimposed on a northwest-southeast-trending 
regional field with a low of -4 m. 


METHOD 


Data processing procedures and the waveform cross-correlation technique for 
estimating differential travel time are as in Chapter 2. 


RESULTS 


Spatial Pattern of SS-S Residuals 

The SS-S residuals are shown plotted at the SS bounce point in Figure 4.5. Much 
scatter is observed, and a clear signal associated with the rise is not obvious. The 
residuals far to the west of the rise (from longitude 130* to 140*W) are more negative 
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(earlier SS) than the residuals in close vicinity to the rise, but a gradual decrease from slow 
at the rise to fast off the rise is not apparent. Residuals to the south of the Galapagos 
Spreading Center tend to be less by 1 or 2 s than those from bounce points to the north of 
the Galapagos Spreading Center. The values of the SS-S residuals with bounce points on 
the Cocos plate vary widely, though some coherent patterns emerge. There are groups of 
early and late residuals with a characteristic wavelength of approximately 1200 km, 
patterns which are made more obvious when the data are smoothed and gridded. The 
residuals in the vicinity of the Galapagos hotspot are positive. Positive residuals are also 
found about 100-200 km to the northeast of the Galapagos hotspot 

Lithospheric Aging 

Cooling and thickening of the lithosphere should yield a tendency toward an increase 
in seismic velocity with increasing lithospheric age. A linear regression experiment was 
performed to examine the correlation of the SS-S residuals with seafloor age. First a 
gridded map of seafloor ages was constructed for the East Pacific Rise from a compilation 
of data from several sources. Gridded age data from the map of Larson et al. [1985] were 
provided by S. C. Cande (personal communication, 1991). Since these data only included 
the age range of 0 - 6 My in our study area, we had to supplement them with data from 
other sources. These additional sources included digitized isochrons from Sclater et al. 
[1981] and magnetic anomalies from Klitgord and Mammerickx [1982], Atwater [1989], 
and Atwater and Severinghaus [ 1989] and from the tectonic map of the Circum- Pacific 
Map Project [1981]. The magnetic isochrons are assigned ages according to the polarity 
reversal time scales of Kent and Gradstein [1986] and Berggren etal. [1985]. Plate age is 
not well determined on parts of the Cocos plate due to its proximity to the magnetic 
equator [Schouten, 1971], and thus no ages have been assigned to such areas [Atwater 
and Severinghaus, 1989]. We place constraints on seafloor age in such regions lacking 
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identified isochrons by assuming symmetric spreading and using the ages from isochrons 

at an equal distance on the opposite side of the rise. 

To obtain a reprsentative age value for the region spanning approximately one 
horizontal wavelength of the incident (SS) wave, an average seafloor age was estimated 
for a 1 # x 1* box centered on each SS bounce point To reduce scatter, measurements 
whose bounce point depths differed by more than 1200 m from the depth predicted by the 
regional depth-age relation were excluded from the age regression. 

The SS-S residuals for the East Pacific Rise are shown plotted against the square 
root of seafloor age in Figure 4.6a. The solid line represents the Unear travel time - age 1 / 2 
relation derived from the data, but the figure indicates that the data are not well fit by a 
straight line. The age dependence derived in Chapter 2 for the plate cooling model is 
shown as a dashed line. The predicted slope is consistent with the longest wavelength 
trends in the data and suggests that perhaps there is an age dependence to the travel times 
but that several other competing processes are also contributing. The observed SS-S 
residuals increase from 0 to 2.5 MyW partly because the Galapagos hotspot region is 
slow but not on zero-age lithosphere. Large positive residuals at 4.5 My 1/2 include SS 
bounce points off the coast of Central America, perhaps due to unmodeled sediments or 
increased crustal thickness. 

We have c alculated the “effective age” of the Uthosphere by taking the depth and 
estimating the age of seafloor of that depth predicted by the plate cooUng model. This 
procedure is consistent with the “Uthospheric rejuvenation” hypothesis of Menard and 
McNutt [ 1982]. The term “rejuvenation” refers to the observation that many midplate 
swells cool and subside at the same rate as normal (younger) Uthosphere at the same 
depth. Thus, thermal rejuvenation appears to produce lithosphere with many of the 
properties of normal Uthosphere of a younger age. The Galapagos hotspot and the 
shallow Cocos and Carnegie ridges, for example, are regions of anomalously shallow 
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depth in our study area, and are presumably related to some sort of hearing event For 
other anomalously shallow regions included in this study (such as the majority of the 
Cocos plate), the mechanism for shoaling is less clear. The slope of the SS-S travel time 
residual versus square root of predicted age is -0.26 ± 0.04 s/My 1 ^ for the age range 0-33 
My (Figure 4.6b). This value is at least of the same sign as the relation predicted by plate 
cooling and with that observed in the Atlantic (Chapter 2), but the magnitude of the slope 
is still significantly lower than expected from the model. 

To look for other systematic variations in the residuals, we correct for age by 
removing the linear dependence on square root of age predicted by the plate cooling 
model, approximately -0.7 s My 1 ^ 2 (Chapter 2). This value for slope is also in good 
agreement with observations in the north Atlantic but admittedly is steeper than what we 
have observed with our limited range of east Pacific data. Future work is warranted to 
study further the relation between travel time and plate age in the east Pacific. 

The travel time residuals in the region of the East Pacific Rise are about 2 s larger 
at a given bounce point age than those in the Atlantic. Part of this difference may be due to 
the fact that since the Mid-Atlantic Ridge is slow spreading the zone of lithosphere of 
young age is extremely narrow, whereas in the fast spreading Pacific it extends over a 
much greater width. So even if the age at two SS surface reflection points is the same in 
both of these studies, the upper mantle portion of the path in the Atlantic region will 
include a greater proportion of older lithosphere than the Pacific region. The magnitude of 
this effect would be about 1 s at young lithospheric ages (0-15 My) and less at older ages. 

Anisotropy 

Another systematic velocity variation that has been suggested as a possible 
contributor to residual SS-S travel times is azimuthal anisotropy. As noted in Chapter 2, 
the SS-S data from the north Atlantic do not show a consistent anisotropic effect We 
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performed least squares inversions on the age-corrected SS-S residuals from the Pacific to 
determine best-fitting 48 patterns. The anisotropy indicated by our regression experiments 
differs significantly from those found in the Atlantic. Our results indicate that for the 48 
fit to the SS-S residuals uncorrected for age the slow directions few SS-S are N44 W and 
N46*E and the peak-to-peak magnitude of the effect is 2.7 s (Figure 4.7a); for the 48 
model fit to the age-corrected SS-S data the slow directions are N48*W and N42*E and the 
magnitude is 3.4 s (Figure 4.7b). Our inversion for a 48 pattern of anisotropy provided a 
variance reduction of 36% for both the uncorrected and the age-corrected SS-S residuals, a 
reduction that is significant considering the number of free parameters involved. 

These results should be interpreted with caution, however, as our modelled 
azimuthal pattern may be partly the result of unmodelled upper mantle heterogeneity, and 
we lack data at many azimuths. The 48 model fits the data well in the azimuth range -40* 
to +20*, but the fit is not good at -45* (where we have a peak in the predicted 48 pattern), 
and the correlation of model to data is also poor in the vicinity of the +45* peak. The 
azim uthal distribution of the data is shown in Figure 4.8. Inspection of Figure 4.8 
indicates that the distribution of azimuths is not uniform, and that data of a given azimuth 
tend to be clumped together geographically. Thus it is difficult to determine whether the 
observed azimuthal signal is due to anisotropy or due simply to the geographic distribution 
of lateral heterogeneity. The majority of the data at east-west azimuths come from Central 
American events propagating to Hawaii, with bounce points in the region 20 N, -120 E. 
The data at north-south azim uths come predominantly from Chile Rise events propagating 
to North America. Most of these residuals are negative (early) and contribute to the 
negative trough in the azimuthal patterns (Figures 4.7a, b). The data most inconsistent 
with the 48 model come from the azimuth N60*W. These data are mainly from bounce 
points west of the rise on seafloor of approximately 50 Myr age. 



152 


Spatial Pattern of Age-corrected Residuals 

A plot of age-corrected SS-S travel time residuals at the SS surface reflection point 
(Figure 4.9) shows several interesting features. The pattern of alternating fast and slow 
residuals, with a characteristic wavelength of 1200 km, is still visible in the Cocos plate 
region. The Galapagos hotspot is marked by positive residuals. The residuals south of 
the Galapagos Spreading Center are for the most part more negative than those to the 
north. Residuals in the region 15*-30*N, 1 15*-130*W are positive. Far from the rise, at 
130*- 140* W longitude, the residuals show much scatter. Residuals along the coast of 
Central America are positive. 


Analysis of profiles 

The sparse sampling of our data does not permit us to contruct a profile along the 
axis of the rise as was performed for the north Atlantic in Chapter 2. However, we have 
been able to construct several profiles across the East Pacific Rise and Galapagos 
Spreading Center which allow us to examine more closely the relation of the observations 
to plate age, and, after removal of the age dependence, to analyze anomalies orthogonal to 
the ridge. The locations of the profiles are shown in Figures 4.3 and 4.4. 

Profiles of SS-S travel time, geoid, and bathymetry values were constructed by 
searching for all SS surface bounce points lying within 250 km (or 200 km for Profile 2) 
of the profile line and then applying a weighted averaging scheme to produce a single 
value every 50 km. To ensure complementarity of data sets, bathymetry and geoid height 
values are obtained at each SS bounce point, and all are corrected for cooling and 
subsidence with seafloor age. In this manner we effectively normalize all observations to 
zero age. The age correction for SS-S travel time data ( -0.7 s My 1 ^) was given in an 
earlier section. The age corrections for depth and geoid are described below. 
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The age correction for depth is based on an empirical depth-age relationship from 
Schroeder [1984], which accounts for the thermal subsidence of oceanic lithosphere away 
from spreading centers, and is appropriate for this part of the Pacific. The Schroeder 
[1984] age-depth relation is d(r) = 2967 - 305 r 1/2 , where d is depth in meters and t is time 
in My. We also examined depth anomalies obtained using the Parsons and Sclater [1977] 
depth-age relation, which was determined using depth data from outside of this area. The 
general trends observed using the Parsons and Sclater [1977] relation are not markedly 
different from those obtained using the Schroeder [1984] relation, but we prefer the 
Schroeder [1984] relation as more appropriate for this region. 

The geoid-age correction is done in two different ways. In the first method the 
theoretical geoid-age relation from the plate cooling model [Parsons and Richter, 1980] is 
removed to obtain the geoid anomaly. The second method, after Hager [1983], takes into 
account the fact that we have already removed the low order signal from the geoid. Hager 
[1983] made use of a global regionalization of plate ages [Mauk, 1977] and calculated the 
predicted effect of plate cooling on the geoid. The predicted field was then expanded in 
spherical harmonics up to degree and order 20. Our geoid-age correction after Hager 
[1983] then consists of the following: we remove the full predicted field up to degree and 
order 20, and then add the low order field (up to degree and order 7) back in. This 
prevents us from removing the low order part of the geoid-age dependence twice 
(empirically by the low degree and order reference geopotential field, and explicitly with 
the theoretical geoid-age relationship). This is especially critical in the very fast-spreading 
region of the Pacific included in this study, as the geoid-age dependence is dominated by 
the low order signal. This sort of correction is not critical in the slow-spreading Atlantic, 
where the geoid-age correction is dominated by wavelengths less than 4000 km, shorter 
than those removed to obtain our residual geoid [Cazenave et al ., 1986]. 

Error introduced into depth and geoid anomalies by isochron mislocation is difficult 
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to estimate precisely, but for an error in age of 2 My, depth and geoid errors at 80 My 
would be about 30 m and 0.2 m, respectively, while at 2 My, an error in age of 2 My 
would have a much larger effect, giving depth and geoid errors of 350 m and 0.3 m. The 
magnitude of this error highlights the importance of accurate age determination, especially 
at young ages. The age correction is relatively more important for bathymetric data than 
for die geoid, because the largest broad-scale variations in seafloor depth results from 
lithospheric cooling effects, unlike the geoid. 

Profile 1, East Pacific Rise south of triple junction 

Profile 1 starts to the west of the East Pacific Rise on seafloor approximately 45 My 
in age, passes over the ridge just south of the equator at 0.5 *S, 103.5 *W, and continues 
to the east of the rise to approximately 92 *W (Figure 4.3, 4.4). Figure 4. 10a shows the 
seafloor age along the profile, with the rise axis marked by a line. The geoid high (Figure 
4. 10b) is centered to the west of the rise and is asymmetric, characterized by a steeper 
gradient over the east flank of the rise than to the west The bathymetry (Figure 4. 10c) 
correlates well with plate age, with the bathymetric high centered at the rise and a gentle 
slope to the wesL The SS-S travel time residuals along the profile are shown in Figure 
4. lOd. The data along this profile are all at similar azimuths (approximately N-S), so 
azimuthal anisotropy can be discounted as a source of the variations. The travel time 
residuals of Profile 1 do not show a strong correlation with age. In fact, the rise is marked 
by a local minimum in the travel time residuals, with larger residuals (later arrivals) 
immediately to the east and the west of the rise. This pattern is similar to the one found by 
Schlue [1981] in a study of SS-S differential attenuation in the same region. To the west 
of the rise the residuals decrease gradually, with short wavelength signals superimposed. 
To the east of the rise the residuals drop off rapidly with age, in a manner similar to the 
geoid signal, but then increase again at about 700 km east of the rise axis. 
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Figure 4. lOe shows the residual geoid signal after correcting for plate cooling with 
the method of Hager [1983]. In comparison, the geoid corrected according to the method 
of Parsons and Richter [1980] is shown in Figure 4.10f, and the higher order geoid 
(reference field up to degree and order 10 removed) is shown in Figure 4. lOi. The high 
older geoid does not have an age correction applied, as we assume that the age dependence 
is dominated by the low order signal [Cazenave et al ., 1986]. This predicted dominance 
of the low-order geoid-age signal leads to a very small geoid-age correction if only the 
high order portions of the predicted field are removed (Figure 4.10e). The age-correction 
in Figure 4. lOf reduces the geoid bulge in the immediate vicinity of the rise, but since the 
geoid data far to the west of the rise show little age dependence, the age-correction 
increases the geoid signal and produces a positive anomaly of about 4 m. The geoid to the 
east of the rise still falls off with distance, as it did before the age-correction was applied, 
but the slope is more gradual in the residual profile. The age-corrected bathymetry (Figure 
4.10g) is marked by a large negative residual with a width of about 700 km at the nse 
axis, indicating that the rise is 400 m deeper than predicted by the plate cooling model. 
Positive depth anomalies of approximately 200 m are found on the east and west flanks of 
the rise axis. This steep residual at the rise axis is produced using either the Parsons and 
Sclater [1977] or the Schroeder [1984] depth-age relation and was also found by Menard 
and Dorman [1977]. The age-corrected travel time residual profile slopes monotonic ally 
from west to east, with a full range of 6 s. A short-wavelength residual low at the rise is 
also present in this profile and in the high order geoid profile (Figure 4. lOi) but is not 
readily apparent in the lower order geoid profiles in Figures 4.10e and 4.10f. 

Profile 2, across the Galapagos Spreading Center 

Profile 2, approximately orthogonal to the Galapagos Spreading Center, starts to the 
south of the spreading center at 5*S, passes to the east of the Galapagos hotspot at 1 S, 
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and crosses the Galapagos Spreading Center at l’N, 89* W. The profile continues north of 
the spreading center onto the Cocos plate, stopping near the Middle America Trench at 
8.7'N, 87.3*W. A comparison of the bathymetry for profiles 1 and 2 shows their very 
different character. In profile 1 , the area west of the rise increased gently up to the rise, 
whereas in profile 2 bathymetry (Figure 4. 1 lc) is much more variable and seems to be 
strongly affected by the Galapagos hotspot and the Middle America Trench and less related 
to simple plate cooling. Because Profile 2 crosses the most tectonically disturbed area, it 
is expected to present the most difficulty for application of simple plate cooling models. 
Once again, data on this profile are all of similar azimuth (approximately N20*W) so the 
effects of azimuthal anisotropy on the shear wave travel time residuals can be neglected. 

We observe a positive correlation between both the age-corrected residual depth (4.1 lg) 
and travel time residual (Figure 4. 1 lh) profiles. The geoid signal is dominated by a 
monotonic increase from south to north of about 12 m (Figures 4. 1 lb,e,f), a pattern 
which is presumably at least partially related to the Middle America subduction zone. 

Hager [1984] has suggested that the signal of subducting slabs is strong up to degree and 
order 9. It is therefore worth examining the high -order geoid (Figure 4. 1 li) as we are 
mainly interested in processes unrelated to subduction. The high-order geoid shows some 
correlation with the bathymetry and travel time, mainly at long wavelengths and in the 
vicinity of the Galapagos hotspot. The Galapagos Spreading Center appears in the middle 
of a south- to- north, negatively sloping feature on the bathymetry, travel time, and high- 
order geoid. 

Profile 3, East Pacific Rise north of triple junction 

Profile 3 stretches from 10*N, 131*W to 9*N, 86*W, crossing over the East Pacific 
Rise at 9.5*N, 103.5*W. To the west of the rise the depth (Figure 4.12c) increases 
gradually with age. The geoid (Figure 4. 12b) rolls off steeply on the west side of the rise 
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and increases slightly at a distance of about 1200 km from the rise axis. The SS-S travel 
time residuals are low (comparatively early) at the rise axis, but in the region from 300 km 
to 2500 fon to the west of the rise, the residuals follow the predicted relation between age 
and travel time residual. The data to the east of the rise are all located on the Cocos plate, 
and, after initially decreasing with age to the east of the rise, increase significantly as the 
coast of central America is approached. The depth anomaly (Figure 4.12g) west of the 
East Pacific Rise is near zero and shows little deviation from the depth predicted by 
Schroeder [ 1984]. The main feature in the age-corrected depth and travel time anomaly 
profiles [Figures 4. 12g and 4. 12h] is the region of high values at the eastern end of the 
profile. The geoid profiles [Figures 4.12e, f, and i] are not weU correlated with the travel 

time and bathymetry. 

Prefile 4, East Pacific Rise north of triple junction 

Profile 4 stretches WNW - ESE from 23*N, 127*W to 7*N, 83*W, crossing over the 
East Pacific Rise at approximately 13.5*N, 104*W. As with the other E-W profiles 
(profiles 1 and 3) the bathymetry (Figure 4.13c) and travel time residuals (Figure 4.13d) 
follow the predicted age dependence to the west of the rise, whereas the residual geoid 
(Figure 4.13b) falls off more steeply than predicted to the west and then increases about 
1300 km west of the rise axis. For this profile the travel time residuals are relatively large 
(late) at the rise axis, as expected from plate cooling models, but contrary to observations 
in the previous profiles. The travel time, geoid, and bathymetry all fall off steeply to the 
east of the rise and then increase as the Central American coastline is approached. The 
age-corrected residuals of Profile 4 do not differ markedly from those of Profile 3. The 
fact that the travel time residuals (Figures 4. 1 3d, h) increase at the eastern end of the 
profile and then decrease sharply as the shallowest depths on the bathymetry profile are 
reached might be partly due to an inaccurate bathymetric correction to the travel times. 



158 


Inversion Experiments 

Inversion experiments were conducted to assess whether the geoid, topography, and 
travel time anomalies presented in Profile 2 can be explained by reasonable variations in 
temperature and/or composition in the upper mantle. Profile 2 was singled out for further 
analysis in order to search for thermal and/or compositional variations related to the 
Galapagos hotspot, and also because of the qualitative correlation between bathymetry and 
travel time at a variety of wavelengths (700 - 1500 km). The geoid, admittedly, is 
problematic in this region due to the signal from the Middle American subduction zone. A 
total of 24 inversions were run, with various combinations of inversion parameters 
(temperature and/or composition), the thickness of the layer in which temperatures or 
composition were allowed to vary (0-150 km, 0-300 km, or 0-650 km), viscosity 
structures (constant viscosity mantle or mantle with a shallow low-viscosity zone), and 
geoid corrections. Results of the inversion experiments are summarized in Table 4.3 and 
Figures 4.14-4.21. 

We find that compositional variations alone (parameterized in terms of Mg#) are 
inadequate to model simultaneously the geoid, travel time, and bathymetry observations. 
Temperature variations alone provide reasonably good fits to the data, with variance 
reductions ranging from 66% for the model with a low viscosity zone and temperatures 
constrained to vary only over the depth range 0-300 km, to 72% for the model with a low 
viscosity zone and temperature perturbations confined to the depth range 0-150 km. For 
the latter model, the range of predicted temperature variations is high (± 150*C), so the 
models with temperatures constrained to extend deeper are more reasonable, as they 
predict smaller temperature variations and thus lesser differences in melt production and 
crustal thickness. The most interesting feature in the Profile 2 inversions is the Galapagos 
hotspot region, marked by slow SS-S travel times, shallow topography, and a slight geoid 
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high- The ocher dominant feature in the profile is a long-wavelength southward slope, 
especially strong in the geoid, which may be partially due to the Middle American 
subdue tion zone. In an attempt to remove the slab signature, we also carried out 
temperature inversions using the high-order geoid (reference field up to degree and order 
10 removed). The high-order geoid is of much smaller amplitude than the low-order 
geoid, and we had to adjust our data covariance matrix to ensure proper weighting 
between the various data. The fits with the high-order geoid were worse than those with 
the low -order geoid, as only a small temperature variation in the upper mantle was 
required to fit the geoid but this variation was insufficient to match the observed travel time 
residual and bathymetry. It may be worthwhile to remove the slab signal from the geoid 
explicitly. 


DISCUSSION 


We found that the SS-S differential travel time is about 2 s larger in the Pacific at a 
given seafloor age than in the north Atlantic. The Pacific has been found to be seismically 
slower at a given age than the Atlantic (or Indian Ocean) in a number of studies [Stark and 
Forsyth, 1983; Zhang and Tanimoto, 1990b; Woodward and Masters, 1991]. In addition, 
heat flow in the Pacific is higher than in the Atlantic [Chapman and Pollack, 1975], 
suggesting that there may be some differences in the thermal structure beneath these two 
ridges or that there are differences in the fraction of heat flow earned by hydrothermal 
circulation. 

Observations of upper mantle anisotropy in the Pacific have been reported in many 
surface wave studies [e.g., Forsyth, 1975; Schlue and Knopoff, 1977; Tanimoto and 
Anderson, 1984, 1985; Cara and Leveque, 1988; Nishimura and Forsyth, 1988, 1989; 
Montagner and Tanimoto, 1990], and the lithospheric portion of this signal has been 
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attributed to the preferential alignment of olivine crystals with their a-axis parallel to the 
fossil direction of plate motion. Nishimura and Forsyth [1989] and Cara and Leveque 
[1988] have suggested that deep-seated (sub-lithospheric) anisotropy is related to the 
direction of present shear rather than the fossil seafloor spreading direction. In general, 
few fast-moving oceanic plates, the direction of horizontal shear in the underlying 
asthenosphere is expected to be similar to that of absolute plate motion [Hager and 
O'Connell , 1979]. 

Our study area encompasses a large range of absolute plate motions, as three plates 
are involved. Not surprisingly, surface wave studies also show a range of anisotropy 
patterns in the region. The study of Montagner and Tanimoto [1990] shows anisotropy in 
the region with the fast axis for Rayleigh waves at N80*E in the Cocos plate, 
approximately east- west in the Nazca plate, and approximately N80*W in the very eastern 
part of the Pacific Plate (within our study area). The absolute plate motion vectors from 
Gripp and Gordon [1990] have azimuths of approximately N35*E for the Cocos plate, 
N90*E for the Nazca plate, and N75*W for the Pacific plate. Our azimuthal distribution of 
data and corresponding 48 models are consistent with alignment of the olivine a axis at 
either N45*W or N45*E. The olivine a axis direction is slow for SH incident at an angle 
of 45* from the vertical (similar to the incidence angles of SS waves in this study), though 
fast for horizontally propagating P. Since most of our bounce points lie on the Cocos 
plate, it is not surprising that our results are most consistent with the absolute motion of 
this plate. 

Synthesizing the results on anisotropy from Chapter 2 and the east Pacific, there is 
not a significant pattern of anisotropy in the slow spreading north Atlantic, but there may 
be significant anisotropy in the faster spreading east Pacific. Travel time data at a better 
distribution of azimuths on the Cocos plate would help us to distinguish between upper 
mantle anisotropy and lateral heterogeneity as the mechanism of our observed azimuthal 
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distribution in travel time residuals. McKenzie [1979] andRibe [1989] have suggested 
that anisotropy in both the lithosphere and the asthenosphere will not be as pronounced 
beneath slow spreading ridges as at fast spreading ridges for several reasons. In the 
lithosphere at slow spreading rates the crystals freeze in too quickly to be affected by 
strains which at fast spreading rates can produce a preferred alignment. In the 
asthenosphere, the lattice preferred orientation of olivine crystals is suggested to become 
aligned with the flow direction when the flow is progressive simple shear, which is likely 
to occur in the asthenosphere beneath a fast spreading ridge. 

We have performed inversions of travel time, geoid, and bathymetry along a profile 
orthogonal to the Galapagos Spreading Center and just to the east of the Galapagos 
hotspot . The results of our inversions indicate that compositional variations alone are 
inadequate to produce the observed variations in travel time, geoid, and bathymetry, but 
that several different parameterizations of temperature variations produce generally good 
fits to the data (variance reductions are approximately 70%, see Table 4.3). The predicted 
excess temperature anomaly at the Galapagos hotspot ranges from 50* to 150’C. We also 
observe a long wavelength signal attributed to the slab subducting along the Middle 
America Trench. Our techniques have not been designed to deal with the strong lateral 
viscosity contrasts that a slab would presumably produce. Further work is needed to 

remove the slab signal from our geoid anomalies. 

Several authors [Watts et al ., 1985, Haxby and WeisseU 1986] have argued for the 
presence of small scale convection in the central and eastern Pacific. Haxby and Weissel 
[1986] observed linear gravity anomalies of wavelength 150 to 500 km oriented 
orthogonal to the ridge. Their observations were from an area to the south and west of our 
study area, so a direct comparison is not possible. They suggest that these anomalies 
mark the locus of longitudinal convective rolls aligned by shear imparted by a fast moving 
plate [Richter and Parsons , 1975]. A denser sampling of data along the axis of the East 
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Pacific Rise might allow us to test this theory, but our present sampling is inadequate. In 
addition, even a 500-km-wavelength signal might be smaller than we can resolve with 
teleseismic techniques. Stark and Forsyth [1983] found a periodic variation of 640- km 
wavelength in shear wave travel time residuals in the Indian Ocean, and they used this 
result to argue for the presence of linear convective rolls in the upper mantle in that region. 


CONCLUSIONS 

We have measured 600 SS-S differential travel times for paths in the East Pacific 
Rise region. The SS-S travel time residuals are only weakly correlated with seafloor age. 
The travel time - age correlation is stronger if rather than using the actual plate age we use 
the age predicted by the depth of the seafloor at the SS bounce point 

We have examined the azimuthal distribution of the SS travel time residuals, and 
although not conclusive, our results are consistent with the presence of anisotropy in our 
study area. The sense of anisotropy is consistent with the fast axis of olivine oriented 
approximately parallel to the absolute plate motion vector for the Cocos plate. It is also 
possible that lateral heterogeneity rather than azimuthal anisotropy is producing our 
observed azimuthal pattern. Sampling at a more uniform distribution of azimuths should 
make this result less ambiguous, and as more seismic stations are deployed at new 
geographic locations our chances of resolving this issue improves. 

We note several differences between the north Adantic and the east Pacific. The 
most obvious difference is that the Pacific travel time residuals are significandy larger than 
those for the Atlantic, even at a fixed age. The travel time - age relation is weaker in the 
Pacific, though this may be partially attributable to the fact that we do not sample a large 
range of plate ages in the Pacific. In the Atlantic our results are not consistent with the 
presence of a simple pattern of azimuthal anisotropy, while in the Pacific our data are 
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consistent with the presence of weak anisotropy in the upper mantle. It has been 
suggested that anisotropy may be more pronounced at fast spreading rates than at slow 
spreading rates due to a rate-dependence of the mechanism for orienting olivine grains, 
and our results are consistent with this suggestion. There is substantial ambiguity in our 
Pacific anisotropy measurements, however, due to a poor sampling of azimuths. 

We have applied a joint inversion of travel time residuals, geoid, and bathymetric 
anomalies for lateral variation in upper mantle temperature and composition to a north- 
south profile orthogonal to the Galapagos Spreading Center. We find that compositional 
variations alone are inadequate to match all of the data simultaneously. Temperature 
variations alone, however, significandy reduce the variance in all of the data. The models 
predict excess temperature in the vicinity of the Galapagos hotspot in the range of 50 to 
150'C. Further analysis is needed to remove the effects of slab structure and possible 
crustal thickening in the east Cocos plate region. 


TABLE 4.1. Digital seismograph stations used. 


Station 

Network* 

AFI 

DWWSSN 

ALQ 

DWWSSN 

ANMO 

SRO 

BOCO 

SRO 

CMB 

DWWSSN 

COL 

DWWSSN 

GAC 

CAN 

HON 

DWWSSN 

JAS 

DWWSSN 

JAS1 

DWWSSN 

KIP 

GEOSCP 

LON 

DWWSSN 

RSCP 

RSTN 

RSNT 

RSTN 

RSNY 

RSTN 

RSON 

RSTN 

RSSD 

RSTN 

SCP 

DWWSSN 

SCZ 

GEOSCP 

WFM 

GEOSCP 

ZOBO 

ASRO 


Latitude, *N 

Longitude, *E 

-13.910 

-171.777 

34.942 

-106.458 

34.946 

-106.457 

4.587 

-74.043 

38.035 

-120.385 

64.900 

-147.793 

45.70 

-75.47 

21.322 

-158.008 

37.947 

-120.438 

37.947 

-120.438 

21.420 

-158.022 

46.750 

-121.810 

35.600 

-85.569 

62.480 

-114.592 

44.548 

-74.530 

50.859 

-93.702 

44.120 

-104.036 

40.795 

-77.865 

36.600 

-121.400 

42.610 

-71.490 

-16.270 

-68.125 


t Network abbreviations as in Table 2. 1 
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TABLE 4.2 Earthquakes used in this study. 


Year 

Month 

Day 

Hour 

Min 

Sec Depth, km 

mb 

Lat., *N Long., *E 

1980 

1 

14 

21 

51 

1.8 

14 

5.6 

-33.19 -69.46 

1980 

1 

20 

1 

3 

33.2 

33 

5.0 

-22.73 -114.20 

1980 

1 

24 

19 

0 

9.5 

11 

5.3 

37.85 -121.82 

1980 

1 

27 

16 

38 

1.1 

10 

5.7 

-35.38 -105.87 

1980 

1 

28 

16 

59 

22.9 

33 

5.3 

-44.% -76.11 

1980 

2 

14 

14 

7 

1.6 

45 

5.3 

-31.77 -71.42 

1980 

3 

29 

6 

41 

50.5 

33 

5.6 

-43.08 -75.20 

X 7UV 

1980 

3 

29 

17 

49 

10.9 

33 

5.2 

-37.93 -73.50 

X 7uv 

1980 

4 

9 

8 

17 

55.1 

10 

5.5 

-31.65 -67.48 

X 7 t/V 

1980 

4 

9 

19 

56 

26.0 

10 

5.5 

-44.61 -80.10 

1980 

4 

14 

23 

57 

27.6 

100 

5.4 

-21.29 -68.73 

X 7 w 

1980 

5 

2 

19 

9 

7.8 

163 

5.4 

-24.13 -66.94 

1 70v 

1980 

5 

26 

18 

41 

36.8 

62 

6.0 

-19.36 -69.29 

1980 

5 

30 

16 

56 

22.2 

33 

5.3 

-23.11 -70.93 

X 70V 

1980 

6 

9 

3 

28 

18.9 

5 

5.6 

32.22 -114.99 

1980 

6 

11 

14 

21 

18.9 

35 

5.5 

-32.72 -71.65 

X 7 VA/ 

1980 

6 

16 

5 

45 

6.9 

87 

5.5 

-22.03 -68.47 

1980 

7 

13 

6 

20 

30.3 

103 

5.6 

-33.47 -70.15 

1980 

8 

3 

3 

0 

49.7 

151 

5.3 

-35.25 -69.99 

1980 

8 

3 

13 

42 

35.4 

10 

5.1 

-35.58 -104.63 

1980 

9 

2 

22 

18 

41.6 

10 

5.0 

-26.79 -112.92 

1980 

9 

26 

20 

26 

34.8 

10 

4.9 

-35.90 -102.94 

1980 

11 

8 

21 

35 

43.1 

103 

5.4 

-24.34 -67.65 

1 70v 

1980 

12 

11 

18 

15 

3.5 

80 

6.1 

-21.27 -68.15 

1980 

12 

20 

20 

26 

47.2 

571 

5.2 

-24.37 -63.42 

1981 

1 

7 

16 

26 

42.5 

37 

5.5 

-23.75 -70.62 

1981 

1 

7 

20 

31 

12.1 

10 

5.1 

-49.91 -114.14 

X 70 X - 

1981 

3 

2 

21 

19 

56.0 

246 

5.1 

-22.64 -65.99 

1981 

3 

23 

19 

28 

10.7 

46 

5.8 

-33.66 -71.89 

X 70 X 

1981 

3 

26 

18 

4 

44.7 

138 

5.8 

-19.37 -68% 

1981 

4 

1 

18 

3 

36.5 

554 

5.9 

-27.31 -63.32 

1981 

4 

6 

14 

34 

1.3 

88 

5.3 

-35.37 -71.05 

1981 

4 

16 

22 

5 

53.0 

33 

5.1 

-20.16 -70.70 

1981 

5 

17 

17 

4 

58.0 

33 

5.4 

-27.22 -71.83 

170 1 

1981 

5 

24 

11 

44 

3.5 

10 

5.2 

-22.16 -114.20 

X TO x 

1981 

6 

1 

19 

50 

13.3 

330 

5.0 

-20.49 -65.19 

1981 

6 

4 

7 

39 

40.4 

33 

5.1 

-34.20 -78.83 

X7u X 

1981 

6 

10 

4 

6 

10.2 

10 

4.8 

-35.80 -102.17 

X TO X 

1981 

6 

16 

5 

41 

49.1 

76 

5.4 

-21.49 -68.34 

X TO X 

1981 

6 

21 

10 

30 

1.1 

36 

5.2 

-20% -70.45 

1981 

6 

22 

17 

53 

21.3 

24 

5.1 

-13.17 -74.52 

I TO X 

1981 

7 

10 

18 

2 

4.9 

10 

5.2 

-37.22 -95.37 

X TO X 

1981 

7 

18 

11 

15 

18.1 

246 

5.0 

-22.68 -66.24 

1 70 1 

1981 

7 

28 

3 

3 

21.6 

43 

5.5 

-41.57 -73.20 

X TO X 

1981 

8 

17 

2 

18 

59.9 

37 

5.5 

14.52 -93.77 

1981 

8 

21 

22 

52 

40.6 

10 

5.1 

-26.51 -114.76 

X TO X 

1981 

8 

22 

23 

47 

41.4 

10 

5.2 

-35.83 -103.30 

X TO X 

1981 

9 

6 

16 

43 

19.9 

10 

5.4 

-36.17 -100.70 
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1981 

9 

19 

11 

40 

56.9 

30 

5.6 

- 39.08 

- 74.81 

1981 

9 

20 

10 

48 

20.3 

234 

5.1 

- 23.08 

- 66.63 

1982 

1 

21 

21 

52 

41.0 

10 

5.5 

19.22 

- 155.60 

1982 

1 

21 

22 

29 

13.0 

13 

5.5 

19.20 

- 155.57 

1982 

2 

10 

20 

00 

ro 

3.5 

226 

5.9 

- 22.79 

- 66.66 

1982 

3 

22 

21 

44 

9.3 

126 

5.1 

- 19.88 

- 68.74 

1982 

4 

18 

11 

30 

55.3 

10 

5.9 

- 28.16 

- 114.07 

1982 

4 

19 

11 

54 

27.5 

250 

5.3 

- 22.60 

- 66.02 

1982 

5 

3 

7 

26 

2.5 

107 

5.5 

- 23.47 

- 68.81 

1982 

5 

17 

2 

39 

26.1 

10 

5.0 

- 35.35 

- 103.84 

1982 

5 

27 

1 

39 

19.5 

10 

5.3 

- 13.38 

- 111.62 

1982 

5 

29 

3 

26 

58.7 

116 

5.2 

- 22.35 

- 68.37 

1982 

5 

29 

17 

22 

24.4 

47 

5.1 

- 24.99 

- 70.25 

1982 

6 

1 

4 

14 

15.3 

32 

6.1 

- 41.63 

- 74.99 

1982 

6 

5 

0 

17 

0.3 

33 

5.4 

- 43.10 

- 75.21 

1982 

6 

15 

1 

20 

35.2 

10 

5.2 

- 36.52 

- 98.65 

1982 

7 

11 

2 

13 

37.7 

39 

5.3 

- 16.65 

- 73.21 

1982 

8 

5 

9 

16 

41.3 

40 

5.4 

- 26.68 

- 70.66 

1982 

8 

23 

21 

35 

5.9 

10 

5.1 

- 29.04 

- 111.78 

1982 

8 

27 

3 

56 

49.0 

10 

5.5 

- 35.52 

- 104.61 

1982 

9 

9 

6 

40 

21.4 

10 

5.0 

- 35.55 

- 102.42 

1982 

9 

26 

4 

2 

6.6 

10 

5.2 

- 42.70 

- 82.26 

1982 

10 
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- 49.92 - 113.59 

5.4 

- 45.71 

- 75.99 

5.1 

- 27.94 

- 67.06 

6.5 

- 24.49 

- 70.17 

6.1 

0.13 

- 77.67 

6.5 

0.15 

- 77.83 

5.4 

- 38.89 

- 92.19 
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1987 

3 

15 

6 

3 

3.1 

33 

5.1 

- 24.30 

- 70.55 

1987 

3 

22 

3 

23 

57.7 

42 

5.9 

- 24.08 

- 70.08 

1987 

3 

23 

20 

47 

57.2 

33 

4.9 

- 34.07 

- 72.23 

1987 

4 

8 

17 

42 

36.0 

49 

5.1 

11.66 

- 86.36 

1987 

5 

10 

15 

16 

20.2 

160 

5.6 

- 30.91 

- 65.42 

1987 

5 

12 

16 

12 

42.1 

121 

5.4 

- 21.79 

- 68.32 

1987 

5 

14 

16 

4 

26.5 

29 

5.6 

- 5.62 

- 81.39 

1987 

5 

15 

13 

49 

14.2 

10 

5.5 

- 49.98 

- 115.23 

1987 

5 

19 

12 

56 

25.0 

39 

5.6 

- 30.33 

- 71.59 

1987 

5 

23 

5 

1 

8.1 

41 

5.2 

- 20.19 

- 70.71 

1987 

6 

9 

22 

46 

11.7 

10 

5.2 

- 35.24 

- 106.68 

1987 

6 

11 

5 

8 

52.8 

33 

5.0 

- 20.51 

- 70.61 

1987 

6 

15 

6 

23 

8.9 

44 

5.5 

- 12.46 

- 76.73 

1987 

6 

19 

19 

0 

8.3 

119 

5.4 

- 21.37 

- 68.38 

1987 

6 

27 

9 

9 

5.5 

61 

5.9 

- 14.13 

- 76.08 

1987 

7 

6 

1 

6 

7.5 

10 

6.3 

- 27.02 

- 108.25 

1987 

7 

8 

11 

50 

14.6 

10 

6.1 

- 26.97 

- 108.20 

1987 

7 

9 

4 

7 

34.3 

64 

5.2 

- 20.48 

- 68.74 

1987 

7 

15 

7 

16 

13.6 

67 

6.0 

17.56 

- 97.18 

1987 

7 

17 

1 

57 

38.7 

10 

5.4 

1.52 

- 85.30 

1987 

7 

21 

13 

27 

13.2 

10 

5.2 

- 36.29 

- 97.24 

1987 

7 

22 

8 

23 

3.7 

98 

5.2 

15.76 

- 93.44 

1987 

7 

29 

15 

22 

33.0 

141 

5.4 

- 28.87 

- 67.17 

1987 

8 

4 

15 

4 

40.4 

38 

5.8 

- 40.46 

- 73.10 

1987 

8 

6 

15 

15 

34.7 

10 

5.6 

- 5.40 

- 105.02 

1987 

8 

8 

15 

48 

57.4 

82 

6.4 

- 19.19 

- 70.14 

1987 

8 

9 

8 

24 

20.0 

10 

5.0 

- 35.23 

- 103.95 

1987 

8 

13 

15 

23 

11.4 

74 

6.1 

- 17.83 

- 71.08 

1987 

8 

15 

9 

34 

47.7 

104 

5.1 

- 23.10 

- 68.57 

1987 

8 

15 

18 

4 

22.1 

32 

6.3 

- 28.15 

- 70.89 

1987 

8 

21 

23 

0 

52.5 

10 

5.3 

23.85 

- 108.77 

1987 

8 

24 

6 

9 

43.2 

33 

5.1 

- 20.15 

- 70.58 

1987 

8 

26 

13 

1 

23.0 

10 

5.1 

- 37.27 

- 95.22 

1987 

9 

8 

2 

58 

52.8 

25 

5.4 

6.52 

- 82.55 

1987 

9 

11 

0 

34 

51.8 

128 

5.4 

- 22.32 

- 68.36 

1987 

9 

11 

4 

2 

3.1 

95 

5.1 

- 31.46 

- 70.73 

1987 

9 

18 

8 

43 

25.9 

71 

5.6 

- 24.41 

- 69.10 

1987 

9 

19 

9 

21 

41.7 

35 

5.1 

- 14.97 

- 75.66 

1987 

9 

19 

21 

18 

25.9 

69 

5.7 

- 9.19 

- 79.07 

1987 

9 

22 

13 

43 

39.7 

23 

6.1 

- 0.96 

- 78.09 

1987 

9 

22 

16 

21 

38.5 

33 

5.8 

- 1.08 

- 78.11 

1987 

10 

3 

3 

35 

6.1 

119 

5.8 

- 18.22 

- 69.31 

1987 

10 

4 

8 

15 

17.8 

45 

5.4 

10.74 

- 85.94 

1987 

10 

6 

14 

39 

54.6 

10 

5.2 

- 35.78 

- 101.02 

1987 

10 

7 

0 

51 

37.0 

117 

5.6 

- 23.07 

- 68.12 

1987 

10 

20 

21 

10 

59.9 

10 

5.3 

0.91 

- 87.09 

1987 

10 

21 

6 

11 

33.3 

10 

5.4 

- 23.10 

- 114.48 

1987 

10 

21 

23 

25 

51.7 

76 

5.4 

- 21.14 

- 69.69 

1987 

10 

27 

21 

58 

17.7 

616 

5.9 

- 28.73 

- 62 .% 

1987 

10 

31 

5 

11 

17.0 

66 

5.1 

- 25.88 

- 69.80 

1987 

10 

31 

8 

52 

30.8 

33 

5.4 

- 6.98 

- 80.57 

1987 

11 

6 

18 

47 

34.9 

538 

5.8 

- 22.84 

- 63.63 

1987 

11 

9 

17 

46 

21.0 

62 

5.4 

- 22.07 

- 69.47 

1987 

11 

17 

3 

40 

8.9 

76 

5.8 

12.57 

- 86.98 

1987 

12 

3 

11 

4 

32.1 

60 

5.6 

- 21.40 

- 68.26 

1987 

12 

8 

14 

47 

59.3 

10 

5.6 

- 32.59 

- 112.16 

1987 

12 

28 

8 

16 

21.3 

52 

6.0 

- 28.08 

- 70.61 
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1988 

1 

19 

1988 

2 

5 

1988 

2 

6 

1988 

2 

22 

1988 

3 

28 

1988 

3 

30 

1988 

4 

12 

1988 

5 

5 

1988 

5 

5 


7 

30 

29.7 

14 

1 

2.2 

18 

3 

53.6 

19 

13 

17.4 

18 

36 

26.4 

23 

50 

56.2 

23 

19 

57.3 

10 

4 

13.5 

22 

32 

48.1 


21 

6.3 

- 24.75 

33 

6.2 

- 24.77 

283 

6.0 

- 18.04 

71 

5.9 

- 20.93 

10 

5.7 

- 36.03 

37 

5.8 

- 24.91 

54 

6.1 

- 17.32 

10 

6.1 

- 26.91 

10 

5.5 

- 26.76 


- 70.60 

- 70.37 

- 66 .% 

- 69.80 

- 102.81 

- 70.41 

- 72.40 

- 113.39 

- 113.70 



TABLE 4.3. Inversion Models and Variance Reduction 


Model: Temperature variations only 


Variance reduction, % 

Layer 

Viscosity 






thickness 

structure 

AT range 

total 

bathvm geoid SS-S 

0-150 km 

cvm 

415 K 

71 

57 

72 

76 

0-150 km 

lvz 

560 K 

72 

59 

51 

92 

0-300 km 

cvm 

175 K 

69 

55 

75 

70 

0-300 km 

lvz 

285 K 

66 

63 

28 

96 

0-650 km 

cvm 

85 K 

70 

47 

76 

74 

0-650 km 

lvz 

155 K 

70 

54 

46 

93 

Model: Compositional variations only 

Variance reduction, % 

Layer 

Viscosity 






thickness 

structure 

AM g# range 

total 

bathvm 

geoid SS-S 

0-150 km 

cvm 

1.7 

35 

50 

90 

-12 

0-150 km 

lvz 

2.7 

23 

69 

52 

-15 

0-300 km 

cvm 

0.7 

37 

44 

95 

-9 

0-300 km 

lvz 

2.0 

23 

23 

44 

7 

0-650 km 

cvm 

0.3 

37 

34 

96 

-6 

0-650 km 

lvz 

1.8 

35 

-6 

69 

25 

Model: Thermal and compositional variations in same layer 

Variance reduction, % 

Layer 

Viscosity 

AT AMg# 





thickness 

structure 

range range 

total 

bathvm 

geoid 

SS-S 

0-150 km 

cvm 

525 K 1.9 

89 

57 

90 

100 

0-150 km 

lvz 

555 K 1.9 

81 

78 

58 

100 

0-300 km 

cvm 

275 K 1.1 

89 

51 

94 

100 

0-300 km 

lvz 

290 K 1.9 

75 

87 

37 

100 

0-650 km 

cvm 

140 K 0.6 

89 

43 

95 

100 

0-650 km 

lvz 

160 K 1.5 

84 

72 

69 

100 

Model: Temperature variations only, high-order geoid 

Variance reduction, % 

Layer 

Viscosity 






thickness 

structure 

AT range 

total 

bathvm 

geoid SS-S 

0-150 km 

cvm 

80 K 

30 

18 

62 

27 

0-150 km 

lvz 

275 K 

47 

52 

5 

67 

0-300 km 

cvm 

35 K 

29 

12 

78 

22 

0-300 km 

lvz 

195 K 

43 

59 

-10 

55 

0-650 km 

cvm 

20 K 

30 

10 

83 

26 

0-650 km 

lvz 

80 K 

23 

31 

-46 

56 


cvm = constant viscosity mantle 
lvz = mantle with low viscosity zone 
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Figure 4. 1 . 


Figure 4.2. 


Figure 4.3. 


Figure 4.4. 


Figure 4.5. 


Figure Captions 

Location map of the East Pacific Rise, after Macdonald [1989]. Spreading 
centers, transform faults, propagating rifts, and overlapping spreading 
centers are shown. 

Distribution of earthquakes (triangles) and seismograph stations (circles) 
used to measure SS-S differential travel times. Stations are from the GDSN 
and GEOSCOPE digital arrays. Earthquakes are from the Harvard CMT 
catalogue (generally mb > 5.0) from the years 1980-1988, and from a listing 
of intraplate events from the years 1980-1987 (E. A. Bergman, personal 
communication, 1991). Plate boundaries are from DeMets et al. [ 1990]. 
Lambert equal area projection with pole of projection at 5'N, - 1 10*W. 

Bathymetric map of East Pacific Rise and Cocos plate region, from 
DBDB5C [US. Naval Oceanographic Office , 1985]. Contour interval 
1000 m. Straight lines show location of profiles examined in this study. 

Residual geoid map for East Pacific Rise and Cocos plate region, derived 
from gridded altimetric data of Marsh et al. [1986] with the low order (to 
degree and order 7, tapered to degree and order 11) portion of GEM-T2 
reference field [Marsh et al., 1990] removed. Contour interval 2 m. 

Straight lines show location of profiles examined in this study. 

Map view of SS-S residuals relative to PREM [Dziewonski and Anderson, 
1981], corrected for Earth ellipticity and seafloor bathymetry. Residuals are 



174 


plotted at the SS bounce point. The size of each symbol scales linearly with 
the magnitude of the residual. A mean of +4. 15 s has been removed from 
the data for plotting purposes. Negative residuals indicate early SS (or late 
S). Plate boundaries are from DeMets et al. [1990]. Mercator projection. 

Figure 4.6. (a) SS-S travel time residual versus square root of seafloor age. Each 

plotted point represents the weighted mean of 25 adjacent data points. 
Horizontal and vertical bars are standard error s of the means of the travel 
time residuals and (age) 1 / 2 . Linear regression yields a slope of -0.01 ± 

0.08 s / (My) 1 / 2 for a 0- 50 My age range (solid line). Age dependence of 
travel time residual predicted by the plate cooling model is shown as a 
dashed line. 

(b) SS-S travel time residual versus the square root of the age predicted by 
the plate cooling model [Parsons and Sclater, 1977] for the given depth. 
Each plotted point represents the weighted mean of 25 adjacent data points. 
Linear regression yields a slope of -0.26 ± 0.04 s / (My) 1 / 2 for a 0-33 My 
age range (solid line). Age dependence of travel time residual predicted by 
the plate cooling model is shown as a dashed line. 

Figure 4.7. (a) SS-S residual versus azimuth 9. Each plotted point represents the 

weighted mean of 25 adjacent data points. The solid curve shows the best- 
fitting 40 variation derived from these data. 

(b) Age-corrected SS-S residual (see text) versus azimuth 6. Each plotted 
point represents the weighted mean of 25 adjacent data points. The solid 
curve shows the best-fitting 40 variation derived from these data. 
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Figure 4.8. 


Figure 4.9. 


Figure 4. 10. 


Figure 4.11. 


Map view of sampling azimu ths. Lines indicate the wave path azimuth at 
the SS bounce point Mercator projection. 

Map view of age-corrected SS-S residuals. Residuals are plotted at the SS 
bounce point The size of each symbol scales linearly with magnitude of the 

residual. 

Profile 1 across the East Pacific Rise south of the Pacific-Cocos-Nazca 
triple junction, (a) Seafloor age at SS bounce points. Data corresponding 
to SS bounce points of “A” and “B” quality are designated by crosses. “C” 
quality data are designated by small dots. The rise axis is marked by a 
straight line, (b) Bathymetry at SS bounce points (see Figure 4.3). (c) 
Residual geoid at SS bounce points (see Figure 4.4). (d) SS-S travel time 
residuals, (e) Age corrected geoid, using the method of Hager [1983]. (0 
Age corrected geoid, using the method of Parsons and Richter [1980] (g) 
Age corrected bathymetry, using the method of Parsons and Sclater [1977]. 
(h) Age corrected SS-S travel time residuals (see text), (i) High-order 
portion of the geoid (GEM-T2 reference field to degree and order 10 
removed). No age correction applied. The residuals shown in b - i are 
weighted moving averages (such that each point is used twice) of 4 adjacent 

data points. 

Profile 2, across the Galapagos Spreading Center. See Figure 4.10 for 
explanation. Location nearest present Galapagos hotspot marked with an 


asterisk. 
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Figure 4. 12. Profile 3 across the East Pacific Rise north of the Pacific-Cocos-Nazca triple 
junction. See Figure 4.10 for explanation. 

Figure 4. 13. Profile 4 across the East Pacific Rise north erf the Pacific-Cocos-Nazca triple 
junction. See Figure 4.10 for explanation. 

Figure 4. 14. Results of combined inversion of geoid, bathymetry, and SS-S travel time 
residuals for upper mantle temperature variations along Profile 2. The 
viscosity structure is taken to consist of a 40-km-thick high-viscosity lid 
overlying a constant- viscosity halfspace. Location nearest Galapagos 
hotspot marked with an asterisk. 

(a) Three solutions for lateral temperature variations: Dotted line: 
Temperature perturbations constrained to be uniform over 0-150 km depth. 
Long-dashed line: Temperature perturbations constrained to be uniform 
over 0-300 km depth. Short-dashed line: Temperature perturbations 
constrained to be uniform over 0-650 km depth. 

(b) Observed (solid line) and predicted profiles of SS-S travel time residual. 
The “observed” profile is actually a filtered version of the observations, 
containing only the wavelengths used in the inversion (850 to 5000 km). 

Line types correspond to those of the temperature models. 

(c) Observed and predicted geoid profiles. Same treatment as in (b). 

(d) Observed and predicted bathymetry profiles. Same treatment as in (b). 

Figure 4.15. 


Results of combined inversion of geoid, bathymetry, and SS-S travel time 
residuals for upper mantle temperature variations along Profile 2. The 
viscosity structure includes a zone extending from the base of the high- 
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viscosity lid to a depth of 200 km with a viscosity equal to 0.01 that of the 
underlying mantle. Other aspects as in Figure 4.14. 

Figure 4. 16. Results of combined inversion of geoid, bathymetry, and SS-S travel time 
residuals for variations in upper mantle composition (Mg#). The viscosity 
structure is taken to consist of a 40-km- thick high-viscosity lid overlying a 
constant-viscosity halfspace. 

(a) Three solutions for lateral composition variations: Dotted line: 
Composition perturbations constrained to be uniform over 0-150 km depth. 
Long-dashed line: Composition perturbations constrained to be uniform 
over 0-300 km depth. Short-dashed line: Composition perturbations 
constrained to be uniform over 0-650 km depth. 

(b) Observed (solid line) and predicted profiles of SS-S travel time residual. 

(c) Observed and predicted geoid profiles. 

(d) Observed and predicted bathymetry profiles. 

Figure 4.17. Results of combined inversion of geoid, bathymetry, and SS-S travel time 
residuals for variations in upper mantle composition (Mg#). The viscosity 
structure includes a zone extending from the base of the high-viscosity lid to 
a depth of 200 km with a viscosity equal to 0.01 that of the underlying 
mantle. Other aspects as in Figure 4.16. 

Figure 4. 18. Results of combined inversion of geoid, bathymetry, and SS-S travel time 
residuals for both upper mantle temperature and composition variations. 

The viscosity structure is taken to consist of a 40-km-thick high-viscosity 
lid overlying a constant-viscosity halfspace. 
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(a) Three solutions for lateral temperature variations: Dotted line: 
Composition perturbations constrained to be uniform over 0- 150 km depth. 
Long-dashed line: Composition perturbations constrained to be uniform 
over 0-300 km depth. Short-dashed line: Composition perturbations 
constrained to be uniform over 0-650 km depth. 

(b) Three solutions for lateral composition variations: Dotted line: 
Composition perturbations constrained to be uniform over 0- 150 km depth. 
Long-dashed line: Composition perturbations constrained to be uniform 
over 0-300 km depth. Short-dashed line: Composition perturbations 
constrained to be uniform over 0-650 km depth. 

(c) Observed (solid line) and predicted profiles of SS-S travel time residual. 

(d) Observed and predicted geoid profiles. 

(e) Observed and predicted bathymetry profiles. 

Figure 4. 19. Results of combined inversion of geoid, bathymetry, and SS-S travel time 
residuals for both upper mantle temperature and composition variations. 

The viscosity structure includes a zone extending from the base of the high- 
viscosity lid to a depth of 200 km with a viscosity equal to 0.01 that of the 
underlying mantle. Other aspects as in Figure 4. 1 8. 

Figure 4.20. Results of combined inversion of high-order geoid (see Figure 4. lOi), 
bathymetry, and SS-S travel time residuals for upper mantle temperature 
variations along Profile 2. The viscosity structure is taken to consist of a 
40-km-thick high-viscosity lid overlying a constant-viscosity halfspace. 
Other aspects as in Figure 4.14. 
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Figure 4.2 1 . Results of combined inversion of high-order geoid (see Figure 4. lOi), 
bathymetry, and SS-S travel time residuals for upper mantle temperature 
variations along Profile 2. The viscosity structure includes a zone extending 
from the base of the high- viscosity lid to a depth of 200 km with a viscosity 
equal to 0.01 that of the underlying mantle. Other aspects as in Figure 


4.15. 
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Figure 4.5 
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Figure 4.9 
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Chapter 5 


Conclusions 


In this thesis we have examined lateral heterogeneity in the lithosphere and 
asthenosphere beneath oceanic regions through a combined analysis of seismic body 
wave differential travel times and attenuation, geoid anomalies, and residual depth 
anomalies. Our focus has been to assess variations at the 1000 km scale, comparable to 
the thickness of the upper mantle. We have examined the data to determine whether 
convection on this scale can produce observable temperature and compositional 
anomalies. We have developed an inversion technique which allows shear wave travel 
time, geoid, and bathymetry data to be combined and jointly inverted for lateral variations 
in upper mantle temperature and composition structure, and we have applied it to data sets 

from both the north Atlantic and east Pacific. 

Several general conclusions can be drawn from our study of SS-S differential travel 
times. In the Atlantic, the SS-S travel time residual decreases linearly with square root of 
age, in general agreement with the plate cooling model to an age of 80- 100 My [Parsons 
and Sclater, 1977]. In the Pacific, SS-S travel time residuals are only weakly correlated 
with seafloor age, although this may be partially atrributable to the fact that we do not 
sample a large range of plate ages in the Pacific. The travel time - age correlation in the 
Pacific is stronger if rather than using the actual plate age we use the age predicted by the 
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depth of the seafloor at the SS bounce point, in the same sense as that predicted by the 
lithospheric rejuvenation hypothesis [Menard and McNutt , 1982]. The Pacific travel time 
residuals are significandy larger than the Adantic data, even at a fixed age. 

Azimuthal anisotropy is not clearly resolved in the Atlantic but may be present in the 
east Pacific. Although not conclusive, our east Pacific data are consistent with the 
presence of anisotropy, with the fast axis of olivine oriented approximately parallel to the 
absolute plate motion vector for the Cocos plate. It has been suggested that anisotropy 
may be more pronounced at fast spreading rates than at slow spreading rates both in the 
lithosphere (due to rate-dependence of the mechanism for orienting olivine crystals) and 
in the asthenosphere (due to the flow regime), and our results are consistent with this 
suggestion. There is substantial ambiguity in our Pacific anisotropy measurements, 
however, due to a poor sampling of azimuths. 

We have observed lateral variations in age-corrected SS-S residuals with a 
dominant wavelength of 1000 to 2000 km. Comparable variations are also observed in 
bathymetry and geoid height. We have formulated a joint inversion of geoid, travel time, 
and bathymetry for lateral variations in temperature and composition as a means to assess 
more quantitatively the observed correlations. In the north Atlantic, inversion for 
temperature favors the presence of an upper mantle low viscosity zone and temperature 
anomalies concentrated at depths less than 300 km. We are unable to match travel time 
residuals simultaneously with geoid and bathymetry solely with lateral variations in bulk 
composition (Mg#). Joint inversions for temperature and composition provide good fits 
to both travel time and and geoid regardless of viscosity structure or layer depth and 
thickness, but the best fits to bathymetry come from models with a low- viscosity zone 
and thermal or compositional variations confined to shallow depth. The variations in 
Mg# predicted in the joint inversion for temperature and composition are comparable to 
those found by Michael and Bonatti [1985] in a study of dredged peridotites along the 
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Mid-Atlantic Ridge and may be related to variations in melt production along the ridge. In 
the Pacific we applied the inversion techniques to a north-south profile orthogonal to the 
Galapagos spreading center. In this region we also find that compositional variations 
alone are inadequate to produce good fits to travel time and geoid and bathymetry 
simultaneously. Temperature variations alone, however, produce significant variance 
reduction. The models predict excess temperature in the vicinity of the Galapagos hotspot 
in the range 50 to 150 K. 

As a complement to the travel time studies, we also measured SS-S differential 
attenuation in the north Atlantic region. As with the travel time residuals, we find that 8t* 
decreases with increasing seafloor age. We do not observe evidence for a narrow region 
of very low Q lying beneath the Mid-Atlantic Ridge as reported by Sobmon [1973], 
although this may be due to the limited spatial resolution of the long-period shear waves 
mtized in our study. We derive empirical Q' 1 - temperature relations by predicting St* as a 
function of lithospheric age by means of the plate cooling model [Parsons and Sclater, 
1977]. We find that our St* variations versus age are well fit with a simple Q-temperature 
relation, but the relation generally predicts smaller variations in Q" 1 for a given 
temperature change than the laboratory-derived Q-temperature relation of Sato and Sacks 
[1989]. Along-axis St* profiles are constructed by combining the temperature variations 
along the axis of the Mid- Atlantic Ridge derived in Chapter 2 and the empirical Q- 
temperature relation derived for the lithosphere. 

The techniques developed in this thesis should be of general utility for both global 
and regional studies of the upper mantle. The Hawaiian hotspot would be an interesting 
new area in which to apply these inversion techniques. Mid-plate swells have been 
hypothesized to result from several effects, including elevated temperature in the 
convecting mantle [Sandwell and Poehls, 1980], thermal expansion within the conducting 
portion of the lithosphere [Crough, 1978; Detrick and C rough, 1978; Menard and 
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McNutt, 1982], and underplating of the lithosphere with depleted mantle material 
[Jordan, 1979]. Each of these processes involves a different source of the anomaly 
causing the swell, and each predicts a correspondingly different geophysical signature. 
The methods proposed here, which distinguish between thermal and compositional 
mechanisms, offer great promise for distinguishing between these competing models. 
Woodward and Masters [1991], for instance, have observed negative SS-S travel time 
residuals with SS bounce points in the vicinity of the Hawaiian swell, a result which 
would favor a lateral variation in composition in the upper mantle beneath the swell. 

Application of these techniques and measurements to a global travel time data set 
will allow general constraints to be placed on upper mantle viscosity structure, the 
presence of partial melt, and the mechanism of lateral heterogeneity. The viscosity 
structure of the mantle is poorly known yet plays a key role in models of mantle 
convection. The viscosity structures we employ in this thesis are quite simple but have 
been chosen to represent two models widely invoked in other studies - a constant or 
nearly constant viscosity mantle [e.g., Peltier, 1989] and a mantle with a thin low 
viscosity layer [e g., Craig and McKenzie, 1986; Robinson et al., 1987]. The viscosity 
structure of the Earth may be temperature and pressure -dependent or vary laterally, but 
we have not considered viscosity structures of this type [e.g., Revenaugh and Parsons, 
1987]. Much work remains to be done to determine ways to incorporate lateral variations 
in viscosity into quantitative treatments of these problems. 

Several other improvements could be made in future studies of the type presented 
here. Our models thus far have been limited to simply parameterized one-dimensional 
variations in temperature and composition within a single layer. It is likely that these 
lateral variations are not constant within a given layer and that there are two-dimensional 
lateral variations independent of lithospheric aging. The techniques outlined in this thesis 
can be generalized to a multilayer system and to two-dimensional wavenumber, although 
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wc feel that an areally denser distribution of seismic data would be necessary to justify 
these more complicated models. Kernels for seismic surface waves are strongly peaked 
in the upper mantle, and such data would provide a useful constraint independent of 
differential body wave travel times. The inclusion of surface wave data would help to 
distinguish between lithospheric and asthenospheric effects and may allow for two or 
more layers to be independently resolved. Extension erf the modelling to three 
dimensions would permit an assessment of the degree to which mantle anomalies beneath 
the ridge extend off axis. Implicit in our age-correction is the assumption that the 
anomalous properties of the ridge mantle are steady state on a time scale of tens of 
millions of years. Examination of the cross-axis profiles in the east Pacific suggests that 
this assumption might not always hold. 

We have not modelled the effects of partial melting which could accompany the 
temperature variations we predict. The effect of retained melt on the physical properties 
of the mantle depends critically on the melt fraction and geometry, characteristics 
presently poorly known. Sato et al. [1988, 1989] downplay the importance of partial 
melt and suggest that most mantle seismic velocity anomalies can be explained by 
temperature variations at subsolidus conditions. Our attenuation measurements can also 
be explained with a thermal mechanism. Combined analysis erf compressional as well as 
shear wave differential travel times and attenuation may help to resolve whether 
significant partial melting is required to explain the travel time residuals in the oceanic 
upper mantle. 

Further analysis of the East Pacific Rise data set is warranted to refine the 
conclusions from that study. We hope to improve our azimuthal distribution of data by 
obtaining seismograms from new stations which have begun operating only recently. 
Since we have gathered data only through 1987, there is currently considerable data to 
which we have not yet obtained access. Sampling at a more uniform distribution of 
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azimuths should make our conclusions regarding anisotropy in the region less 
ambiguous, and as more seismic stations are deployed at new locations our chances of 
resolving this issue improves. Additional data would also be of use in extending the 
analysis of age dependence to older lithosphere, and would be of use in attempts to assess 
further why the upper mantle in the east Pacific region is slower than in the north Atlantic. 

The data from the Cocos plate show interesting small-scale ( X from 800 to 1200 
km) patterns which might be best examined in two dimensions. We may find it necessary 
to refine the numerical techniques for obtaining kernels at die smaller wavelengths, 
however, as we employ a predictor-corrector algorithm which becomes unstable when the 
derivative of the kernels approaches zero (as is the case for geoid kernels at short 
wavelength for a mantle with a thin low viscosity zone). 

Other complicating factors for the East Pacific Rise study include the effect of the 
subdue tion zone along the Middle America Trench and possible crustal thickness 
variations in the eastern Cocos plate. Crustal thickness variations can be removed from 
the bathymetry and the geoid by simple isostatic mass balance, and from the travel time 
residuals by calculating the additional delay due to excess crust relative to an equivalent 
thickness of mantle. The predicted slab effect can presumably be removed from the low- 
order geoid in the manner of Hager [1984]. The effect of the slab on the differential 
travel times is more difficult to assess. 

An attenuation study in the east Pacific would be a natural extension of this 
work, as the relevant seismograms have already been collected for the travel time study in 
Chapter 4. From their study of Rayleigh wave phase propagation, Canos and Mitchell 
[1981] found that the east Pacific region has lower Q than the Atlantic. An analysis of 
body wave differential attenuation, with horizontal resolution superior to that achievable 
with surface waves, would constitute an interesting complement to our travel time 
analysis. In another study Schlue [1981] examined differential attenuation between the S 
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and SS phases with SS bounce points in a small area close to our Profile 1 of Chapter 4. 
The results of the Schlue [1981] suggested that Q- 1 does not decrease with age, at least in 
the small range of plate ages (0-15 My) that were included in his study. 

In conclusion, we feel that the procedures developed in this thesis represent an 
important first step towards the discrimination of the global extent and mechanism of 
lateral heterogeneity in the upper mantle on a regional (1000 to 4000 km) scale. 
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Appendix II. SS-S Relative Attenuation (5t*) Values for the North Atlantic Region 


Origin time 


Stn 

Yr Mo 

Da Hr Mn Sec 

GRFO 

80 

6 

25 12: 4:56.9 

GRFO 

80 

7 

26 12:53:40.7 

ZOBO 

80 

10 

10 12:25:23.5 

ZOBO 

80 

10 

10 15:39: 9.8 

ZOBO 

80 

10 

13 6:37:39.2 

GRFO 

80 

11 

26 17:35:39.1 

BOCO 

80 

12 

5 13:32: 5.9 

KONO 

81 

4 

22 23:16:54.1 

SCP 

81 

7 

7 21:10:57.7 

SCP 

81 

8 

13 2:58:11.2 

GRFO 

81 

10 

18 4:31:2.7 

KEV 

81 

11 

28 19:24: 5.6 

KONO 

81 

12 

6 14:54:29.6 

SCP 

82 

1 

3 14: 9:50.4 

BER 

82 

1 

30 2:35:11.0 

TOL 

82 

4 

6 19:56:53.7 

ANMO 

82 

4 

10 9:47: 4.4 

BOCO 

82 

5 

2 7:12:44.5 

KONO 

82 

5 

10 1:25:56.8 

SCP 

82 

8 

12 8:41:50.5 

ANMO 82 

8 

12 8:41:50.5 

GAC 

82 

10 

1 5:35:47.3 

ZOBO 

82 

12 

22 15:40:13.3 

GAC 

83 

1 

17 12:41:29.3 

RSNT 

83 

1 

17 12:41:29.3 

RSON 

83 

1 

17 12:41:29.3 

COL 

83 

1 

17 12:41:29.3 

RSCP 

83 

1 

17 12:41:29.3 

GAC 

83 

1 

19 0:2:13.6 

RSCP 

83 

1 

19 0: 2:13.6 

RSNY 

83 

3 

19 21:41:44.3 

GAC 

83 

3 

19 21:41:44.3 

RSCP 

83 

3 

19 21:41:44.3 

RSCP 

83 

3 

24 4:17:31.6 

BER 

83 

3 

31 13:12:52.3 

KONO 

83 

4 

10 1:16: 5.8 

GRFO 

83 

4 

11 8:18:10.3 

BCAO 

83 

4 

11 8:18:10.3 

TOL 

83 

5 

9 15:53:2.7 

KONO 

83 

5 

12 10:51:49.6 

KEV 

83 

5 

12 10:51:49.6 

RSCP 

83 

5 

14 23:13:46.5 

GRFO 

83 

5 

20 17:40:14.4 

GAC 

83 

7 

14 2:54:18.6 

RSNY 

83 

8 

6 15:43:52.6 

RSCP 

83 

8 

6 15:43:52.6 

SCP 

83 

10 

13 13: 6:38.2 

RSSD 

83 

10 

13 13:6:38.2 

GRFO 

83 

10 

13 23:27:33.2 

RSON 

83 

10 

13 23:27:33.2 

SCP 

83 

10 

21 20:34:49.1 

SCP 

83 

10 

30 4:12:27.7 

TOL 

83 

11 

22 14:20:58.8 


Epicenter 

Dep mb Lat Long Dist 

151.1 5.7 4.70 -75.35 84.2 

15.0 5.2 7.47 -33.89 56.4 

10.1 6.5 36.25 1.36 83.8 

10.0 6.2 35.72 1.31 83.5 

15.0 5 2 36.53 2.07 84.4 

40.0 5.0 7.96 -72.64 80.0 

15.0 5.0 35.87 1.68 75.8 

15.0 5.1 7.25 -36.06 62.6 

10.0 5.7 -0.07 -18.59 67.2 

18.0 5.5 44.86 17.33 65.9 

26.0 5.4 8.45 -72.81 79.8 

15.0 5.1 24.22 45.19 61.3 

10.0 5.5 8.25 -38.25 62.6 

10.0 5.8 -0.96 -21.83 65.6 

34.6 6.0 16.71 -61.47 64.2 

42.6 5.9 13.79 -91.95 79.7 

15.0 5.0 7.61 -35.81 69.8 

15.0 5.4 44.01 -28.97 55.7 

79.8 5.3 10.50 -62.56 72.0 

10.1 5.2 -1.13 -24.34 64.0 

10.1 5.2 -1.13 -24.34 84.2 

16.0 5.0 -1.03 -15.35 70.4 

15.0 5.3 46.25 -27.54 72.2 

10.1 6.0 38.13 20.38 67.5 

10.1 6.0 38.13 20.38 73.4 

10.1 6.0 38.13 20.38 74.2 

10.1 6.0 38.13 20.38 76.9 

10.1 6.0 38.13 20.38 79.6 

15.0 5.1 37.88 20.90 68.0 

15.0 5.1 37.88 20.90 80.1 

65.0 5.6 34.75 24.89 72.5 

65.0 5.6 34.75 24.89 72.5 

65.0 5.6 34.75 24.89 84.6 

25.0 5.3 38.58 20.54 79.5 

56.8 5.4 2.45 -76.81 84.0 

15.0 5.2 16.48 -46.45 59.0 

19.2 6.0 10.08 -62.61 71.9 

19 2 6.0 10.08 -62.61 80.5 

11.0 5.5 8.06 -82.87 76.3 

10.1 5.6 17.04 -46.52 58.6 

10.1 5.6 17.04 -46.52 68.4 

15.0 5.0 38.06 20.55 79.8 

15.0 5.2 18.32 -70.20 70.8 

25.0 5.2 36.51 21.37 69.2 

10.0 6.3 39.89 24.66 69.0 

10.0 6.3 39.89 24.66 812 

10.0 5.6 -0.66 -21.73 65.4 

10.0 5.6 -0.66 -21.73 84.9 

10.0 5.3 7.25 -34.35 56.9 

10.0 5.3 7.25 -34.35 65.3 

15.0 5.1 40.54 30.05 75.9 

11.6 6.1 40.47 42.05 82.4 

35.2 6.3 0.31 -79.99 79.0 


Bouncepoint 

Lat Long Azim 5t* grade 


34.43 

-43.80 

51.8 

0.8 A 

30.43 

-16.33 

39.3 

1.8 A 

12.07 

-36.84 

228.4 11.1 B 

11.76 

-36.72 

228.8 

8.7 A 

12.30 

-36.59 

228.4 

1.1 A 

35.97 

-41.46 

53.2 

3.1 C 

24.% 

40.75 

243.1 

0.8 C 

35.36 

-21.04 

29.8 

1.3 A 

23.05 

-43.72 

309.8 

3.7 B 

53.95 

-32.31 

265.4 

0.9 C 

36.31 

-41.50 

53.5 

0.5 A 

51.56 

-27.27 

33.5 

2.1 A 

36.09 

-22.51 

31.2 

2.1 C 

22.26 

45.63 

311.7 

0.0 C 

43.00 

40.02 

41.5 

4.5 C 

34.89 

-54.55 

67.5 

3.9 C 

25.46 

-67.28 

295.1 

3.8 B 

26.01 

-55.33 

223.2 

3.3 B 

40.00 

-39.90 

40.6 

3.9 C 

21.94 

47.10 

313.0 

1.9 C 

21.88 

-60.46 

298.3 

4.6 A 

25.27 

-39.54 

314.6 

8.5 C 

15.92 

-51.27 

208.3 

1.1 B 

53.20 

-23.77 

278.5 

3.9 A 

69.42 

-14.97 

310.1 

2.3 B 

60.68 

-27.05 

285.6 

4.3 A 

75.88 

7.28 

343.2 

9.5 A 

51.23 

-33.84 

267.5 

3.5 A 

53.22 

-23.37 

278.8 

4.7 A 

51.27 

-33.45 

267.8 

3.1 A 

51.90 

-19.67 

280.7 

2.3 C 

52.70 

-19.37 

282.0 

3.8 C 

51.02 

-29.53 

271.1 

5.6 C 

51.51 

-33.97 

267.0 

1.0 C 

37.69 

-52.33 

38.5 

3.5 C 

41.19 

-27.80 

38.7 

3.8 B 

35.38 

-34.55 

52.3 

1.4 B 

9.52 

-21.74 

94.4 

2.0 A 

29.78 

49.42 

62.3 

5.5 B 

41.49 

-27.80 

39.0 

7.5 A 

48.08 

-29.05 

30.8 

3.2 B 

51.25 

-33.72 

267.6 

6.1 A 

41.30 

-38.75 

59.0 

3.8 C 

52.65 

-22.49 

280.1 

2.7 C 

54.43 

-22.44 

275.3 

1.1 A 

53.52 

-32.82 

265.6 

3.8 B 

22.43 

45.58 

311.5 

2.3 B 

27.65 

-54.71 

306.1 

3.2 C 

30.35 

-16.63 

39.5 

0.0 B 

32.38 

-56.80 

315.3 

1.3 C 

55.59 

-23.73 

270.3 

2.0 C 

59.74 

-17.66 

270.4 

2.8 A 

24.75 

47.94 

56.8 

5.4 A 
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Origin time Epicenter Bouncepoint 


Stn 

Yr Mo 

E* 

Hr Mn Sec 

Dep 

mb 

Lat 

Long 

Dist 

Lat 

Long 

Azim 

St* 


RSSD 

83 

12 

22 

4:11:29.3 

10.0 

6.4 

12.02 

-13.64 

82.0 

36.79 

-50.04 

297.6 

3.3 

B 

ANMO 83 

12 

22 

4:11:29.3 

10.0 

6.4 

12.02 

-13.64 

85.5 

32.10 

-54.74 

288.4 

5.3 

C 

SCP 

84 

2 

11 

8: 2:50.0 

15.0 

5.3 

38.11 

21.86 

72.6 

51.92 

-26.67 

272.9 

4.0 

C 

TOL 

84 

2 

26 

8:18:16.6 

115.6 

5.9 

-17.31 

-70.65 

84.1 

13.14 

-41.69 

44.1 

2.8 

C 

KEV 

84 

3 

20 

17:19: 6.4 

10.0 

5.1 

10.75 

-41.96 

72.8 

44.28 

-25.70 

27.7 

4.2 

B 

BER 

84 

3 

30 

7:59:53.6 

10.0 

5.8 

17.65 

-59.84 

62.6 

43.30 

-38.72 

41.5 

1.3 

B 

GAC 

84 

4 

22 

6:14:21.7 

10.0 

5.8 

-0.14 

-19.59 

67.0 

25.32 

-42.14 

315.8 

1.0 

B 

GAC 

84 

5 

6 

9:12: 1.7 

10.0 

5.0 

38.66 

25.57 

70.2 

54.88 

-21.07 

277.9 

-0.3 

C 

GAC 

84 

6 

21 

10:43:40.5 

34.3 

5.7 

35.74 

23.80 

71.2 

52.93 

-20.70 

280.9 

6.2 

C 

RSCP 

84 

6 

21 

10:43:40.5 

34.3 

5.7 

35.74 

23.80 

83.3 

51.15 

-30.87 

269.9 

8.6 

C 

NE13 

84 

6 

24 

11:17:12.0 

16.4 

6.0 

18.09 

-69.23 

59.7 

32.61 

40.27 

68.0 

1.5 

B 

TOL 

84 

6 

24 

11:17:12.0 

16.4 

6.0 

18.09 

-69.23 

59.7 

33.26 

40.56 

66.7 

2.5 

A 

GRFO 

84 

6 

25 

18:45:10.3 

42.0 

5.1 

18.13 

-69.35 

70.3 

40.98 

-38.31 

58.7 

1.5 

B 

GAC 

84 

7 

19 

5:22:15.7 

10.0 

5.4 

-6.68 

-12.48 

76.5 

22.44 

-37.87 

316.4 

7.7 

A 

GDH 

84 

7 

19 

5:22:15.7 

10.0 

5.4 

-6.68 

-12.48 

80.9 

32.53 

-22.93 

343.8 

7.5 

B 

RSON 

84 

8 

30 20:12:59.7 

10.0 

5.0 

7.19 

-33.72 

65.7 

32.42 

-56.39 

315.1 

3.4 

A 

ANMO 84 

8 

30 20:12:59.7 

10.0 

5.0 

7.19 

-33.72 

71.8 

25.51 

-66.08 

294.9 

2.4 

B 

GAC 

84 

10 

21 

3:51:47.0 

10.0 

5.1 

-6.71 

-12.38 

76.6 

22.44 

-37.81 

316.4 

1.3 

B 

KONO 

84 

11 

1 

4:48:49.9 

10.0 

6.5 

8.36 

-38.82 

62.7 

36.20 

-22.90 

31.5 

1.7 

A 

ANMO 84 

11 

1 

4:48:49.9 

10.0 

6.5 

8.36 

-38.82 

66.9 

25.49 

-69.04 

295.3 

0.9 

C 

GAC 

84 

11 

5 

4:17:32.9 

10.0 

5.1 

-11.54 

-13.79 

79.4 

19.60 

-38.91 

319.3 

3.5 

A 

GAC 

84 

12 

5 

13:39:23.6 

10.0 

5.1 

-0.46 

-13.36 

71.2 

25.81 

-38.31 

313.5 

2.6 

C 

SCP 

84 

12 

8 

12:24:41.9 

10.0 

5.7 

-0.86 

-22.95 

64.7 

22.21 

46.29 

312.2 

3.5 

C 

RSON 

84 

12 

8 

12:24:41.9 

10.0 

5.7 

-0.86 

-22.95 

78.6 

29.44 

49.20 

315.6 

5.7 

B 

RSSD 

84 

12 

8 

12:24:41.9 

10.0 

5.7 

-0.86 

-22.95 

84.2 

27.32 

-55.49 

306.5 

6.7 

B 

ANMO 84 

12 

20 23:31:18.0 

10.0 

52 

6.94 

-33.32 

72.2 

25.42 

-65.83 

295.0 

2.0 

B 

KONO 

85 

1 

10 

17:47:56.1 

10.0 

5.8 

10.83 

-43.44 

62.6 

37.93 

-26.02 

34.6 

2.6 

A 

KEV 

85 

1 

10 

17:47:56.1 

10.0 

5.8 

10.83 

-43.44 

73.2 

44.50 

-26.91 

28.1 

3.7 

A 

NE16 

85 

3 

16 

14:54: 1.1 

14.0 

6.2 

17.10 

-62.34 

60.9 

35.86 

-35.35 

59.1 

8.5 

C 

NE06 

85 

3 

16 

14:54: 1.1 

14.0 

6.2 

17.10 

-62.34 

62.3 

38.29 

-36.30 

54.4 

2.9 

C 

GRFO 

85 

3 

16 

14:54: 1.1 

14.0 

6.2 

17.10 

-62.34 

66.5 

39.16 

-33.77 

56.6 

5.8 

B 

ANTO 

85 

3 

16 

14:54: 1.1 

14.0 

6.2 

17.10 

-62.34 

83.1 

38.60 

-21.61 

70.5 

7.3 

C 

GRFO 

85 

4 

20 

18:23:48.2 

12.3 

5.6 

;oi 

-77.35 

82.3 

37.59 

44.56 

54.5 

3.5 

A 

RSNY 

85 

6 

4 

12: 6: 3.6 

10.0 

5.1 

08 

-19.07 

66.0 

24.79 

41.76 

314.8 

2.4 

A 

ZOBO 

85 

6 

5 

1:41:42.3 

14.5 

5.1 

56.72 

-34.72 

77.9 

20.99 

-56.09 

198.6 

4.8 

C 

RSNY 

85 

6 

6 

2:40:12.8 

10.0 

6.3 

1.00 

-28.08 

59.7 

24.49 

47.18 

318.8 

4.0 

C 

TOL 

85 

6 

10 

3:23:31.8 

26.0 

5.6 

3.24 

-78.99 

76.4 

26.34 

47.28 

58.3 

1.4 

C 

BCAO 

85 

6 

26 

17:10: 1.9 

26.6 

5.6 

19.25 

-64.73 

82.2 

15.67 

-21.74 

101.4 

2.0 

C 

ANTO 

85 

6 

26 

17:10: 1.9 

26.6 

5.6 

19.25 

-64.73 

83.5 

40.50 

-22.72 

72.1 

6.1 

c 

GAC 

85 

7 

1 

9:53:36.1 

10.0 

5.0 

1.13 

-26.48 

61.7 

25.37 

46.34 

318.4 

5.0 

c 

RSON 

85 

7 

1 

9:53:36.1 

10.0 

5.0 

1.13 

-26.48 

74.9 

30.07 

-51.55 

315.7 

1.7 

c 

RSNY 

85 

7 

22 21:32:27.9 

15.0 

5.4 

34.16 

28.40 

75.0 

52.65 ‘ 

-17.70 

280.9 

4.3 

B 

RSON 

85 

9 

22 

18:23:12.2 

10.0 

5.7 

12.46 

-44.47 

55.2 

34.01 

-63.46 

316.6 

0.7 

B 

RSNY 

85 

9 

23 

17:28:41.8 

10.0 

5.2 

-18.00 

-13.75 

83.3 

15.24 

-39.34 

321.7 

5.0 

C 

RSNY 

85 

9 

27 

16:39:48.6 

43.8 

5.5 

34.05 

26.94 

74.2 

52.15 

-18.39 

281.1 

0.7 

B 

ANTO 

85 

10 

5 

15:24: 2.2 

10.0 

6.5 

62.53 - 

123.93 

76.2 

75.73 

4.81 

144.2 

2.8 

C 

SCP 

85 

10 

12 22:20:37.6 

10.0 

5.3 

0.99 

-29.15 

59.2 

22.69 

49.92 

314.0 

1.1 

C 

ANMO 85 

10 

12 22:20:37.6 

10.0 

5.3 

0.99 

-29.15 

79.0 

22.49 

-63.27 

298.0 

7.5 

C 

RSON 

85 

10 

27 

19:34:57.0 

10.0 

5.5 

36.46 

6.87 

68.7 

55.92 

-35.15 

286.7 

7.1 

A 

GAC 

85 

11 

10 

12:39:50.8 

10.0 

5.2 

4.57 

-31.96 

55.7 

26.73 

49.69 

319.4 

2.2 

B 

RSON 

85 

11 

10 

12:39:50.8 

10.0 

5.2 

4.57 

-31.% 

68.9 

31.23 

-55.18 

315.9 

6.8 

C 

RSSD 

85 

11 

10 

12:39:50.8 

10.0 

5.2 

4.57 

-31.% 

74.0 

29.06 

-61.25 

305.7 

3.6 

A 

KEV 

85 

11 

10 12:39:50.8 

10.0 

5.2 

4.57 

-31.% 

75.3 

40.02 

-17.80 

23.6 

1.2 

C 

GAC 

85 

11 

21 

21:57:14.9 

10.1 

5.4 

41.49 

19.69 

65.0 

54.67 

-25.68 

274.9 

2.6 

B 

RSNT 

85 

11 

21 

21:57:14.9 

10.1 

5.4 

41.49 

19.69 

70.1 

70.79 

-18.08 

306.6 

7.7 

C 
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Origin time 


Stn 

Yr Mo 1 

Da Hr Mn Sec 

ANTO 

85 

12 

23 19:37:54.8 

GDH 

85 

12 

24 4:13:21.2 

KONO 

86 

2 

12 23:41:38.5 

GAC 

86 

3 

6 0: 5:38.3 

RSON 

86 

3 

6 0: 5:38.3 

RSNT 

86 

5 

5 3:35:38.8 

NE16 

86 

5 

7 20:43:32.6 

NE16 

86 

5 

7 22:47:10.2 

NE10 

86 

5 

7 22:47:10.2 

RSON 

86 

5 

13 8:44: 2.0 

RSSD 

86 

5 

21 1:45:24.8 

GAC 

86 

5 

22 19:52:19.5 

RSNY 

86 

5 

22 19:52:19.5 

GRFO 

86 

6 

9 2:17:38.3 

NE13 

86 

6 

11 13:48:3.3 

TOL 

86 

6 

11 13:48:3.3 

NE11 

86 

6 

11 13:48: 3.3 

BCAO 

86 

6 

11 13:48: 3.3 

RSCP 

86 

6 

15 23:58:44.1 

RSON 

86 

6 

.7 8:48:18.6 

WFM 

86 

6 

24 6:56:53.0 

RSNY 

86 

6 

24 6:56:53.0 

GAC 

86 

6 

24 6:56:53.0 

RSCP 

86 

6 

24 6:56:53.0 

GAC 

86 

7 

10 6:53: 4.3 

RSON 

86 

7 

10 6:53: 4.3 

RSSD 

86 

7 

10 6:53: 4.3 

GRFO 

86 

7 

18 17:22:41.6 

GRFO 

86 

8 

7 22:32:50.9 

GAC 

86 

9 

13 17:24:33.7 

RSNT 

86 

9 

13 17:24:33.7 

ANMO 86 

9 

20 1:31:14.0 

RSNY 

86 

10 

2 10:12:39.8 

ZOBO 

86 

10 

27 0: 9:31.9 

KONO 

86 

10 

27 14:11:58.0 

GAC 

86 

12 

7 14:17:10.5 

ZOBO 

86 

12 

8 3: 3:26.1 

GRFO 

87 

1 

13 13:23:59.7 

GRFO 

87 

1 

13 19:30:10.9 

GDH 

87 

2 

1 6:56: 1.2 

NE14 

87 

3 

6 4:10:41.9 

GAC 

87 

3 

12 23:10:31.4 

GAC 

87 

4 

14 0:13:13.1 

SCP 

87 

5 

5 10:50:55.3 

GAC 

87 

5 

5 10:50:55.3 

GDH 

87 

5 

5 10:50:55.3 

GAC 

87 

7 

28 1:44: 9.7 


Epicenter 


Dep 

mb 

Lat Long 

Dist 

10.0 

5.5 

62.24 -123.95 

76.5 

10.0 

5.0 

7.53 -33.88 

63.0 

15.0 

5.2 

17.30 -62.41 

66.1 

35.0 

62 

40.33 51.60 

82.1 

35.0 

6.2 

40.33 51.60 

84.2 

15.0 

5.8 

37.72 37.70 

77.7 

23.1 

6.1 

51.60-174.65 

83.0 

31.3 

6.5 

51.33 -175.43 

83.3 

31.3 

6.5 

51.33 -175.43 

85.9 

15.0 

5.6 

41.03 43.92 

81.3 

15.0 

5.2 

14.90 -20.02 

75.5 

33.2 

5.1 

34.12 26.72 

74.0 

33.2 

5.1 

34.12 26.72 

74.1 

49.5 

5.0 

54.59 -168.10 

76.1 

20.0 

6.0 

10.70 -62.93 

59.1 

20.0 

6.0 

10.70 -62.93 

59.4 

20.0 

6.0 

10.70 -62.93 

61.7 

20.0 

6.0 

10.70 -62.93 

80.8 

15.0 

5.0 

1.16 -26.46 

64.5 

15.0 

5.4 

4.55 -32.55 

68.5 

15.0 

5.7 

0.53 -17.41 

63.9 

15.0 

5.7 

0.53 -17.41 

66.8 

15.0 

5.7 

0.53 -17.41 

67.8 

15.0 

5.7 

0.53 -17.41 

72.0 

15.0 

5.4 

4.31 -32.49 

55.6 

15.0 

5.4 

4.31 -32.49 

68.8 

15.0 

5.4 

4.31 -32.49 

73.8 

15.0 

5.8 

10.84 -69.24 

75.7 

15.0 

5.4 

7.53 -81.06 

85.8 

15.0 

5.8 

36.80 22.64 

69.8 

15.0 

5.8 

36.80 22.64 

75.3 

15.0 

5.4 

0.72 -29.00 

79.3 

15.0 

5.3 

34.65 29.16 

75.2 

15.0 

5.3 

46.06 -27.27 

72.2 

15.0 

5.2 

7.26 -36.20 

62.7 

15.0 

5.1 

43.01 26.01 

61.5 

15.0 

5.1 

47.72 -27.70 

73.2 

15.0 

5.5 

5.68 -78.90 

85.8 

15.0 

5.3 

5.87 -78.86 

85.6 

15.0 

5.6 

0.55 -17.47 

72.7 

15.0 

6.5 

-0.06 -77.84 

77.5 

15.0 

5.4 

0.36 -17.81 

67.7 

15.0 

5.3 

-14.96 -14.74 

81.5 

15.0 

5.0 

-0.11 -18.78 

67.1 

15.0 

5.0 

-0.11 -18.78 

67.5 

15.0 

5.0 

-0.11 -18.78 

73.1 

15.0 

5.0 

1.02 -25.62 

62.3 


Bounce point 


Lat 

Long 

Azim 

5t* 


75.81 

4.32 

143.6 

9.4 

B 

38.65 

-39.01 

350.2 

2.3 

A 

43.78 

-39.03 

44.2 

1.6 

A 

64.38 

-6.78 

277.1 

5.3 

B 

73.02 

4.11 

289.4 

6.6 

B 

73.68 

8.39 

321.7 

2.4 

C 

87.01 

-15.01 

160.9 

2.5 

A 

87.69 

-12.41 

163.8 

3.4 

B 

85.29 

-29.17 

149.2 

6.2 

C 

70.26 

-11.96 

285.8 

4.3 

A 

37.05 

-54.45 

296.6 

6.3 

B 

52.94 

-18.29 

282.4 

0.9 

B 

52.13 

-18.61 

281.0 

0.5 

A 

88.50 

1.47 

169.9 

1.4 

B 

27.77 

-37.23 

59.9 

6.2 

B 

28.41 

-37.48 

58.6 

0.6 

B 

29.75 

-36.91 

57.4 

3.2 

C 

9.94 

-21.85 

94.8 

2.8 

B 

20.87 

-52.67 

304.1 

4.4 

C 

31.15 

-55.56 

316.1 

1.4 

C 

23.86 

-40.00 

313.4 

3.4 

C 

25.19 

-40.74 

313.8 

4.4 

B 

25.91 

-40.81 

314.5 

6.9 

B 

21.44 

47.49 

301.4 

10.2 

A 

26.56 

-50.00 

319.9 

6.3 

A 

31.03 

-55.52 

316.2 

7.1 

A 

28.83 

-61.57 

306.0 

5.9 

B 

36.97 

-38.90 

54.2 

1.2 

B 

37.55 

47.27 

54.3 

1.0 

C 

53.15 

-21.91 

279.8 

3.9 

C 

69.40 

-11.53 

312.3 

3.9 

B 

22.34 

-63.18 

298.2 

1.6 

B 

53.13 

-17.46 

280.5 

7.1 

B 

15.83 

-51.12 

208.5 

3.2 

C 

35.38 

-21.13 

29.8 

1.5 

B 

57.08 

-23.10 

273.2 

3.9 

A 

16.68 

-51.61 

207.3 

1.4 

C 


35.99 -45.85 53.0 5.0 B 

36.10 -45.81 53.1 1.0 B 

35.92 -26.67 344.3 1.3 B 

22.76 -45.63 58.5 2.2 B 

25.77 -41.04 314.8 6.4 C 

17.59 -39.72 321.2 3.6 C 

23.01 -43.84 309.9 3.6 A 

25.42 -41.65 315.5 4.3 A 

35.50 -27.66 344.9 5.0 A 

25.39 -45.83 318.0 0.6 A 


gt* — Differential attenuation between the S and SS phases 
Other abbreviations as in Appendix 1 


Appendix III. SS-S Travel Time Data for the East Pacific Region 
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